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ABSTRACT 


The controlled behavior of a fully instrumented 
eight tray, nine inch diameter glass distillation column, 
separating a mixture of methanol and water, has been 
studied using feedback, feedforward and combined feed- 
forward-feedback control. Single point overhead product 
composition control was implemented by manipulation of 
the reflux flow rate compensating for disturbances in the 
feed flow rate. The dynamics of the column, determined 
using the pulse testing techniques, were successfully 
represented by first order plus time delay models, while 
the process gains were determined using the transient 
response testing procedure. Due to the nonlinear process 
responses exhibited by the column, two process models 
were determined, one for increasing and the other for 


decreasing overhead product composition responses. 


The feedforward operating characteristics, measured 
experimentally, were found to be linear, resulting in a 
constant feedforward gain over the range of feed flow 
rates studied. The feedforward gain was subsequently 
found to exhibit a slight functional relationship with 
the feed flow rate during the on-line evaluation of the 


feedforward controllers. 
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Calculation of the frequency response of the feed- 
forward controller using the frequency responses deter- 
mined from the open loop pulse tests showed that the 
process responses could be characterized by simple trans- 
fer function models. However, the feedforward controller 
frequency response values calculated from the simple models 
used to predict the transient response did not agree with 


those calculated directly from the pulse tests. 


Gain, gain plus time delay and gain plus time lag 
feedforward controllers; gain and gain plus time lag 
feedforward controllers with feedback trim were implemented 
using an IBM-1800 digital computer. The effectiveness of 
these control schemes was compared to that obtained using 
a conventional feedback controller for feed flow rate dis- 
turbances. Both increases and decreases away from and 


returning to a reference steady state were employed. 


ine vesults indiceterthat:: 

(a) gain feedforward control is not as effective 
as feedback control 

(b) dynamic feedforward control in the form of 
either a time delay or first order time lag 
results in an improved control behavior 
compared with that obtained using only 


feedback control | 


-bast ait to cellent super | 
Pwr Ce. 


tated esenoq2e7 yonouper? aitd pnteusy 
gag ted’ baworle otead ‘setug goal | 


il 
-ensys slamre yd bos resoesda ow 09 2 


ieffordnos byawio tees’ one snevewall 
ofsbom atamte eat mort botstuolso 20 ay avg 


djfw se vps ton brb Seog 2o1 ascot sai ¥ 


pst ‘emtt 2ulq wfsp bas ysiteb emt evtg atop i 
cont 4 

est emtt evfq atsp brs ntsp vat Toxdnoo bam 
wats 


batnemafqmt avow mive Aosdbset nitTw pra f'Tortnos bre ed 
. 


to 22onevittoottes eft .vetuqmos Tathpro 008f-MaT ns Ba 
. { al ; J 
pnfew bentssdo Fens oF boveqmos ote ay Oe 


atb offs wol? beat vo? velfoxrtaod soadbast reno tia 08 


bis mov? yews 25289%99b bis esessvont Asod > vepanse 
RYT yy 

.bsyolqma 979w etete vbsstea sons 1etoy 5 of anti 

~ y 


ke a 


7 


Font sistent es Tuest @ 

svitostts 26 tow @f Tovtnoo byéswrotbsst nye 
ford nos josdvest 26) ©) | 

to myot add af fowsnoo ial ah stash 
ce, HA 


| | ps ale 
osf amts tsbv0o tevt? vo ysteb omts phinishish- 
7 


re 


rotveded Tovtmos bevorgatinn nb ad tues 

vino pntey bentssdo todd date boveqnes 

: | Pires, tovtnon Aaed bes) 
. oo | 


the addition of feedback trim to a feed- 
forward controller, besides eliminating 
offsets in the controlled peenonee due to 
errors in the gain determination or im- 
measurable disturbances, also improved the 
CF ec tiveness Of sene Combine dy.comtro | 
scheme over that obtained using only the 


same feedforward controller. 
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CHAPTER 1 
INTRODUCTION 


Distillation columns are one of petrochemical 
industry's most versatile and important unit operations. 
There are relatively few processes available which do not 
involve at least one distillation column. For a large 
number of these processes the distillation column is at 
the heart of the process, where the ability to maintain 
close control over the distillation column reflects the 
Obie yercOn ma lntannie economic: \ con trol over; that) proces’s. 
This is exemplified in such processes as catalytic crack- 
ing, where even small increases in the recovered gasoline 
yield and the ability to maintain it can markedly increase 


ChempVordcd blac yyor the unit. 


Despite the extensive use of distillation columns in 
industry, they remain among the most difficult unit opera- 
tions to. control; due to the long characteristic  time..con- 
stants associated with the response of the column outputs. 
Since distillation columns are also generally among the 
final processing stages of a product, they, too, are sub- 
ject to all the disturbances which may occur throughout the 


process. 


The quality of the products produced from distillation 


a 
> 


vena Petry 
j saat ‘ a ci ‘ % at pr 
My iy i ‘eran | 


[so tnistsovtsg ‘yo 9m ons 2 
.2nofesersqo thru Snedrodet bei 2 
ton ob dotdw sfdsfteva esees ‘ vat 

spel 6 YOU -Amuloo norset raat 2 ‘b 
ts at nmulos norte rl rsetb eat 2 | : 
nistatsm of ydt hide. ‘ott ovale a 9d 0 S9604 
ot etoaTtey nmwfhos corres 9A weve” Tortmoo § 
.e2gso0%q tent yovo forioos atmonase nteta hsm ot ving 
-19612 Ditytledeo 26 292299019 © doue wt bot ti fqnaxs et 2 rw 
entlozsp bevavoosy odd at 252n979NF Tit ome nove o+oriw 
s269vont ylbsivem mo at atesarem os ystitrds ont bas is 


tm 9d Fo aie an 


nt 2amuloo hottehttseth to 920 aviensdxs ne 
-s1ego stau tIyortith Feom sis pnoiis ntamay yodt nese 
-oo 9mts ottetvetos iano paol ans of sub ,forvsnoo one 
.etuqivo amufoo sat To senoqesy 9dd ddttw. vosetaoees : 
odd proms yl fsvensp oats one enmufos notsetitaetb 
-dve 918 ,00t .ysdt .doubovg 5 to espede | entee2soo A 
odd juatpyorns t0990' Vem, yi wisi boas a wit % 


TOTP ; wt 


ps 
eo) o : , : 
al y é : | 
n 4 a 
ie Di al te a i) 
; me} a ey 


nottstitteth mort beow! 


~ 


columns is normally controlled using feedback controllers, 
which operate either on a direct measurement of the product 
quality or on an intermediate measurement, such as tem- 
perature, from which the product quality may be inferred. 
Many methods have been developed to select which inter- 
mediate variable gives the best control. However, the 
control of this intermediate variable does not generally 
guarantee that the product quality will remain constant, 
unless the exact relationship between the inferred and 
controlled variables is known exactly. For example, in 

a number of processes, temperature is used as an inferen- 
tial variable supposed to reflect the composition. How- 
ever, since temperature will only be an exact reflection 
of the composition for a binary system at constant pres- 
sure, control of a multicomponent product using a tem- 
perature to infer composition will not guarantee a con- 
stant composition product. The best measurement to use 

in a feedback control scheme is the measured value of the 
property of the product stream which it is desired to 


maintain constant. 


It is impossible to completely eliminate a transient 
deviation in the controlled variable using only a feedback 
controller because, as the name implies, a deviation from 
the set point must occur before any corrective action can 


be introduced. The amount of the deviation that will result 
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with a feedback controller, as measured by the integral of 
the absolute error (IAE), is proportional to the product of 
the proportional © band, the reset time and the disturbance 
size (115). “This product is generally large for most in- 
dustrial columns due to the long time constants associated 


with these columns. 


The feedback controller does not know, a priori, what 
control action to implement for a disturbance, but contin- 
ually adjusts its manipulative variable until the measured 
value of the controlled variable and its set point are in 
agreement. In essence, the feedback controller eliminates 
the effect of a disturbance using a trial and error solu- 
tion, characterized by the oscillatory behavior which 


often results. 


The advantages of using feedback control are that 
essentially no knowledge of the process dynamics is required 
and that corrective action can be implemented no matter 
which input variable caused the disturbance. The major 
disadvantages include the slowness of response, the os- 
cillatory control action and the large integral of the 


absolute error (IAE) obtained. 


From the preceding discussion, it would seem obvious 


that if the corrective action could be introduced when the 
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disturbance entered the column rather than after the dis- 
turbance has made itself evident in the product quality, a 
much more satisfactory control system would result. This 
is precisely the purpose of a feedforward controller. 
Feedforward control can, theoretically at least, completely 
eliminate the effect of a disturbance on the product qua- 
lity, if the effect of the disturbance on the product qua- 
hity pcan .besdescribed exactly... Since, in most cases, 
these relationships are very difficult to model, very 
Simple: funetions,esuch as ayfirsteorders lags(/9)4.a-first 
order lead-lag (115), a second order lag plus dead time 
(66) and a first order lead plus a second order lag (82) 
have been used to approximate these relationships. The 
most stringent constraint on these model descriptions is 
placed on the steady state, since small inaccuracies in 
the feedforward gain can be amplified into quite sizable 
offsets in the controlled variable due to the relatively 
large gains associated with the open loop responses. The 
values of the parameters involved in these approximate 
controller models are chosen to minimize the deviation 
(IAE) in the product quality. If the parameters of the 
approximate models are not chosen carefully, the controller 
may exhibit too much lead or lag action, resulting in a 
controlled response that would be less satisfactory than 


that obtainable using only a feedback controller. 


The main advantage of feedforward control is its speed 
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of response, with the lead action being obtained by intro- 
ducingpthe icorrectivewaction, onthe bas Tsof the measured 
disturbances entering the column rather than on the dis- 
turbances leaving the column. The major drawback of feed- 
forward control is that it has no indication whether or 
not it has taken the correct control action, making the 
effectiveness of the control dependant upon the accuracy 
of the feedforward model and the absence of unmeasured 


disturbances. 


Since their advantages seem to complement each other, 
it is obvious that the feedforward and feedback schemes 
Should be combined. The combined control scheme should 
exhibit the fast initial corrective action due to the lead 
action of the feedforward controller, and also the long- 
term removal of any offsets by feedback action no matter 
what their source. The presence of the feedforward control- 
ler does not make the process less difficult to control 
using feedback control; the gain in effectiveness is ob- 
tained because the feedback controller must. do less work 
to incrementally trim the corrective action taken by the 
feedforward controller. In order to eliminate any offset, 
integral feedback control is the minimum necessary. 
Whether proportional or derivative control are beneficial 
depends upon the need to correct for only offsets due to 


errors in the model description or for deviations caused 
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by immeasurable disturbances. 


For a simple distillation column, the main measurable 
disturbances include the feed flow rate, the feed composi- 
tion, the feed enthalpy and the reflux enthalpy. Special- 
purpose controllers have been proposed to control both of 
the latter (57). The disturbances which are generally not 
measured include the steam enthalpy, the cooling water 
enthalpy and the heat loss to the atmosphere. Both the 
reflux flow rate and the steam flow rate are considered 
as the major Rue ae ate variables, although, depending 
upon the control scheme chosen, the overhead and bottoms 
flow rate and the feed plate location have also been con- 
sidered. The overhead and bottoms product compositions 


are generally considered as the main controlled variables. 


The majority of the published industrial effort (57, 
66,68,82,90,115) appears to have been directed to develop- 
ing feedforward controllers. which compensate for feed 
flow rate disturbances, while most of the academic atten- 
Pion 13.30. 56,..6 9.7.0,./1.579,105.133) has been focused. on 
developing feedforward controllers. compensating for feed 
composition disturbances. Feed flow rates tend to vary 
much faster and more often than does the feed composition 
and, an fact, both. Shinskey, (115), and Lupfer, (57:68). suggest 


that feedback control designed on the basis of feed flow 
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rate disturbances should be sufficient to compensate for 
the gradual feed composition changes. They do stress, 
however, that no standard solution can be applied to each 
situation, resulting in an individual analysis and probably 


a unique solution for each case. 


Although industry seems to have accepted feedforward 
control as a standard control algorithm, surprisingly few 
results have been published. While Foxboro (34) has re- 
ported the installation of over 100 feedforward control 
Systems, of which more than half have been applied to dis- 
tillation columns, Scrimgeour (113) has recently reported 
that of over 150 feedforward control systems installed on 
distillation columns throughout Canada and the United 
States, only four applications have been described in the 
literature of the past eight years. The economic results 
reported (66,68,82,90,115) indicate that the gains possible 
in areas such as closer composition control, increased 
throughput or decreased utility usage can pay for the in- 
stallation of a feedforward controller in an extremely 


short period of time. 


The present study was initiated in order to obtain 
an experimental comparison of the effectiveness of a 
single-point feedforward controller with that obtained 


using a standard feedback controller subjected to distur- 
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bances in the feed flow rate. This information would hope- 
fully fill an apparent void in the literature. It was also 
hoped that a simple method for estimating the approximate 

feedforward controller parameters could be established, as- 
suming only a minimal knowledge of the dynamic behavior of 


the column. 


Although this study was not carried out on the com- 
plex multicomponent, multiplate, multiproduct columns 
characteristic of industrial columns, it was hoped that the 
principles outlined would help bridge the gap between 


theory and application. 
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CHAPTER 2 
LITERATURE SURVEY 


2.10flatreducétion 


The literature survey considered pertinent to this 
study will be presented under the following topics: 
i) distillation dynamic models 
ii) feedback control 
iii) feedforward control 
iv) pulse testing techniques 


v) complex curve fitting techniques. 


2.2 Distillation Dynamic Models 


The various methods which have been employed to re- 
present the complex relationships involved in distillation 
column dynamics will be reviewed. Williams (137) and 
Rademaker and Rijnsdorp (99) have presented excellent des- 
criptions of the complete set of difference-differential 
equations necessary to describe plate dynamics. These 
relationships result from the combination of equations de- 
scribing the material and energy balances, the phase equili- 
brium, the plate efficiency, the plate hydraulics, the 
plate mixing and the physical properties. These authors 


discuss the validity of some of the assumptions which are 
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commonly used in order to simplify the set of equations 
so that they can be solved more readily. In addition, 
Gould (37) and Rijnsdorp (100) emphasize the importance 
of some of the secondary effects such as hydraulic de- 
lays and mixing on the ultimate dynamic response. Other 
review articles available include those of Archer and 
Rothfus (2), Rosenbrock (107), and the annual reviews of 
distillation literature published in Industrial and 


Engineering Chemistry (51). 


Examination of the voluminous literature pertaining 
to the dynamics of distillation columns indicates that 
the information may be classified according to the assump- 
tions made and the method of solution utilized. These 


topics include 


(a) analytical solution of the linearized 
material balance equations 

(b) analog computer solution of the linearized 
material balance equations 

(c) frequency response stepping solution of the 
linearized material balance equations using 
a digital computer 

(d) numerical integration of the nonlinear material 
and energy balance equations using a digital 


computer 
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(e) approximate models where the parameters 
are determined from the steady-state 
conditions 

(f) approximate models where the parameters are 
determined experimentally 


(g) miscellaneous methods. 


The articles of Armstrong et al (3,4,5,140,143) 
illustrate the early attempts to obtain analytical solu- 
tions to the linearized material balance equations. A 
number of simplifying assumptions, including the exis- 
tence of a linear vapour-liquid equilibrium relation, 
must be made in order to obtain the analytical solution. 
These authors have compared the results of their analyses 
to experimental transients with moderate success. The 
largest deviations appeared in the portion of the re- 


Sponse for short times. 


Lamb et al (60,105) and Gerster et al (6,7,81,129) 
have presented a method of solving the linearized material 
balance equations developed at the University of Delaware 
using an analog computer. This procedure allows the use 
of a vapour-liquid equilibrium relation which has been 
linearized for each tray. The major disadvantage of this 
approach is the large analog computer which is generally 


required to solve even relatively small problems. The 
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results predicted by this model compare quite favourably 
with some experimentally determined transients. Lamb and 
Rippen (105) have also described a method of solution of 
the same set of equations in the frequency domain using 

a stepping procedure on a digital computer. These cal- 
culations result in a series of Bode diagrams which il- 
lustrate the relationships between the input and output 
variables. The frequency response can then be fit using 
a transfer function of the appropriate degree of complex- 
ity or simplicity. A more recent description of this 


procedure has also been presented by Bollinger (11). 


This model has been used extensively to describe 
the plant transfer matrix describing the response of a 
distillation column in subsequent control systems ana- 
lyses by Luyben (69,70,71,72,73,80,123), Janis (56) and 
Wardle and Wood (132) among others. Wood (141) and 
Cadman et al (20,21,22) have also adapted this model to 
describe the dynamic relationships involved in multicom- 


ponent distillation. 


Huckaba and co-workers (30,48,49) have presented a 
detailed theoretical and experimental evaluation of a 
nonlinear model. This procedure uses a specially de- 
veloped predictor-corrector numerical integration algo- 


rithm to solve the set of differential equations, in- 
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cluding the nonlinear material and energy balances, and 
nonlinear phase equilibrium and enthalpy relations. 
Similar models have been used in other dynamic studies 
by Svrcek (128), Davison (28) and Peiser and Grover (39, 
95). These latter workers have included a description 
of the plate hydraulics, plate holdups, reboiler and feed 
preheater dynamics to analyze an industrial column which 
was experiencing unstable bottoms product control. The 
Suspected cause of the problem, namely flooding, could 
not be established using a steady state model, but was 
subsequently verified by solution of the dynamic model. 
The modifications proposed, evaluated using the dynamic 
model, improved the situation and resulted in stable 


bottoms product contro]. 


The calculation of the dynamic responses of dis- 
tillation columns using the models outlined previously, 
generally requires a large general-purpose digital com- 
puter. In many cases, this approach to control systems 
cannot be justified nor can the large general-purpose 
computers be used for real-time on-line process control. 
It seems obvious then, that methods of either predicting 
or determining approximate models for the column is very 
necessary. A number of workers have attempted to corre- 
late the parameters of simplified transfer functions based 


on responses calculated from linearized models similar to 


ae). sei aatow bas: (08): 
notsqtroeeb 5 babulont. svar: sd 

best bas t9lhodex ,2qublod steiq. nent Ble 
do tew naw 10.2: istyteubat as sehen (20% 

ant forties toubord 2mos3 od aided enn onto 
biwoo »pntboolt yfemen. emaldord ond to, 32 
26w Jud ,.febom stete ybsote 6 ey, beniar 
-Tobom atmsayh edd to notsuloe, vd bet Titey ar 

stmenyb ott pateyv bodsulsve .bseogorq anordeott bom 


“~ iE : 


ofdste at betiugen bas nora suate: ont bovorqmh fetem 
Tondwon Souborg enosted 


; oa : 
a 


@ %y ; . Lb : 
-2fb to egsenogesy afmsanyb sit to nottetuabas edt 7 
-vlewotverq bent fiuo efsbom add enfew cio ma 7 
-moo fettprb seoqrug-Isvense Spi1st 6 2s tupsy Vi horse § | 


amesey2 fortnos 03 dosotggs eid, 120269 sit pore 
szoqruq-Isionep spre! eft 59 lianas; 
.fovsno2 22920%q antl-no smis-I597 ‘0% Bess me 
enttotberq isAldts to ae | 
v79v ef nmutoo edt vor oer ixe Pee (en en en i 
-oy102 ot betqmotts onsale noivow to yedmum A .yyezeeost . 
bozsd enotdonut etensyd bektt lambs to exesemnmed | sal ew 
Pe PEA hie i 


- 


C1 Ti imental 
mane Cael: 
Fr. | 


r 7 f 


ety 


ge 


ib ; i . is i 
\ Ay : m - an 
7 Q [on 


iii 
7 
ol 
; 
) 


that presented by Gerster et al, generally with limited 
success. These studies include the work performed by 
Gilliland and Mohr (36), Mohr (87), Moczek, Otto and 
Williams (85,86), Bhat and Williams (10) and Wahl (133, 
134). Wahl's approximate model is based on the solu- 
tion of. the model of Genstermetual.cusing.ithecstepping 
procedure of Rippen et al for a large number of initial 
steady states. Examination of all the Bode diagrams 
indicated that the column responses could be represen- 
ted by simple transfer functions of the ote 
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d) reflux flow 


for all plates X K 
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The time constants (T)sT3>T7»Ty) calculated from these 
Simple transfer functions have been presented as graphi- 
cal correlations based on the initial steady state condi- 
tions. Both Wahl (133) and Beaverstock (8) have used 
this method to model a distillation column for an ana- 
lysis of a feedback control system. Williams et al (10, 
85,86) assumed the response of the column could be repre- 
sented by a second order plus time delay transfer function. 
They have developed a relationship between these time 
constants and an overall time constant, called the inven- 
tory time, tor the distillation column with its plate 
holdups lumped into an effective holdup. This method was 
used to analyze the control of a large industrial column 


producing very pure products (86). 


Osborne et al (92) have proposed dividing the columns 
into sections and writing a set of linearized equations 


describing the material balance by lumping all the 
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individual tray holdups in each section into a single 
overall value. The vapour and liquid material balance 
equations for each section are coupled by the mass trans- 
fer occurring between the phases. This mass transfer 
constant can be determined from the steady state rela- 
tionships while the holdup parameter must be adjusted 

to fit the experimental response. By adjusting parame- 
ters, this method could predict experimentally determined 


responses as well as any other. 


Approximate models determined using various standard 
dynamic testing procedures have been evaluated for a 
number of different industrial columns (50,98,112,114, 
136,142). The models developed have been presented as 
Bode diagrams with no particular attempts made to fit 
these responses to an approximate transfer function. 
Rademaker (98) subjected a tall turbogrid column to both 
frequency and transient response tests to determine the 
response of the column to various input disturbances. 
These measured responses compared favourably with the 
responses predicted from a model outlined previously by 
Rademaker and Rijnsdorp (99). Wood and Robins (142) 
have compared the response of a 30 tray industrial co- 
lumn measured using step, sinusoidal and stochastic 
variations in various input variables to the responses 


predicted using the linearized perturbation model of 
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Gerster et al. The comparison was considered sufficiently 
successful to justify its use in a subsequent control 
study (132). Janis (56) also arrived at a similar con- 
clusion by comparing the responses of a distillation 
column determined using the pulse testing technique with 


those predicted from the linearized model of Gerster et al. 


Numerous miscellaneous methods have been proposed 
to determine the dynamic response of a distillation co- 
lumn. These include, among others, solution of the 1li- 
nearized material balance equations using flow graph 


analysis (52) and matrix algebra (64). 


Levy et al (64) have presented an interesting com- 
parison of the major time constants of three different 
model formulations based on an eigenvalue analysis. These 
models consisted of equations established from 

a) component, mass and enthalpy balances 

b) component and mass balances 


c) component balance. 


It is interesting to note that for all cases the 
dominant time constant is associated with the accumula- 
tion of chemical species which appears to be equally 
disturbed throughout the column. The second largest 


time constant of the more comprehensive model (a) indi- 
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cates that the mode is associated with the thermal be- 
havior of the accumulator. The temperature dependence of 
this mode disappears as the condenser and reboiler hold- 
ups are made smaller. This analysis indicates under what 
type of conditions the solution of the linearized compo- 
nent balance equation might fail to provide an adequate 


representation of the response of the column. 


Although most of the column descriptions developed 
to date have concerned binary systems, a few articles 
have been presented for the multicomponent case (19,111, 


141). 


Examination of the results presented previously 
indicate that the numerical methods are advanced enough 
to solve the problem for either the linear or nonlinear 
case. Although the set of equations describing the dyna- 
mics can be easily developed, their solution is made 
difficult by the lack of sufficient knowledge concerning 
the evaluation of the secondary effects, such as dynamic 
plate efficiencies, plate hydraulics and plate mixing. 
Under most conditions the contributions of these effects 
are small, but can become important limiting conditions. 
The majority of the studies reported here involved the 
solution of linear models. However, a great number of 


industrial columns, ie product purification columns, 
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operate in a nonlinear range where extremely pure products 
are produced. More work should be done to determine the 
effects of these nonlinearities. The large models are 

not generally useful for control system analysis due to 
the substantial amount of computer time required to 

obtain a solution. More work should be done to provide 
workable methods to predict approximate dynamic relations 
for the column, which could be used in either the ana- 
lysis of, or as models in, an on-line computer control 


system. 


2.3 Feedback Control 


Since the theory of feedback control has been fully 
developed in most standard process control texts (16,27, 
39,40,115), only those references of specific interest 


to distillation column control will be mentioned. 


Numerous authors (44,74,101,108,109,138,139) have 
reviewed the various available methods which are commonly 
used to maintain control of the material and energy balances 
in the column. It is generally accepted that the product 
quality is much more sensitive to changes in the material 
balance rather than changes in the energy balance. Two 
methods are available to control the material balance. The 


first method, designated as direct, involves controlling 
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the overhead product flow rate and allowing the reflux 
flow rate to be reset by the level in the reflux accumu- 
lator. .The second method, designated as indirect, in- 
volves controlling the reflux flow rate and allowing the 
level in the reflux accumulator to reset the overhead 
puoduct flow rate: Both Shinskey (115.116,117,118), and 
Nisenfeld (89,90) have advocated the direct approach, whiie 
Buckley (17) seems to prefer the indirect approach. Using 
both a steady state and dynamic simulation, McCune and 
Gallier (84) have presented a comparison of both approa- 
ches to the material balance control. They concluded 
that the direct approach would be the most preferred con- 
figuration. Despite these conflicts, the literature 
contains many references illustrating the popularity of 
both the direct and indirect material balance control 


schemes. 


Lupfer et al (67) have outlined the development of 
various special-purpose controllers to control the in- 
ternal reflux, feed enthalpy’ and the reboiler heat and 


bottom level. 


If a direct analysis of the overhead or bottoms compo- 
sition cannot be made, then control may still be obtained 
using any other measurable variable to which the controlled 


variables can be related. The control of product composi- 
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tion by control of an intermediate tray temperature has 
been analyzed by Jafri et al (55), Luyben 7 5e 6) 20N)k. 
Shoneman and Gerster (119), Harriott (39) and Chanh (23). 
Chanh has compared the predicted controlled response of 

an eight tray pilot distillation column with that measured 


experimentally using various control plates. 


The simultaneous control of both overhead and bot- 
toms product compositions using only feedback control is 
Complicated by the interactions experienced between the 
two manipulative variables, usually reflux and steam flows. 
Rijnsdorp (102,103,104), Luyben (73) and Berry (9) have 
illustrated various methods available, which are capable 
of decoupling the interaction, making the two-point feed- 
back control scheme feasible. A recent article by Masel- 
li and Miller (83) suggests that the interactions existing 
in a two-point feedback control scheme can be minimized by 
utilizing a control scheme at one end of the column, con- 
sisting of an intermediate plate temperature controller 
whose Cascade set point is reset from the overhead compo- 


sition controller. 
2.4 Feedforward Control 


Feedforward control has been successfully applied 


to numerous unit operations other than distillation columns, 
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IicMudIng reactors (425777807131) e driers ae,’ neutral 
izers (34), boilers (34) and evaporators (34,53,88). 
This section, however, will briefly describe only a few 


of the applications related to distillation column control. 


The general theory concerning feedforward control has 
been described in a number of recent control texts (38, 
45,79,115). Shinskey (116,117) has derived a simplified 
feedforward controller, based on a direct material 
balance control scheme, which maintains the ratios of D/F 
and S/F constant for disturbances in the feed flow rate 
and composition. Adequate dynamic compensation was oe 
tained using a first order lead/second order lag dynamic 
compensator. This model was used by MacMullen and 
Shinskey (82) and Nisenfeld and Stravinski (90) as the 
basis for developing feedforward controllers for a large 
industrial superfractionator and an industrial azeotropic 


distillation column. 


Lupfer et al (57,58,65,66,68) developed a simple 
feedforward controller which was used to control two 
industrial butane splitters. The model calculates the 
bottom product flow as a function of the feed flow and 
composition from the overall material balance. Steady 
state simulations were performed to obtain a relation 
between the reflux flow rate and the feed flow rate and 


composition. The reflux was then incorporated into the 
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model using a regression equation which provided the 
best fit to the simulation results. A dynamic compensa- 
tor consisting of a second order lag plus a time delay 
was added to the feed flow. Since the changes in feed 
composition generally occur with a low frequency, no 
dynamic compensation was considered necessary for a feed 
composition disturbance. Bornard et al (14) have applied 
a first order lead/lag feedforward controller to a large 
industrial superfractionator based on a model determined 
experimentally using a transient response analysis. 
Other industrial applications of feedforward control to 
fractionators have been presented by Skillern and Wil- 
liams (121), Roach (108), Svrcek and Wilson (129) and 
Maselli and Miller (83). 


Rippen and Lamb (105) have utilized their stepping 
procedure in conjunction with the linearized model of 
Gerster et al to determine the two-point feedforward 
controllers required to maintain both overhead and bot- 
toms product compositions constant by manipulation of the 
reflux and steam flows to correct for disturbances in the 
feed flow rate and the feed composition. The feedforward 
controllers, presented as Bode diagrams, were fit with. 
combinations of first and second order lag terms. The 
control of the distillation column, described by Gerster's 


model, was simulated on an analog computer using the appro- 
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ximate feedforward controllers. King et al (59) have 
obtained an expression for the steady state gain based 
on the solution of the algebraic equations of Rippen's 


stepping procedure. 


Bollinger and Lamb (12) have formalized the method 
required to synthesize an ideal feedforward controller 
given the plant transfer matrix. This procedure was 
later extended to include feedback trim of the feedforward 


controller (13). 


Luyben and his co-workers (69,70,71,72,73,/78,/9,124) 
have made extensive use of the procedures presented by 
Rippen and Lamb and Bollinger and Lamb throughout their 
numerous analyses of feedforward control systems. These 
studies have dealt with the determination of feedforward 
controllers required for two-point composition control by 
manipulation of reflux and steam flows (69,79), reflux 
flow and feed plate location (72,123), as well as the de- 
termination of feedforward controllers for columns which 
exhibit inverse responses (71). Luyben has fit the feed- 
forward controllers used in these studies with simple 
combinations of first and second order lags. Simulation 
results obtained by Stull (127), using the linearized 
model of Gerster et al, indicate that very simple dynamics, 


such as first order lag for feed flow disturbances and no 
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dynamic compensation for feed composition disturbances, 

are adequate for controlling most columns. Stiso (125), 
using the comprehensive model of Huckaba et al, has also 
Shown the advantage of using simple dynamic compensation 


over static feedforward control.. 


Janis (56) has compared the feedforward controllers 
predicted using the linearized model of Gerster et al and 
calculated using Rippen and Lamb's stepping procedure with 
those calculated from the open loop responses determined 
experimentally using the pulse testing procedure. The 
two methods gave good agreement for low frequencies; with 
major deviations occurring at frequencies above approxi- 
mately 1.6 radians/dimensionless time. (The dimensionless 


time was defined as total holdup/reflux flow). 


Cadman et al (19,20,21) have also used the stepping 
procedure of Rippen and Lamb to calculate the open loop 
dynamic response of a multicomponent distillation column 
from the linearized model of Gerster et al. These open 
loop models were then used to calculate feedforward con- 
troLlers 10M two-poing. control. and single point overhead 
and bottoms composition control using the method outlined 
by Bollinger and Lamb. The effectiveness of the various 
feedforward control schemes were studied on a simulated 


five plate distillation column separating a ternary mixture. 
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A number of other attempts have been made to experi- 
mentally illustrate the effectiveness of feedforward 
control as app lted to a distri tatron colummr I,'fo a0, ol , 
79,89). During these studies, the feedforward controller 
was not directly interfaced with the process. The control 
strategy was calculated for a specified disturbance off- 
line, usually on a large digital or analog computer, 
using the various Feedforward models proposed, and then 


implemented manually. 


A noteworthy exception was the study presented by 
Wardle and Wood (132). These workers have studied the 
effectiveness of feedforward control as applied to an 
industrial scale column. Both single-point control of 
overhead composition by manipulation of either reflux 
or boilup flow rate, and two-point control of overhead 
and bottoms composition by manipulation of both reflux 
and boilup flow rates were examined for disturbances 
in the feed composition only. The feedforward control- 
lers were calculated from a linearized model similar 
to that of Gerster et al which had previously been veri- 
fied experimentally (142), using the stepping technique 
of Rippen and Lamb. The Bode diagrams which were cal- 
culated were fit with polynomial transfer functions; 
third order lag for the single-point controller; third 


order lag plus first. order lead for the two-point control- 
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ler, using the complex curve fitting procedure outlined 

by Levy (63). The polynomial transfer function control- 
lers were implemented using pneumatic analog elements. 

A comparison of the static and dynamic feedforward control 
was made by comparing the controlled response to the 


column's open loop behavior. 


Under feedforward control, small offsets occurred 
in the overhead composition indicating that the deter- 
mination of the steady state gain of the feedforward con- 
troller was not sufficiently accurate. The large initial 
error obtained under static feedforward control could 
be sharply reduced using the calculated dynamic feedfor- 


ward controller. 


Very few articles have. been published concerning the 
implementation of combined feedforward-feedback control 
systems. Both Shinskey (115) and Lupfer (66) mention 
that feedback trim has been added to some feedforward 
control systems but give very few details. Luyben (70) 
has analyzed a feedforward controller with feedback 
trim from an intermediate plate temperature controller, 
but gives only scant sqinaawude results. Maselli and 
Miller (83) have outlined a philosophy of feedforward 
control with feedback trim which is used at Sun Qil. 


This method uses a feedback trim from an intermediate 
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tray temperature controller whose set point is reset 
from an overhead composition analyzer. This technique 
has proven successful on a number of columns, including 
deethanizers, butane splitters, depropanizers and de- 
butanizers. However, very little operating data is 


given. 


Wardle and Wood (132) added feedback trim to the 
Single-point feedforward controller controlling the over- 
head vapour composition. They found that in addition to 
eliminating the steady state offset caused by inaccuracies 
in the model, it also acted as a dynamic element in the 


control loop, reducing the need for any further dynamic 


compensation. 


In summary, it can be stated that only sparse ex- 
perimental studies are available in the literature des- 
cribing the effectiveness of feedforward control systems. 
Although some industries appear to have accepted feed- 
forward control as a standard control strategy (14,83, 
121), the details of the majority of the systems installed 
have not been published (113). Therefore, reliance must 
be placed on the many university studies which are mainly 
based on computer simulations for the majority of the 
information necessary for the synthesis and implementation 


of such systems. 
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Z.0 Pulse Testing 


The theses of Wildman (135) and Lees (62), the 
recent paper by Pollack and Johnson (97) and the mono- 
graph of Hougen (46) concerning the pulse testing pro- 
cedure were consulted during this study. No attempt 


to survey the literature in this area was undertaken. 
2.6 Complex Curve Fitting Techniques 


Stroble (126) and Bosley and Lees (15) have recently 
surveyed the various methods available to fitetransfer 
functions in the frequency domain. Chanh (23) has also 
made preliminary comparison between several complex 
curve fitting techniques, including that presented by 
Levy (63), and fitting of transfer functions to time 
data using the search procedure proposed by Rosenbrock 
(110). He concluded that curve fitting in the time do- 
main was the preferred method. When compared amongst 
themselves, no obvious advantage could be claimed for 
any one of the frequency domain complex curve fitting 


techniques studied. 
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CHAPTER 3 
PROBLEM FORMULATION 


wei @rntroduction 


Ana's Sections contains arbrief idescription of: 
i) the pulse testing technique and the 
subsequent curve fitting of the 
frequency response values 
ii) the general expression for ideal 

feedforward control and the modification 
necessary for its application 

iii) the use of various loop record functions 
available in the direct digital control 
(DDC) program to implement the different 


forms of the feedforward controllers used. 


3.24 Pulse Testing 


Since a complete mathematical description of the 
pulse testing procedure has been given elsewhere (62,135), 
only a brief discussion will follow. By definition, the 
transfer function of a process, assuming the initial 
conditions are zero, is given by the Laplace transform 
of the output divided by the Laplace transform of the 


input. 
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The frequency response of the transfer function may be 


evaluated by substituting 
Ss = 1w 


into Eqita\t? on a(dezb)idaoti wisig 
co -jwt 
S y(t)e dt 
G(iw) be 
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Substituting the Euler relation for e- 


el fPE ane cos(wt) ~ 


rex(tye tat 


Tos 


sin(wt) 


302) 
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(3.4) 


and replacing the infinite integration limit with a finite 


value T, Equation (3.3) may be reduced to a complex ex- 


pression of the form 


AA + iBB 

Glia) sbtets-as—= 

CC 0D 

iL , 
where AA = JS, y(t)sin(wt)dt 
BB. = sly (t)cos(wt)dt 
cc = sx (t)sin(wt)dt 


DD = s'x(t)cos (ut) dt 
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An on-line pulse testing analysis program (PTAP) 
(135) is available on the IBM-1800, which calculates 
the frequency response over the frequency range of in- 
terest using Equation (3.5) from the transient response 
time records for the input and output variables. The 
program also contains the option of evaluating the in- 
tegrals using three different methods: 

(a) Filon's method 

(b) Trapezoidal rule 

(c) Fast Fourier transform 
as well as.either listing. the frequency response values 
or plotting them as a Bode diagram. Generally, the tran- 
sient pulse test data was also saved on cards for later 


off-line analysis using the pulse test programs. 


Wildman (135) also outlined a series of guidelines, 
which should be followed during the design of a pulse 
testing experiment. His recommendations were: 

(a) obtain an approximate value for the dominant 

time constant (t) of the process response 
to be studied; 

(b) select a pulse duration that is of the same 
order of magnitude as the approximate time 
constant; 

(c) employ a pulse height of such magnitude that 


a significant-measurable response results, 
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yet sufficiently small to cause a minimum 
excitation of any TOMWINCaaT ties in the 
process response; 

uUse@eacsharp pointed pulse 17 possible, 
however, a rectangular pulse shape 

Should be adequate for most industrial 
testing; 

the frequency range to be studied should 
include the decade above and below the 
decade containing the frequency value 
given by the reciprocal of the approximate 
Eime? cons tant: 

the time record of the response should 
contain at least five approximate time 
constants of data in addition to appro- 
ximately one time constant of initial 
steady state data; 

the sample interval for the data is 


given by 
Osler 
ih = 03:96) 


W 
max 


while the total number of data points 


to be ¢otshected is calculated from 


3 T (3.7) 
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3.3 Frequency Response Curve Fitting 


The mathematical development, which follows, is 
based on the outline presented by Chen and Haas (24) of 
the method proposed by Levy (63). The transfer function 
is defined in the frequency domain as a ratio of poly- 


nomial expressions of the form 


2 
ayta, (iw)ta iw) cea (uy) 


( 
G(iw) = Side aid eciadae a (3.38) 


eee ae 
by tb, (iw) +b iw) +...+b (iw) 


| 
In order for the transfer function to be physically 

realizable, the order of the numerator must be less than 
or equal to the order of the denominator. Equation (3.8) 


can be rearranged to express both the numerator and deno- 


minator as a single complex number. 
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2 4 
O = (b,-b,w +b ,w 78.55) (3.9c) 
2 4 
Y = (b,-b,w +b ew era de) (3.9d) 


If F(iw) is a function, having a real part R(iw) and an 
imaginary part I(iw), which is defined to fit the experi- 
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Fd ii wi). 6 ) Rd be +0 PAF dred) (3.10) 


then an error function can be defined as the difference 
between the perfect fit transfer function F(iw) and the 


‘best fit' transfer function Gliw) as 


E{iw) =: Fliw) = Glin) Ress] uh) 
= Fite) = NOM ia) (acl tan 
DEN(iw) 


Multiplying both sides of Equation (3.11) by DEN(iw) gives 
E(iw) DEN(w) = DEN(iw) F(w) - NUM(w) C3n 125) 
which may now be expressed as a complex number 
E(iw) DEN(w) = r(w) + is(w) Ga tlees) 
The magnitude of the modified error term is 


ISEVi OG)? DENUm) |= 9Yore (a) + 64a) (37.14) 


A weighted error function, E, is defined by squaring the 
modified error magnitude and summing it over all frequencies 


of interest. 
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Equation (3.16) can be differentiated with respect to 
each of the unknown coefficients and the results set to 
zero, which will minimize the error function. A set of 


linear algebraic equations result of the form 


= 
bs 
Cc 
Oo 
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WW 
=< 


OMe love) 


————e ———d 


whose solution may be obtained using standard matrix 
procedures, giving the coefficients for the polynomial 
transfer function. A program called LEVY has been written 
to solve Equation (3.17) for the best fit coefficients given 

(a) the experimental amplitude ratio, phase 

angle and frequency values 
(b) the number of experimental points 
(c) the degree of the polynomial model required, 


n, Mm. 


A relatively unsuccessful attempt was also made to 
include the determination of a best fit time delay. The 
time delay does not contribute to the amplitude ratio, 
only to the phase angle of the frequency response. Given 
the magnitude of the time delay, Thies the experimental 


phase angle was modified by 


> = a - Dok (3540) 


The amplitude ratio and modified phase angle values 


were curve fit using LEVY. The program would attempt a 
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Simple search to find the ‘best fit' value of the time 
delay yielding the minimum residual variance in the 
amplitude ratio data. However, the results obtained 
indicated that the evaluation of the time delay by this 


procedure was unreliable. 


The major drawback in using the complex curve fitting 
procedure to estimate transfer function models is the 
evaluation of the significance of the fit. How good must 
the fit be in the frequency domain to yield a good appro- 
Ximation to the actual time response? Shoneman and Gerster 
(119) have recently suggested that fitting the amplitude 
ratio to within 10 decibels and the phase angle to within 
40°, gave an adequate representation of the transient 
response. Although Wardle and Wood (132) considered the 
transfer functions fit using the Levy procedure to be 
sufficiently accurate to use in their control system ana- 
lysis, they did not outline the criteria upon which this 


judgement was made. 


3.4 Ideal Feedforward Control 


Bollinger and Lamb (12) have presented a generalized 
procedure for calculating ideal feedforward controllers 
given the plant transfer matrix. Assuming the distillation 


column dynamics can be represented by a series of linear 
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time invariant equations between the output and input 
variables, the distillation column plant transfer function 


matrix can be represented as 


X G G G F 
D BE. Pk ae EP e ee PS 
X ebapa eh Gen edeel ecu thee an nee mmo Xe 
B PE. PONE “PT "PL. “PR “PS 
. 
D : c G G f 
DEE =, DX) oe Tee Diee OR) US : (3249) 
B G G G G 
Bel BY.) Selmi BRE BS R 
0 G G G G G G S 


Equation (3.19) has the form 


OUTPUT = PLANT TRANSFER MATRIX INPUT (35720) 


Using the control strategy chosen for this project, 
namely to maintain control of the overhead product compo- 


sition, Equation (3.19) may be simplified to 


Substituting the control strategy, (ie Xp = 0) into 


Equation (3.21), and solving for the reflux flow rate gives 
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This equation describes how the reflux should be ad- 
justed for disturbances in various input variables. Since 
only feed flow disturbances were studied during this pro- 
ject, the ideal feedforward controller of interest is des- 


Cribed by 


g |. Ser OPE (3.23) 
ff GppR 


By a similar analysis, the feedforward controllers 
can also be derived, which would maintain both the over- 
head and bottoms product compositions constant. In this 
case, the model description (Equation (3.19)), can be 


simplified to 
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x Seeng die ROME 
D Ap Hoe esp i ; 
2 | | : i (3524) 
XB Spe px pg. 1. 
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Substitution of the control strategy, (ie X = X = 10) 


and solving for the manipulated variables R and S$ gives 
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The ideal feedforward controller describing how both 


reflux and steam flow should be adjusted for feed flow 
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disturbances is 
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3.5 Implementation of the Ideal Feedforward Controller 


The block diagram of the feedforward control scheme 
implemented during this project is given in Figure 3.1. 
The response of the reflux control loop to a change in 


its setpomims given by 


G G 
R 
1 ees ee (3.227) 
1+G G 
sp RC RV RM 
In order ‘to’ satisfy ‘the @ontrol criteria, Cie X. =%0), 


D 
it is obvious that the effect of a disturbance in feed flow 


on the overhead composition, given by 


X = 
D Gor F (orice) 


must be eliminated by an equal and opposite reaction from 


the controller 
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G R 


VM hd tal dais ina Ae iniy (3.29) 
FR sp 


This may be expressed as 


r G G G G F 
Gor F = - FM ee R PR (3.30) 


Solving Equation (3.30) for the feedforward controller, 


Gee, gives 


G ‘ PE RM (3-531) 


This expression can be further simplified using the 
following assumptions based on the observations of the 
response of the reflux control loop and the flow measure- 
ment transducers: 
(a) the gain of the closed loop transfer 
function of the reflux flow control loop 
LO) ae Set. point. changes 1s 1.05 
(b) the dynamics of the closed loop response 
of the’ reflux flow control o0p- 1s fast 
compared to the response of the column, 
(c) the dynamics of the flow i Wikentrns 
transducers are also fast with respect 
to the response of the column. 
Equation (3.31) can then be reduced to 
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In order to implement the ideal feedforward controller, the 
steady state gain must be modified to reflect the different 
calibration ranges exhibited by the reflux and feed flow 


rate transducers. 


3.6 Loop Record Functions Available in Direct Digital 


Control (DDC) Program 


This section will describe some of the functions 
available in the direct digital control (DDC) program 
used to implement the feedforward and feedback control- 
lers. Since a complete description of the control pro- 
gram exists elsewhere (29), only those functions of direct 


interest will be outlined. 


3.6.1 Exponential Filter 
OUT = (1-FF) INPT + FF MEAS (3133) 
where QUT - filtered value of measurement at this 


sample instant 


INPT - raw unfiltered measurement at this 
sample instant 


MEAS - filtered value of measurement at the 
previous sample instant 


FF - filter constant (O<FF<1) 


3.6.2. Proportional iplus wlntegral..Controller 


ERROR = (MEAS - SETPT) + BIAS (S730) 
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KI - integral controller constant 


RESTI:REST2 - sum of past integral contri- 
butions to the control output (reset) 
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desired ratio 

miscalibration factor, ratio .of 
gains of wild measurement to 
controlled measurement 

current filtered wild measurement 


ratio of zero offset: to gain of 
wild measurement (generally zero) 


ratviovof zero offset. to gain, of 
controlled measurement (generally zero) 


set point of controlled measurement 


3.6.5 Engineering Units Conversion 


OUTP secg 
where QUT - 


MEAS. - 
A x 
B z 


MEAS +B U353.9)) 


value of the current measurement in 
engineering units 


current value of filtered measurement 
twice span of calibration 


calibration offset 


3.6.6 Data Accumulation Loop Record 


The DDC program will allow the accumulation of histo- 


rical data in a special type of loop record, with the 


Oldest Valueo-in the. last position. This loop can, then be 


used to give a 


time delay. 
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The DDC program can obtain the input for a loop 
record from any of the specified locations: 

a) multiplexer 

b) measurement of another loop record 

c) output of another loop record 

d) oldest historical measurement from a 


data accumulation loop record. 


3.7 Implementation of Control Schemes Using 


Loop Record Functions 


This section will present a brief description of 
how the various loop record functions are combined to 
implement the various control schemes employed. Table 
3.1 contains the various loop record functions used by 
each loop record, while a complete listing of all loop 


records defined is available in Appendix A. 
3.7.1 Feedback Control 
The expression describing an ideal proportional 


plus integral controller is given by 
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ts. 


CABUES3< 1 


loop 
record 
number 


0170 
ND 
0250 


0251 


0252 
0260 


0261 
0262 


ALGORITHMS FOR LOOP RECORDS USED 
TO IMPLEMENT THE CONTROL SCHEMES 


loop record functions used 


data acquisition of feed flow 

data acquisition of overhead composition 

input from measurement of 0170 

time delay algorithm 

input from either measurement of 0170 or 
oldest value in 0250 

proportional plus integral control 
algorithm 

exponential filter algorithm 

input from output of 025] 

ratio algorithm 

INPUT HOMeowLpUtLOT Ucod 

proportional plus integral control 
algorithm 

input from measurement of 0175 

proportional plus integral control 

input from output of 0260 

proportional plus integral control 
algorithm 

Output to current output station 
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This function is implemented on the IBM-1800 using the 


proportional plus integral controller function represented 


by 


Gi(t) = .KC(1+KIAtze) (3.41) 


Figure 3.2 illustrates the values of the important loop 
record parameters for the feedback controller. The compo- 
sition measurement is obtained from the multiplexer. The 
calculated output was applied to the current output 
station (COS), whose current output acts as the process 
set point. A feedback controller requires only a single 


loop record. 


3.7.2 Gain Feedforward Control 


The transfer function describing a static feedfor- 


ward controller is given by 
G = K (3.42) 


Four loop records were required to implement this control- 
ler using the loop record functions illustrated in Figure 
seo.) ine firsts loop record calculates» the deviation of the 
feed flow rate disturbance from its reference steady state 
value using the proportional control algorithm with a unit 
gain. The second loop, set up using the ratio algorithm, 


calculates the necessary reflux correction required. The 
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FIGURE 3.2 LOOP RECORD ORGANIZATION FOR FEEDBACK CONTROL 


FIGURE 3.3 LOOP RECORD ORGANIZATION FOR GAIN 
FEEDFORWARD CONTROL 
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thira loop.) setup as,a proportional controller ,.addsethe 
reriux correction. to thesreflux set point. senis loop 
also adjusts the value of reflux prior to sealculating tne 
square of the reflux in the fourth loop. The output of 
the fourth loop, a signal which is proportional to the 
reflux squared, is cascaded as the set point of the 


reflux flow controller. 


3.7.3 Gain plus Time Delay Feedforward Control 


The transfer function describing this controller is 


given by 


= cépe , 
See Kee® ah Si 


Five loop records are required to simulate this controller 
using loop record functions. The first loop is a special 
data accumulation loop record. The length of the loop, 

N, 18S specified as the closest integer value which satisfies 


the following relation 
Ne oydkts «¢ Fon T (3.44) 


The remaining loop records are identical to those described 
previously for the static feedforward controller, except 
that the second loop, which calculates the deviation of 
the feed flow disturbance, obtains its measurement from 


the oldest historical value in. the data accumulation loop 
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record, rather than directly from the multiplexer. This 


arrangement OW loop recowds ts Tllustratedyinsrigure 3.4. 
34/.4) Gain plus Time LagsFeedforward Control 


The transfer function describings thast cohtyoller is 


expressed a's 


K 


tstl 


Four loop records are required to implement this controller 
using standard loop record func tions.....i[hese. 1oop.cecords 
are identical to those described previously for gain feed- 
forward controll except ine this caseathee-first loop app! tres 
an exponential filter to the measurement it obtains from 
the multiplexer. The Pt an between the filter factor 
(FF) and the time constant (t) is given by Jacobson 
and Fisher (54) to be 

At 
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Thisvexpression is valid forsfi iterffactors in the range 
Ded Soh Fe) 0 ( Sy) 


Figure 3.5 illustrates this arrangement of loop records. 
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3.7.5 Feedforward-Feedback Control 


Combined feedforward and feedback control is imple- 
mented with a gain and a gain plus time lag feedforward 
controllers. The loop records required are illustrated 
in Figures 3.6 and 3.7 respectively. These loop record 
configurations result from the combination of those re- 
quired for feedback roWeyo | and those required for the 
feedforward control. The only change required is that 
instead of applying the output of the feedback loop di- 
rectly to the current output station, it is used to bias 


the output of the feedforward controller. 
3.7.6 Qther Feasible Controllers 


A first order lead/lag feedforward controller could 
be employed by combining loop record functions for ex- 
ponential filtering and a proportional plus derivative 
controller. The time constant of the lag element is still 
given by Equation (3.46), while the time constant of the 


lead is given by 
Pye Ka it (3.48) 


Some preliminary tests were performed using this control- 
ler, but when the results using the much simpler time lag 


model proved satisfactory, further testing of this control- 
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LOOP RECORD ORGANIZATION FOR GAIN FEEDFORWARD 


CONTROL WITH FEEDBACK TRIM 


FIGURE 3.6 


54 


he 


Le eit fone 

= ° 3 
ey 2 | 
jglslzlstslais) 


4 
pov 
i 
i] 


oo fu fa iow 
me fat | Oo | oe 


GSAWAOTOTI32 UIAD HOF WOLTANIMAOAO GHOIZA GOs a.€ 
MIST NOAGCIIA HTIW JONTHOD 


0261 


N 
Ke) 
NI 
Oo 


025] 


SP 


0262 


. 
wn 
oO 


Wi QF WE FF WT O 
=f Wi DAF oF; wut] 


Mi rey]uy oO 
milrlojy wi =~ 


LOOP RECORD ORGANIZATION FOR GAIN PLUS FIRST ORDER 


FIGURE 327 


55 


TIME LAG FEEDFORWARD CONTROL WITH FEEDBACK TRIM 
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ler was discontinued. 


It appears that a general transfer function of the 


form 


S ee ee ee (3.49) 


could be obtained by using a series of n_ proportional 
plus derivative controllers in series with m exponential 
filters. This controller would be restricted to those 


transfer functions in which no complex roots appear. 
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It is important that the loop records be calculated 
in the proper order, and not be allowed to get out of 
phase. This can be accomplished by adjusting the phase 
time parameter of the loop record. The proportional 
controller algorithm could also have been used in place of 
the ratio algorithm, which, in some cases, would have 
saved a loop record. This method was not ysed because use 
of the ratio algorithm allowed the storage of the feed- 
forward gain and the miscalibration constant parameters 


separately. If the proportional control algorithm were 
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used, it would necessitate lumping these terms into a 


single parameter. 


Although the digital computer is actually a sampled 
data controller, the analysis which has been performed 
has assumed that the computer acts as a continuous control- 
ler. This assumption is justified since the sample time 
of 16 seconds is small compared to the dominant time 
constant of the column, which is of the order of 600 


seconds. 
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CHAPTER 4 
COLUMN DESCRIPTION AND OPERATION 


4.1 Introduction 


This chapter contains a brief description of the 
distillation column and its ancillary equipment, a 
tabulation of the process variables that were monitored, 
a discussion tot the data transducing required between 
the process and the computer and a discussion of the 
control schemes employed during this investigation. 

The operating procedures required for column start-up 
and shut-down and the calibration techniques used are 
also outlined. The chapter is concluded with a discus- 
sion of the initial steady state conditions selected 
and the material balance errors of closure typical of 


those obtained throughout the study. 
4.2 Distillation Column and Ancillary Equipment 


A nine inch diameter, eight tray glass distillation 
column separating a binary mixture of methanol and water 
was used to study the effectiveness of the various control 
schemes studied during this project. Each plate contained 
four 2 1/4 inch by 1 7/8 inch bubble caps mounted in a 


square pattern, 1 1/2 inch diameter inlet and outlet weirs, 
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and a one inch diameter downcomer. Each plate was also 
equipped with a feed location inlet, two thermocouple 
slots for measuring the plate temperature, a toggle- 
valve sample line for obtaining a sample of the plate 
liquid and two sample points used for the continuous 
plate composition sample (not installed for the present 


project (128)). The plate spacing was 12 inches. 


Four three foot diameter tanks are available to 
contain the fresh feed and each of the two product 
streams. These tanks are piped to allow either con- 
tinuous or batch mode operation. Normally the methanol- 
rich product is stored in one of the tanks, labelled 
the top product tank, and the water-rich product is 
stored in another tank, labelled the bottom product 
tank. This leaves two tanks available for storage of 
fresh feed. Two fresh feed tanks are normally required 
for studies concerning feed composition disturbances. 
For batch operation, the desired feed composition is 
prepared in one of the feed tanks by mixing the required 
amounts of water and methanol from the product tanks. 
The contents are circulated through the feed pump and 
back to the feed tank to ensure a uniform composition. 
While the column is in a batch operation, the feed to 
the column is pumped from the feed tank, while the pro- 


ducts are withdrawn from the column to their respective 
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product tanks. For continuous operation, the feed is 
withdrawn from the feed tank and the products are re- 
turned to that feed tank. At present the feed can be 
introduced on either the fourth or fifth tray from the 


bottom through a small shell and tube feed preheater. 


The reboiler is a vertical basket type with con- 
densing steam on the shell side. It consists of 38 
1/2 inch tubes, two feet in length with a three inch 
diameter central downcomer to facilitate liquid recir- 
culation. The tube sheet is six inches in diameter. 
The section above the tube sheet, the glass vapour 
return line and an extra glass section below the bottom 
plate, act as the vapour/liquid disengaging volume for 
the reboiler. The bottom plate of the reboiler contains 
facilities for installation of a thermocouple, the low 
level leg of the level differential pressure cell, the 
bottom product withdrawal, a liquid sample toggle valve, 
a capacitance probe and its temperature compensation 
thermistor, (which were not implemented during this pro- 
ject (128)) and the outlet and inlet for a small con- 
tinuous sample to be withdrawn through an on-line process 
gas chromatograph. The bottom product, withdrawn from 
the reboiler base, is returned to the appropriate tank 


through a series of small coolers. 
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The vertical nine inch diameter glass condenser 
contains as many 1/2 inch diameter U-shaped tubes as the 
volume would allow. The cooling fluid is water on the 
tube side. The reflux accumulator consists of a nine 
inch diameter to 1 1/2 inch diameter glass reducer at the 
base of the condenser. The base plate is equipped with 
taps for a capacitance probe and its temperature compen- 
sating thermistor, a thermocouple, the high and low 
legs of the level transducer and a liquid sample toggle 
valve. The overhead product is split with the reflux 
being pumped back to the top plate through a small shell 
and tube reflux preheater and the overhead product being 
pumped back to the appropriate tank through a series of 


small coolers. 


The overhead composition was measured with a small 
parallel plate capacitance probe connected to a Foxboro 
Dynalog capacitance recorder/controller. The plates of 
the probe are approximately 1/4 inch square with a gap 
of about 1/8 inch. The plates are imbedded in a plug of 
teflon approximately 1/2 inch diameter. The capacitance 
analyzer measures’ the difference in the dielectric con- 
stant of the mixture due to changes in the composition. 
Since.the.dielectric; constant.jis;also. sensitive, to.tem- 
perature variations, the recorder is equipped with a 


thermistor which detects the temperature so that compen- 
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sation can cancel this effect. The basic theory of this 
method and detail concerning the design of the capaci- 


tance cell are available from Svrcek (128). 


The bottom composition is analyzed every 150 
seconds using an on-line industrial gas chromatograph. 
A small stream is circulated from the reboiler bottoms 
past the chromatograph sample valve and returned to 
the reboiler. The chromatograph is under complete control 
of the IBM-1800 computer. The computer generates a 
control signal causing the chromatograph to take a sample. 
The computer then monitors the detector current, detects 
the peaks, calculates the compositions, presents the 
report, and then commences another analysis cycle by 
commanding that the next sample be taken at the appro- 
priate time. The chromatograph may also be used to 
obtain an analysis of the feed composition, when necessary. 
Operating information and procedures as well as a de- 
tailed description of the chromatograph system are avail- 


able from Berry (9). 


All flows are measured with flange tap orifice 
plates. The reflux, feed overhead and bottoms product 
flow rates were sufficiently small that quadrant edge 
(130) orifices could be used. The remaining flows, 


steam and cooling water were measured using standard 


$2 . n 1 
2tdt to ysosds pend StT 390088 
_fosqes 90% Ro npFeeb" Say entavgono 


(88ST) doo4ve mort stie 6 
. eo 


. , j - so mrs bs 
Oef yxoeve bosyfens ef nottfeor mos MoLIe 
- oo oa 7 
-t Tie ae ake! 
.tqsypotsmords 260 [sfeteubat antl =n ib onte 
oT] 
emogtod vesltodsy sds mov batstwovto ef maevte 
a er 
oj benvutesy bas osvisv siqmse2 iqsipotsmoids & 
‘a : 
’ . es 
fortnos sistqmed Spiny 2Ff fqsrpotsmoyns siT . sl teds 
6 2estsvenep Yeatuqmos ofl .satvamos 0031 -Me 
5 9361 of daBsipotamorda shtd pafeuso fenghe forge 
aN 
rneyyuo votostsbh sAt evotTAom Hens vodut 7 
: oa - ~ 
ait 2tnes2o7q .2hottieoqmos sat 2etstuotss whe 
vd sfoy reviens weitons 2s9nsmmoo non? brid i 


_ 


-ovags edt +6 nedet od olomee Sxon ont Feds ONT ORem 
7 

9a 02 Yam dasypotsmoims snT -omry ats. 

sis to eteylsns ns! ated 
s flaw 26 zouwheso%q bas notssmvoTHt ent. By 

‘ 

frevi %¥65 masteye AqstposeMo do off to noriqro2sb a 

.(@) yxted mort § 

4 om 


Fthvo ast sonelt Adtw bervesem ove 2wolt rt 
org emosttod bis bestitevo best exultes 9aT -2936T 

bs tasrbsup stadt Iloame yl iasrorttye S19W 20167 
-ewolt oninftemeay ett .be2eu od bivos egattrao | 


sbasie prtev bayvesom Stew todew pntlooo bis mB 


63 


Square edge orifices. The differential pressure created 
at the orifice was transduced to a 3-15 psig pneumatic 
Signal using a Foxboro differential pressure cell. The 
levels in the condenser and reboiler were measured with 
low range (30 inch) Foxboro differential pressure cells. 
The column pressure was measured using a Foxboro Pres- 
sure Transmitter with a range of -10 to +10 inches of 
water. All temperatures were measured using iron-con- 
stantan thermocouples. A more detailed description of 
the column and its associated instrumentation is avail- 
able from Pacey (93). A detailed schematic diagram of 


the column is presented in Appendix B. 


It is worthwhile to note that the column, as used 
by Svrcek (128), was disassembled and reinstalled in a 
new laboratory prior to the testing reported in this 
work. During reinstallation, the column piping was 
modified so that various possible control schemes could 
be studied with a minimum of physical effort expended 
in reconfigurating the control system. Most of the pneu- 
matic instrument lines were installed with a flexible 
plastic tubing, which can very easily and quickly be 
changed. Also, all recorders, indicators and control- 
lers were placed in special modules which were installed 
in a series of cabinets. This modular type of instru- 


mentation permitted the standardization of modules which 
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facilitate the borrowing of a similar instrument from 


other equipment in the event of an instrument failure. 
4.3 Process Variables Monitored 


Throughout this project, the data acquisition capa- 
bility of the IBM-1800 DACS computer facility was used 
to measure and record the process variables listed in 
Table 4.1. All variables are available to the computer 
through the multiplexer, but not all variables were 
supplied with a loop record. A loop record makes the 
measurement of a variable available to the direct digital 
control (DDC) program, which performs the bulk of the 
data acqUisition*and_control functions. Other special- 
purpose data acquisition and control functions were per- 
formed by user written programs by either direct access 
to the multiplexers, or by obtaining the measurements 
from the DDC loop records. In adapt the values of 
most of the process variables are also available from 
standard analog indicating and recording instruments. 
Table 4.1 also indicates that sufficient analog control- 
lers are also available to completely control the column, 
and that only a few of the digital controllers have been 
implemented. The column at present can be operated under 
analog control, but a few extra digital control loops 


must be implemented before it.can be run completely via 
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computer control. 


4.4 Data Conversion Between the Process and the 


Computer 


Before the computer can perform any function, it 
must have the process information available in a form 
which it can use. The process variables are available 
infphystcal units, ie ice psig. in H,0 and Ib/min, but 
the computer requires a discrete digital value, ie 7/209. 
Similarly, to make use of the information generated by 
the computer, this digital data must be converted into 
a physical variable in order to drive the final control 
element. This section will illustrate the methods used 
during this project to convert the physical process 
variables into digital values, and the digital computer 
outputs ito. physical control variables. -Figure 4. 1 
illustrates the various steps involved with the data 


transducing process. 
4.4.1 Analog to Digital Conversion 


The first step in converting the process variables 
is to generate a signal, either pneumatic or etectronic, 
which is proportional to the measurement. The flows are 


converted to a pneumatic 3-15 psig signal using 
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an orifice plate and a differential pressure well. fhe 
levels and column differential pressure are also converted 
to pneumatic 3-15 psig signals using low range differ- 
ential pressure cells. An absolute pressure transmitter 
is used to convert the column pressure to a pneumatic 
3-15 psig signal. The overhead composition is converted 
to a pneumatic 3-15 psig signal using a special pen 
position transducer for use with the Dynalog capacitance 
recorder. A 10-50 milliampere current proportional to 
the pneumatic signal is; produced. from a pressure to 
current (P/I) convertor for each pneumatic measurement. 
The current is further converted to a proportional 1-5 
Volt electrical signal across a 100 ohm resistor. All 
temperatures are converted to a 0-10 millivolt electrical 
Signals using iron-constantan thermocouples. These 
electrical signals are transferred to the computer, 
usually on shielded twisted pair wires, and each vari- 
able is assigned to a unique set of terminals on the 
multiplexer. The multiplexer is essentially a bank of 
relay switches. When the computer requires a measure- 
ment, “the appropriate multiplexer relay 1S ‘Closed, ap- 
plying the voltage signal to the analog to digital (ADC) 
converter. The range of the analog to digital converter 
ts normally 0-5 volts: Thus Tow evel electrical Si1gwals:, 
ie 0-10 millivolt temperature measurements, must be 


applied to the analog to digital converter through an 
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| 
appropriate amplifier to produce a proportional 0-5 
volt signal. The continuous voltage signal is digitized 
by the analog to digital converter in the range 0-32767, 
which is the largest integer contained in the 16 bit 
word length of the IBM-1800. This measurement is now 


available to be read into the appropriate core storage 


location to be used by the computer software. 


4.4.2 Digital to Analog Conversion 


When the computer has completed the calculation 
algorithm and is ready to apply the result to the pro- 
cess, it is available in an appropriate core storage 
location as a binary integer in the range 0-32767. This 
digital value is applied to a digital to analog conver- 
ter (DAC), resulting in a proportional 0-10 volt electri- 
cal signal. The computer activates the appropriate pulse 
output (PO) relay which applies the 0-10 volt electrical 
signal to the current ouput station (COS). The current 
output station is essentially a zero order hold device 
which maintains the input signal until changed the next 
time. The current output station converts the applied 
voltage to a proportional 10-50 milliampere electrical 
signal. The 10-50 milliampere electrical signal is further 
converted to a proportional 3-15 psig pneumatic pressure 


signal, which is used to drive the final control elements. 
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The final control elements in this project consisted of 
the set points of the pneumatic flow control systems. 
These signals could just as easily have been applied 


directly to the control valves. 


4.4.3 Comments 


The major problem associated with signal conversion 
and transmission is noise pick up, which should be avoid- 
ed if possible, or at least minimized. The transmission 
distance of both current and pneumatic signals should be 
minimized, using high level voltage signals, where pos- 
sible, to°transport process information over longer 
distances. The high level voltages are preferred for 
data transmission because of the smaller amount of sig- 
nal deterioration experienced due to noise pick up than 
with a low level signal. When low level signals must be 
used, ie thermocouple voltages, care must be taken not 
to run the transmission lines past noise generators such 


as fluorescent lamps, etc. 


Pneumatic instrumentation is not normally recommen- 
ded for computer control applications, since equivalent 
electronic instrumentation is available which can elimin- 
ate some of the required data conversion, ie pneumatic to 


current. Less opportunities for signal distortion can 
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result from these fewer conversions. However, the pneu- 
matic column instrumentation already existed and the 


data acquisition system was designed using them. 


4.5 Description of the Control Schemes Implemented 


4.5.1 Degrees of Freedom 


According to Smith (122), the number of degrees 
of freedom available in defining the steady state opera- 
tvon Of ea distillation column witha total condenser. and 
a partial reboiler are 

DF = COMP + 2 NSTAGE + 9 (4.1) 
However, Howard (47) has suggested that, for control 
purposes, this expression is inadequate since it neglects 
the additional transient degrees of freedom related to 
the holdups on each stage. This would add an additional 
NSTAGE + 1 degrees of freedom, so that the total degrees 
of freedom would become 


DF = COMP + 3 NSTAGE + 10 (4.2) 


Since the control engineer is concerned with an 
existing column, the following variables are generally 
fixed by the design of the column: 

pressure drop in each stage NSTAGE 


pressure drop in reflux divider ] 
heat loss in each stage (except reboiler) NSTAGE - 1 
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heat loss in reflux divider ] 
holdup in each stage (except reboiler) NSTAGE - 1 

feed plate location 1 
total number of plates ] 

iO tehs: SNS TAGE 42 

In addition to these variables, the reflux temperature is 
generally fixed by the operation of the condenser while 
the feed conditions are usually determined by the opera- 


tion of the downstream processes, fixing the following 


degrees of freedom: 


feed conditions COMP mre 
reflux temperature ] 
Oa i: COMP + 3 


Of the COMP + 3 NSTAGE + 10 degrees of freedom 
available to define the column operation, COMP + 3 
NSTAGE + 5 degrees of freedom have now been fixed by the 
column design and operation, leaving only five variables 


available for specification by the control engineer. 


Of these specifications, one is used to define the 
column pressure, while two are used to define the levels 
in the condenser and the reboiler. The remaining two 
degrees of freedom are used to maintain the material and 
energy balances. The specification of these degrees of 
freedom, as defined in this study, include: 

a) cooling water flow rate (column pressure) 


b) overhead product flow rate (condenser level) 
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c) bottoms product flow rate (reboiler level ) 


d) reflux flow rate (material balance) 
e) steam flow rate (energy balance) 


With these specifications, the column operation is 
completely defined. The control scheme implemented 


is illustrated schematically in Figure 4.2. 
4.5.2 Material and Energy Balance Control 


Control of the distillation process requires the 
maintenance of the material and energy balances. The 
various control schemes proposed to fulfill these re- 
quirements have been reviewed by many authors (16,40,44, 
84,89,115). From the many control schemes suggested, 
only two basic control philosophies emerge, dubbed 
direct and indirect material balance control. Figures 
4.3 and 4.4 illustrate simple examples of indirect and 
direct material balance control schemes. In the indirect 
material balance control scheme, the energy balance is 
manipulated by the steam flow rate and the material 
balance is adjusted indirectly by the reflux flow rate 
through the reflux accumulator level controller. Both 
product flows are then adjusted under level control to 
maintain the column inventory. In the direct material 
balance control scheme illustrated, the energy balance is 


controlled directly by the steam flow rate, while the 
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FIGURE 4.3 INDIRECT MATERIAL BALANCE CONTROL 


FIGURE 4.4 DIRECT MATERIAL BALANCE CONTROL 
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material balance is controlled directly by manipulation 
of the overhead product flow rate. The reflux and bot- 
toms product flow rates are then adjusted under level 
control to maintain the column inventory. This control 
scheme could easily be reversed with the reflux flow rate 
defining the energy balance and the bottoms flow rate 
directly defining the material balance. The steam and 
overhead product flow rates would then be adjusted under 


level control to maintain the column inventory. 


Examination of the many articles suggests a contro- 
versy exists in the use of the direct and indirect control 
schemes. The advantages claimed for the direct material 
balance system (84,89,115) include 

- greater degree of self regulation 

- greater degree of accuracy 

- reduced interactions 

- tighter control 

= sphiprictty? 

Buckley (16,17) has suggested that the main disadvantage 
of the direct method is the large reflux cycles which may 
be set up in the column due to uncontrolled reflux tempera- 


ture fluctuations, resulting in an oscillatory product 


quality. 


Despite the many advantages claimed for the direct 
material balance approach, the indirect method of material 


balance control has remained an extremely popular industri- 
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al control scheme. For this reason, the indirect material 
balance control scheme has been implemented during this 
project. The fact that such a controversy does exist 
suggests that an experimental comparison of both control 


schemes should be made. 


In order to complete the definition of the column, 
the specification of the variables not fixed by the 
column design is accomplished as follows: 

a) the feed and reflux temperatures are 

maintained a few degrees below their 
boiling points by a temperature controller 
which adjusts the steam flow to a small 
heat exchanger; 

b) the feed flow rate, the reflux flow rate 
and the steam flow rate are maintained 
under flow control; 

c) the feed composition is prepared in the 
feed tank; 

d) the overhead and bottoms product flow 
rates are controlled from the levels 
in the reflux accumulator and the reboiler 
respectively; 

e) the column pressure is maintained by adjust- 

ment of the cooling water flow rate to the 


overhead condenser. 
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The overhead product flow rate, the bottoms product flow 
rate and the cooling water flow rate are adjusted by con- 
trol valves, which receive their pneumatic control sig- 
nal directly from their respective controller. This 
method of control is generally considered bad practice 
(115) due to the hysteresis exhibited by the valve. 

The installation of either a valve positioner or a flow 
controller whose set point is cascaded from the primary 
controller will eliminate any effects caused by these 


nonlinearities. 


4.5.3 Composition Control 


Once the operation of the column has been complete- 
ly defined by the material and energy balances, a compo- 
sition control scheme may be specified, which will achieve 
a particular objective. The control objective, chosen 
during this study, was to maintain the overhead product 
composition by manipulation of the reflux flow rate. 

This control strategy is currently implemented using 
either standard analog controllers or the control capa- 


bility of the IBM-1800. 


Various analog controllers have also been provided, 
including an intermediate tray temperature controller 


(23), a differential pressure controller and a ratio 
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relay, although they have not been used during the present 


ScUuGy. 


A block diagram of the analog feedback control 
scheme implemented is shown in Figure 4.5. The feedback 
controller is a special capacitance measuring, pneumatic 
Foxboro Dynalog controller. The output of the composition 
controller is cascaded to the pneumatic setpoint of the 
Foxboro Consotrol controller controlling the reflux flow 
rate. The same control scheme is also illustrated in 
Figure 4.6, except that the function of the analog compo- 
sition controller has been replaced by the IBM-1800 
digital computer. The only difference between these two 
control schemes exists in the control algorithm by which 
each controller calculates its output. The digital com- 
puter uses a digital approximation to the ideal feedback 
controller expression (29), while the analog controller 
approximates this expression with a series of mechanical 
linkages and pneumatic resistances and capacitances. An 
example of a complex expression which has been developed 
to describe the dynamics of the pneumatic Dynalog control- 
ler is available from Svrcek (128). Consequently, a 
correspondence between the control parameters of the two 


feedback control schemes would not be expected to exist. 


In both cases, the feedback control action is 
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FIGURE 4.5 BLOCK DIAGRAM FOR ANALOG FEEDBACK CONTROL 
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FIGURE 4.6 BLOCK DIAGRAM FOR DIGITAL FEEDBACK CONTROL 
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activated by disturbances measured in the overhead product 
quality. The control actions, calculated based on these 
excursions from the set point, are cascaded to the pneu- 


matic set: point..of the reflux-flow rate-control ler. 


The feedforward control is implemented using the 
digital computer as shown in Figure 4.7. The feedforward 
control action is activated by the measured disturbances 
in the feed flow rate as they enter the column. The con- 
trol actions, calculated based on these excursions from 
their set point, are also cascaded to the set point of 


the reflux flow controller. 


Examination of the block diagrams for these two 
control schemes suggests the Rexe Togiical step; the com- 
bination of the two control schemes, which is illustrated 
in Figure 4.8. This control scheme has also been imple- 


mented using the computer. 


The above figures are examples of what is called 
analog supervisory control, in contrast to direct digital 
control (33,53,88). Analog supervisory control allows 
the local analog controllers to control each individual 
loop, while the computer only supplies the set points 
for these local controllers. Thus, if the computer fails, 


the plant can continue to operate by supplying a manual 
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computer 
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FIGURE 4.7 BLOCK DIAGRAM OF DIGITAL FEEDFORWARD CONTROL 


computer 


FIGURE 4.8 BLOCK DIAGRAM FOR COMBINED FEEDFORWARD-FEEDBACK CONTROL 
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set point to each controller. The computer can now spend 
most of its time optimizing the process operation, rather 


than controlling the many individual control loops. 
4.6 Operating Procedures and Calibration 


Since no formal operating procedures existed for the 
column, the start-up and shutdown procedures (93) were 
evolved. It will be noted that the column was not started 
up on total reflux, but rather as soon as the reflux 
accumulator contained a sufficient level, the value of the 
desired reflux flow was set on the reflux controller 
and the column was allowed to come to steady state. 
Operating in this manner insured that the column could 


be tested within two to three hours after initial start- 


up. 


All the column measuring instruments were calibrated 
using the methods outlined in Appendix A.2. The calibra- 
tion results are given in Figures A.1 to A.7 and Table 


A.1 of Appendix A. 
4.7 Column Operation 


The operation of the column was first checked by 


attempting to reproduce some of the steady state data 
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previously reported by Svrcek (128). Typical results of 
these steady state comparisons, which involved taking 
simultaneous samples from all plates, are presented in 


Table A.3 of Appendix A. 


4.8 Selection of the Reference Steady State Conditions 


Since the behavior of the column was to be studied 
for relatively large disturbances (of the order of + 20%) 
in the input variables, feed, reflux and steam flow rates, 
the selection of a suitable initial steady state to allow 
for normal operation after all such disturbances was 
necessary. The initial requirements sought included: 
a) the disturbance should not cause the 
column to flood, weep or exceed the 
condenser or reboiler capacity; 
b) the disturbance should not force the 
column operation through a pinch condition, 
causing extremely nonlinear responses; 
c) the product flow rates should be approxi- 
mately equal since the feed will nominally 
be made up of equal quantities by weight of 


methanol and water. 


The typical steady state conditions used by Svrcek 


(128) were initially employed, but the large increases in 
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the feed, reflux or steam flows caused the column to 
flood. These steady state conditions were adjusted until 
the disturbances did not cause flooding. The reference 
steady state conditions established and the reproducab- 
ility as given by the standard deviation obtained through- 


out the course of this study are presented in Table 4.2. 


The choice of the steady state operating conditions 
represented a compromise based on the requirements pre- 
viously outlined. The responses to the various distur- 
bances exhibit very nonlinear behavior, in that the over- 
nead praduGt composi tic, is dy ivensim.o “Or omt fof ethe 
pinched conditions at the top of the column, while the 
bottoms product composition response is discontinuous 
dt zero percent by weight nfetianol@in-addieitn “to “similar 


pinch proaucing and pinch relieving effects. 


The choice of a steady state, exhibiting nonlinear 
dynamics would, howetes | be expected to give the proposed 
feedforward control schemes a much more severe test than 
would choosing a steady state which exhibits a more li- 


near dynamic response. 
i-J- Material BatancesEyyors of Closume 


Material balance calculations were performed using 
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the data acquisition capabilities of the IBM-1800, in 
conjunction with the user written programs DASS, DATAC 
and BALNC. Table 4.3 illustrates the average materia ' 
balance errors with its standard deviation obtained 
throughout the course of this study, while Table A.4 
of Appendix A gives a specific example of the report 


produced by the material balance programs. 
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CHAPTER 5 
DEVELOPMENT OF AN OPEN LOOP DYNAMIC 
MODEL OF THE COLUMN 


5) pintroduction 


The various methods used in previous studies to 


describe and evaluate the dynamic behavior of the dis- 


tillation system have been outlined in Section 2.2. 


These methods include: 


i) 


ii) 


177) 


deriving the complete set of fundamental 
equations describing the related pheno- 
mena, making as many simplifying assump- 
tions as necessary in order to obtain a 
solution and applying either analytic, 
analog or digital computer techniques 

to obtain a solution 

using an approximate correlation method 

to obtain the parameters of a transfer 
function representation as a function of 

the initial steady state operating conditions 
determining the parameters of an approximate 
transfer function representation using 


standard dynamic testing techniques. 


Since obtaining a solution to the general equations 
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Characterizing the system is not difficult with the 

large digital computers currently available, no longer is 
it necessary to linearize these equations; so consequently 
the nonlinear behavior of the process may be analyzed. 
The lack of sufficient accurate information describing 
such phenomena as plate mixing, plate hydraulics and 
dynamic plate efficiencies, among others, are the major 
limitations to this technique. The ease with which the 
solution can be obtained depends on the number of assump- 
tions made; ie the more assumptions made, the easier the 
solution. The large digital computer required to obtain 
the solution is generally unsuitable for controlling the 
process. Thus, for control purposes, a simplified model 
would still have to be developed, whose parameters could 
be updated periodically based on the off-line solution 


of the large model. 


The advantages of obtaining a reasonably good appro- 
Ximate model based on the initial steady state operating 
conditions are obvious. However, the relationships be- 
tween the column inputs and outputs, as described pre- 
viously, are very complex, suggesting that even if a 
correlation can be found, it too may be equally as complex. 
Although a recent study by Wahl and Harriott (134) has 
introduced a correlation which may have merit, in general, 


the correlation methods developed up to now have not 
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experienced wide spread acceptance. 


Approximate models of a number of columns (50,56,85, 
98) have been successfully determined using various dyna- 
mic testing techniques. The range of the models developed 
experimentally is, by necessity, limited by the restricted 
operating range of the column inputs so as not to cause 
unduly large transients in the column output variables. 
When the initial steady state conditions of the column 
change appreciably, a new determination of the approximate 
model parameters would, ideally, be required. This method 
could, of course, be used with the large fundamental 
model to calculate the parameters of a simple model off- 
line, which could then be used to update the parameters 


of the on-line model for the changing conditions. 


A combination of transient response and pulse testing 
analyses were employed in this study to obtain the steady 
state and dynamic relationships describing the column 
responses. The pulse testing procedure was chosen based 
on the results presented by Lees (62) and Wildman (135). 

A computer program, PTAP (pulse testing analysis program) 
was also available on the IBM-1800, having been developed 
by Wildman (135), which calculated the Fourier transform 
of the input and output variable time data, yielding the 


amplitude ratio and phase angle versus frequency diagram. 
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Due to the nonlinearity of the system, the steady state 
process gains could not be accurately determined from 
the pulse testing data; consequently, a transient response 


analysis was used to determine them. 


In the following sections the data used in the 
determination of the process gains will be presented first, 
followed by the data used in determining the dynamic para- 
meters for a simple first order lag plus time delay ap- 
proximate model. The results of fitting the measured 
frequency response diagrams to more complex polynomial 


transfer functions will also be presented. 


5.2 Determination of the Process Gain 


5.2.1 Procedure 


The feed flow rate, reflux flow rate and steam flow 
rate were the input variables; the overhead product 
composition and flow rate, and the bottoms product compo- 
sition and flow rate were the output variables considered 
during this portion of the study. The column was ini- 
tially allowed to come to steady state with all the input 
variables set to their reference steady state values, as 
listed in Table 4.2. After the column had reached this 


steady state, the values of all the column variables were 
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documented in a material balance report, similar to that 

presented in Figure A.4, prepared by programs DASS/DATAC/ 
BALNC. Details of these programs are given in Appendix G. 
The input variables were then varied, one at a time, over 


the ranges: 


feed flow etc eae OR 0 
reftlTux flow =o | 1957-025 
steam flow =| 2705S 0125 


When the column has reached the new steady state, a se- 
cond material balance report is generated, recording all 
the column measurements. The results obtained using the 
feed composition as an input variable were also measured, 


but were reported elsewhere (94). 


So 2.c ReESUItS 


The experimental operating conditions from which the 
process gains were calculated are given in Tables C.1, 
C.2 and C.3 of Appendix C. The gains, initially calcula- 
ted from the unsmoothed operating data, are presented in 
Tables C.4, C.5 and C.6 of Appendix C. The dependence of 
these open loop gains on the magnitude of the disturbance 
in the various input variables studied is illustrated in 
Figures 6.1, 5.2 ‘and -o.3. A Jeast squares’ fit.or the, pro- 
cess gains, as given by the expressions in Table 5.1, are 


shown in these figures by the dashed line. Examination 
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LEAST SQUARES EQUATIONS FIT TO THE 
EXPERIMENTAL OPEN LOOP PROCESS GAIN DATA 


Reflux Flow 


range. sURGoea Ano 
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of these figures indicates that the scatter about these 
straight lines is appreciable, suggesting that it would 
be more appropriate to smooth the operating data before 
calculating the gains. The operating conditions of the 
column, smoothed by visual examination, are displayed in 
Figures C.1, C.2 and C.3 of Appendix C. The dependence 
of the process gains calculated from the smoothed opera- 
ting data on the magnitude of the input variables is 
ponpared to the dependency calculated previously in 
Pagures. 5..), Sac. and 523. The process. garns, when: calcu- 
lated from the smoothed operating data, exhibit the 
'S'-shaped curve characteristic of nonlinear systems 

when plotted against the input disturbance. It is in- 
teresting to note that, in almost each case, the least 
Squares straight line is a reasonable approximation to 
the calculated nonlinear relation. This is particularly 
true in the case of the process gains for the composition 


Variables. 


Table 5.2 illustrates the reproducability of the 
reference steady state as indicated by the standard de- 
viation obtained throughout the transient tests. These 


errors, although small when compared to the absolute 


value of the measurement, are important when a calculation 


requires the use of the difference between two values of 


the same order of magnitude. Since the process gains are 
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calculated from the ratio of two such differences, small 
errors can cause quite large errors in the calculated 


gains. 


Certainly some of the scatter, particularly in the 
case of the overhead and bottoms product flows, was also 
due to the failure to maintain a constant feed composition. 
Damilatroneot these resulrvs, diven in Tapes C. 1, C.2 
and C.3 of Appendix C, indicated that some of the scatter 
in these process gains appears to correlate with the 
variation in the feed composition used. These gains appear 
to be more sensitive to feed composition changes than do 


the overhead and bottoms product composition gains. 


The overhead and bottom product flow rate process 


gains can be related by the following expressions: 


(feed disturbances) Kye + Kae Sl U5. Va) 

(serial Xy pdrieSsturabia nce S)))> ckea = tikka) =o 0 (5..1b) 
DR BR 

steam disturbances K + K = 0 Saale 

( ) DS BS ) 


The consistency of the overhead and bottom product 
flow rate gains, as calculated from*both the original 
data and a smoothed representation of the same data, have 
been compared using Equations (5.1) in Tables C.1 - C.6 


of Appendix C. Examination of the residuals presented “in 
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these tables indicates that despite the good closure on 
the material balances obtained, generally less than a 3% 
discrepancy, large deviations occur when the gains calcu- 
lated from the original operating data are substituted in- 
to Equation (5.1). Although the variation in the values 
calculated from Equation (5.1) is smaller using the gains 
calculated from the smoothed representation of the opera- 
ting data, the variations are still larger than would be 
anticipated for such close material balance closures. 
These discrepancies reflect some of the difficulties en- 
countered in using the indirect material balance control 


schemes as discussed by Nisenfeld (89) and Shinskey (115). 


5.3 Characterization of the Column Dynamic Behavior 


5.3.1 Procedure 


The pulse testing procedure, designed using the 
guidelines presented by Wildman (135), was used to deter- 
mine the dynamic relationships of the column. The column 
was allowed to come to steady state with the input variables 
set to the values listed in Table 4.2. When the steady 
state conditions were achieved, they were recorded using 
the on-line material and energy balance programs. In 
addition, the system data acquisition and buffering to disk 


programs were initiated to save the data. Collection of the 
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initial steady state data continued for a period of 

about 10-15 minutes in order to define the steady state 
conditions for the program PTAP. A rectangular pulse 

of about 10-15 minutes duration and of sufficient magni- 
tude to cause a measurable response was implemented by 
adjustment of the set point of the input variable flow 
controller. The pulse disturbances were applied to each 
variable in turn, the feed, reflux and steam flow rates. 
The output variables of interest in this portion of the 
study included only the overhead and bottoms product 
compositions. The responses of the plate temperatures 
were also recorded (94) for use during a. subsequent pro- 
ject (23). Each pulse test was repeated a number of times, 
in both the positive and negative directions, in order 

to qualitatively compare the reproducability. The magni- 
tude and duration of the pulse were varied slightly during 
a few tests, but no attempt was made to quantitatively 
analyze the effects, if any. At the completion of the 
pulse test, when the response returns to the initial 
steady state, the steady state conditions were again re- 


corded by the on-line material and energy balance programs. 


The complete set of transient response data was then 
recovered from the system disk files using the program 
GTDAT. A hard copy of the data was obtained on the sys- 


tem printer in addition to a set of punch cards. The 
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Cards, DeSsTdes acting aS a Nistorical y“ecoud, provided 


input to the pulse analysis program. 


The program PTAP was modified slightly from that 
described by Wildman (135) in that the program contained 
an option to punch the calculated frequency response 
on cards, and the frequency response plotting subroutine 
employed by Farwell (32) was substituted for the origin- 


al plotting subroutine. 


Based on the guidelines of Wildman (135), the pulse 
testing procedure was governed by the following character- 
eS CS 

a) the approximate time constant determined 

from the data of Svrcek (128) and by 
preliminary. transient response, sles ts) was 
found to be between 10 and 15 minutes 

b) the corner frequency is between 0.067 

and 0.1 radians/minute 
c) the frequency range of interest is 
0, OD). to 1. 0. radian si/ manne e, 
On the basis of these conditions, it followed that: 
d) the sample interval should be less than 
19 seconds, based on the maximum frequency 
of interest and the Lees-Dougherty cri- 


teria (135) for minimizing the numerical 
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approximation error; consequently the 
closest sample interval available in 
the DDC program of 16 seconds was 
employed 

e) the total time record to be collected 
Should require about 10-15 minutes of 
Steady state data and between 50-70 
minutes of transient data 

f) the Erecuency response values calcu- 
lated below a normalized frequency 
content value of 0.3 are not considered 


reliable. 


Seis. co RESUS 


The reference steady state, about which the pulse 
tests were conducted, along with the pulse characteristics 
are contained in Table D.1 of Appendix D. The complete 
transient response data for all the pulse tests: are 
available elsewhere (94). Examples of typical transient 
responses of the overhead and bottom product compositions 
to both positive and negative pulses in the input vari- 
ables are shown in Figures D.1, D.4, D.7 and D.10. Par- 
ticular note is made of Figure D.7, which illustrates an 
inverse transient response in the overhead product compo- 
sition for decreasing disturbances in the steam flow 
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D.6, D.9 illustrate the frequency response representation 
of increasing overhead and bottom product compositions 
nesponses respectively »awhnke ihigumes eDahl 9.0. ka, CDSL3 
Show the frequency response representation for decreasing 
overhead product compositions. The frequency response 
for a decrease in bottoms product composition was not 
evaluated since the disturbances ae ae the composition 
to pass through zero percent by weight methanol, causing 


a discontinuity in the response. 


Comparison of the frequency response results calcu- 
lated from the repeated pulse test indicates that, by 
visual examination at least, the agreement obtained is 
quite good at the lower frequencies, while major discrep- 
ancies appear to occur in the slope of the resulting 
high frequency asymtotes. These discrepancies are undoubt- 
edly a result of the nonlinear behavior of the column 
response. No special effort had been made to insure 
that an exact replicate pulse would be formed during the 
repeated pulse tests; in fact, for a number of tests the 


pulse duration was deliberately varied. 


The frequency response diagrams, for frequencies 
greater than 0.3 to 0.4 radians/minute, exhibit a great 
deal of scatter characteristic of the pulse testing 


procedure. The reasons for this phenomena have been 


07 


nottsingseorqget s2noqzo% yons 
re 
anorsteoqmoa Jouboxg nodod b a 
er.ag .Sf.0 ,f1.@ esvword me cP svi 
: Pan) 
ontesstoab rot nofisins2stgoy s2nod b ciahe 


Bi a 


senogea?y yoneupet? sat _enota tz0ge 29 99 er 
r$tzoqmod Jouborg zmodiod Pe 


> P * a. 
7]O0N 2SW HO 


ait beeus> 2egonsdwws2Fb ads sonte 


FFe@oqmes on 


we 


on 


ronsdtom tiptow yd sngovs@ Ones cquorts 


.s2noqzot say ar ed tunt so 


feo es fueet g2noqges yonsupayt 94t To ara 
. 7 7 
—_— ~ 7 


iad ,3a6fid artasronr geod 3 


s2fuq batesqey sat mort ft 
ivy: 

ttetdo Tnemeg tes Bans .teseT 36 aotsentmsxe [et 

yztb woOLss ‘dw .2stoneupeyt vswol ont 36 boop 2: 

ss 

cote eft at swoo0 oF 369qg6 29) 

a 

syjzetb o2oHnT .2stodmy25 YoOnsupst 

notveted rssntioen sat Yo tinest 

"Won t oF sbDé6i ty0 TTS [ 663 Sqz ov 3 

ontxub bemvet ed bluow s2iuq stsoriqey t26X9 | 


+0 imun 6 Ot ,t0B8T Nt 22Fees s2fugq. 


bs Yaa 


botvsv ylatervsdiieb 26w nots 54 


“ot ,2mstpsatb sznogzst oneonest ae 
9 

stdtdxea .stuntm\ensrbsy #.0 oF ©. 0 neds * 
ontteat s2fuq eat Te aitatrsioa vents em 


nsod eve stamonsiq etdd tet 2Moenay oat . 1ubs2¢ 


Loy 


discussed in detail by Wildman (135). These unreliable 
frequency response values generally occurred for normal- 
ized frequency content values in the range 0.3 to 0.5, 
but in some cases are as high as 0.7 to 0.9 (94), com- 
pared to the value of 0.3 suggested by Wildman (135). 

The termination of the data collection before the process 
had returned to its steady state and the rounding of the 
corners of the rectangular pulse , imparted by the re- 
sponse of the flow controller to its set point change, 
probably contributed to the excessively. targe values of 


the normalized frequency content cut off point obtained. 


5.4 Approximate Column Model 


5.4.1 Procedure 


The parameters of a first order lag plus time delay 
approximate model were obtained using the following methods: 
a) the transient response was fit using the 
Rosenbrock (110) gradient search technique 
as outlined by Berry (9) 
b) the frequency response was fit using the 
complex curve fitting technique of Levy 
(63) modified such that a best fit time 
delay could be included 


c) the parameters were determined based on the 
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frequency response characteristics 

exhibited by a first order lag and 

time delay functions: 

- the time delay makes no contribution 
to the amplitude ratio 

=ethe tirst order time constant 1s 
defined as the reciprocal of the frequency 
at which the value of the amplitude 
atic: Us «50.407 (20 secipaels ) 

- the time delay is calculated from the 
difference between the measured phase 
angle (¢) and 45 degrees as 


cA (d-45,) desori (5.2) 
Spies 


d) the process gains, based on the magnitude 
of the pulse test, were determined from the 
Gain expressions established bye transrent 


response testing. 


Several of the open loop frequency responses were 
also fit by a series of more complex polynomial transfer 
functions. No attempt was made, however, to determine 
which of these more complex models gave a better approxi- 
mation to the transient response than does the simple 


first order lag plus time delay model. 
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Steer. RESUTtS 


A first order lag plus time delay transfer function 


of the form 


—T.§ 
e D (5.3) 


TStl 


has been assumed to describe the dynamic response of the 
distillation column. This form of approximate model has 
been used successfully in a number of previous control 
studies (9,14,112). The parameters of this model, deter- 
mined for the various pulse tests performed, are given 

in Tables 5.3, 5.4 and 5.5. Average transfer functions 
are presented in Table 5.6, describing the response of 

the composition output variables to each of the flow 

input variables studied based on the values determined 
from the time domain curve fit. The transient response 
calculated from these simple models are compared with 

some actual transient responses measured in Figures D.1, 
D.4, D.7, D.10 of Appendix D. In addition, the frequency 
response of these approximate models is also compared 

with the frequency response calculated from the appropriate 
transfer function model presented in Table 5.6. This 
Comparison as sown invrigures: 0.2.0.5, 70.5, (D-3. 0.6, 
D.9, D.11, D.12 and. DeiS<«wf AppendtxSDwaThe values of the 


parameters determined for the more complex polynomial 
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TAB Ee 5. 6 NORMALIZED TRANSFER FUNCTION MODELS 


Test Designation 


Positive Composition 


Feed + 


Feed B 


Reflux + 


Reflux B 


Steam + 


Steam B 


Negative Composition 


Feed - 


Model 
Responses 
-6.05s 
Oe Mas 
F 2/.0S5 + 1 
x -4.05 
be TN ach eed 
5 H2t5sat d 
-0.8s 
Padi Maree 
R Lose 4 
-8.0s 
abpitade of® 
R 11.05 + 1 
-2.55 
a eS 
S 20S + 4 
-4.05s 
X @ 
weDeene 
S 2nadS) tel 
Response 
, -8.05s 
jotthem ined ai 
fF 10.05 + 1 
X draly- 35 
beDed 
R T5s0s +1 
Xp r SiS 
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20.8- 
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f + 26.%5 


20.4- 
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if + 22.81 


28.0- 
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transfer function representations of the column's dynamic 
response, using the method of Levy (63), are given in 


babies Do lie Dalz. .D..1Su0n -Appendix —D. 
5.92 Discussion 


The process gains and the parameters of a first 
order lag plus time delay transfer function, describing 
the column dynamic responses, were determined from inde- 
pendent tests. The nonlinear behavior of the column 
responses necessitated the determination of the process 
gains as a function of the magnitude of the input dis- 
turbances, and the determination of the model parameters 
as a function of the input disturbance direction. Bor- 
nard et al (14) have also experienced the similar non- 
linear behavior of the process gains. They also con- 
sidered that despite the nonlinear behavior of the column, 
the model parameters of their process model would not 


vary significantly with the magnitude of the disturbances. 


The dynamic parameters determined directly from the 
time domain data, were used to determine an average model 
describing the column dynamics, rather than those calcu- 
lated from the frequency response data. These dynamic 
parameters were considered the best available, since they 


were determined directly from the time domain data 
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rather than from a transformed data set (ie frequency 
response data). The difference between the time delays 
determined from ine two frequency response methods is 
probably due to the insensitivity of the simple search 
technique used to find the 'best' time delay incorpor- 
ated into the method of Levy. The quantitative compari- 
son of time domain and frequency domain curve fitting 


procedures has been left for a future project. 


In general, the gains calculated from the pulse 
testing procedure do not agree with aide determined 
from the transient response analysis. This is because 
the pulse does not excite the nonlinearities to the 
same degree as a step of the same magnitude. Figure 
5.4 represents the response of a hypothetical system to 
both a step and a pulse of the same magnitude. If the 
process illustrated is nonlinear, the response to the 
pulse cannot be expected to reflect the nonlinearities 
in the process that it did not experience (ie those 
nonlinearities in the output which were excited by the 
Step disturbance at times greater than the pulse dura- 
tion). Thus the gains calculated from the pulse test 


would be expected to be biased. 


A comparison of the process gains and the dynamic 
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response 
response 


disturbance 
disturbance 


time time 


FIGURE 5.4 RESPONSE OF SYSTEM TO A STEP 
AND PULSE OF SAME MAGNITUDE 


parameters of the approximate models indicates that both 
product composition responses are much more sensitive 

to the steam flow rate than to either the reflux or feed 
flow rates. The product composition responses are least 


sensitive to the feed flow rate. 


Since both feed and reflux flow rate disturbances 
essentially alter the internal reflux distribution in the 
column, it would be expected that their relative effects 
would be of the same order of magnitude, while those for 
steam flow rate disturbances would be expected to be much 
greater, due to the changes in the vapour boil up rate. 


The response of the column to vapour flow rate changes is 
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much more rapid than that due to liquid flow rate changes. 


The dynamic response of the overhead product compo- 
Sition to feed flow rate resembles that of a feed compo- 
Sition change. Table 5.7 compares the parameters of the 
approximate dynamic model for each of these disturbances. 
The values determined for the time delay and time constant 
are of the Same onder of Naqnttude.  Thts asuggests that 
the overhead product composition responds at a rate rela- 
tive to composition changes, while the bottom product 
composition, responds a@t¥a rate lrelative to-Riquid flow rate 
changes. This would indicate why the bottom product 
compositiom response, is much more sensitive to the feed 
flow rate disturbances than is the overhead product compo- 
sition response. These results tend to agree with the 
discussion presented by Wahl and Harriott (134). The 
relative speed of response, slowest for composition distur- 
bances, fastest for vapour flow rate disturbances, is 
also reflected in the order of magnitude of the time 
delays associated with the response of the overhead and 


bottom product flow rates given in Table 5.8. 


Using the complex curve fitting technique of Levy 
(63), the fit of the frequency response data can be im- 
proved using a more complex polynomial transfer function 


representation. A first order lead/second order lag 


‘ 


f 


.zapnsda ads? wort sit 8 


wn62 (+ bas yalob omry ans 107T bontwieseb eon 


7 
uy 


-ogmo2 tauborq beslieve sa to senogest 
3 $ — i 
-oqfios best 6 to tady eoldmeeo sbilitedihc' . 
sdt to 2eyetoumsyed oft 2esveqmoo VY. @ lds — 
aonedxutetb s25ne to doe9 VOT febom omen b at: 
a 
it ef2one 2td7T .oebwitnpsn to YSbYO Siike an 
"ee —* 
snoqesy nots feoqmos FoueNG baat 
sit aftdw ,esonseds fomreoqmod 
to1.6 +6 ebnogesy nots teog 
beatin te a 
ca . i {oa 
stsatbnt biwow eft s.esene 
afd of . 2 siom daum ef eenogest nots keoe 
sb | rave add ef nada 290ngdzut2bb Soe 
: a3 2tfuesy seed? ..s2n0g25" yam 
ryveH bose (heW yd bodasesrg notea. 


49 
of2 .sznoqea% To besa2 Sve 


\ Prcold 


naa 
3 spyo sAt nf besasftenr 0: 


jolt yvoqsy yo? tesdee? . 260m 
mm 


smodeay oAt ttw boteiooees . a 


deT at novip 2eat61 wolt touboxg nod 


= 


& 


| +0 suptndost onfssFt sviw> xsignos sag mab 
st6b az2neqg2ee7 yonsues7? sft. To st ad 
joftonut rat2ne1s Fatmonylog xsiqmoo stom 6 anteuts 


sel a9ob%0 bacase\bes! yebvo tettt A .nordes 


PAB LR pou COMPARISON OF APPROXIMATE MODEL PARAMETERS 
DETERMINED FOR FEED FLOW RATE AND FEED 
COMPOSITION DISTURBANCES 
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transfer function representation appears to give the 

most significant improvement in the fit above that obtained 
with the simple first order lag plus time delay model. 

The use of additional terms gives only a marginal improve- 


ment in the fit obtained. 
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CHAPTER 6 
FEEDFORWARD CONTROLLER MODEL 


6.1 Introduction 


The feedforward controller model has been calcu- 
lated using the method outlined by Bollinger and Lamb 
(12), with the plant transfer matrix described by the 
open loop process transfer function models, determined 
from the transient response and pulse testing program. 
These models were also used to estimate the parameters 


required for the experimental feedforward controllers. 
6.2 Determination of the Feedforward Gain 
6.2.1 Procedure 

The feedforward controller gain has been estimated 
using the Folvewing two methods: 


a) substitution of the open loop process 


gainssrallustrated in'Figurese5:7 andusezs 


into Equation (3.23) 
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state of the smoothed feedforward 

operating conditions (reflux flow rate 

as a function of the feed flow rate 

at constant overhead product composition). 
dR 


ia Dain yak oo (Ga.23) 
ff Redes 


The values calculated from the above methods were 


also compared with the gains determined during the actual 


on-line experimental control studies. 
6.2.2 Results 


The feedforward controller gains calculated from 
the various methods are displayed in Figure 6.1. The 
large amount of scatter in the gain calculated from the 
unsmoothed feedforward operating data amplifies the 
conclusions reached in the previous chapter, that the 
gains Should be calculated from the smoothed operating 
data rather than the untreated experimental data. The 
feedforward operating data, illustrated in Figure E.1 
of Appendix E, appear to be best fit by a straight line 
of the form 

Rel = ese 6a 005500 F (Gn 3)) 
from which the feedforward gain is a constant 


Kees ee 0.50 (6.4) 
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This value has also been displayed on Figure 6.1. The 
feedforward controller gains, calculated according to 
Equation (6.1), using the open loop process gains, 


Kyp and Kyp, have also been included in Figure 6.1, 


The two point feedforward controller gains required 
to maintain both the overhead and bottoms product compo- 
sition constant by manipulation of the reflux and steam 
flow rates, were also evaluated by substituting the 
values of the open loop process gains at the reference 


steady state into Equation (3.24) as 


we |e (6 


This controller, unfortunately, has not been evaluated. 


6.3 Characterization of the Controller Dynamic Behavior 


6.3.1 Procedure 


The dynamic feedforward controller is calculated by 
substitution of the frequency response representation of 
the open loop process models, displayed in Figures D.2, 
Ure. pre and U.taeor Appendix LD Ineo Equation, (3,2¢3.)< 
These calculations were performed using the computer pro- 
gram FFCI1, which is outlined in Appendix G. The dynamic 


feedforward controller model was also calculated using 
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the frequency response representations of the approximate 


open loop transfer functions listed in Table 5.6. 


Estimates of the dynamic parameters of simple time 
delay and first order time lag transfer functions were 
made, based on these calculated feedforward control 
frequency response representations. The parameters chosen 
during the later experimental evaluation were based on 


these estimates. 


The time delay was determined directly from the phase 
angle of the frequency response representation as 


T = (6.6) 
57 3 0) 


The time constant was determined using the complex curve 


fitting procedure of Levy. 
632 2 Y Resiteltis 


The frequency response representations, calculated 
according to Equation (3.23), using the open loop frequency 
representations, are displayed in Figures E.2 and E.3 of 
Appendix E. Also shown on these figures are the frequency 
response representations of the feedforward controller, 
calculated using the frequency response representations 


of the approximate models. These frequency responses 
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illustrate, qualitatively, the characteristics exhibited 
by the feedforward controller. The values of the time 
delay and the time constant, presented in Tables 6.1, 
6.2 and 6.3, both appear to be less than 12 minutes for 


the existing column. 
6.4 Approximate Feedforward Controller Models 


The parameters of other higher order polynomial 
transfer function representations of the feedforward 


controller of the form 


. 
: Ga? 
Grp —————-— (O27) 


are presented in Tables E.4 and E.5 of Appendix E. 
Ov9 Discussion 


The results presented in Figure 6.1 provide addi- 
tional evidence of the nonlinear behavior of the column. 
If the column actually responded in a linear manner, the 
feedforward controller gains should be those calculated 
using the open loop process gains, Kye and Kyps However, 
the superposition theorem does not appear to be applicable 


in this case, as witnessed by the lack of agreement between 


the gains calculated from the open loop process gains and 
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TABLE 6.1 TIME DELAY DETERMINED FROM THE FEEDFORWARD 
CONTROLLER FREQUENCY RESPONSE DATA FOR 
INCREASING DISTURBANCES IN THE FEED FLOW RATE 


Test Designation 
Frequency Reflux21] Reflux21 Reflux21 Reflux24 Reflux24 
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Feed20 Feed23 Feed24 Feed23 Feed24 
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TABLE 6.2 TIME DELAY DETERMINED FROM THE FEEDFORWARD 
CONTROLLER FREQUENCY RESPONSE DATA FOR 
DECREASING DISTURBANCES IN THE FEED FLOW RATE 


Test Designation 
Frequency Reflux20 Reflux22 Reflux23 Reflux23 
Feed21 Feed22 Feed2] Feed22 
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TABLE 6.3 TIME CONSTANT DETERMINED FROM THE FEED- 
FORWARD CONTROLLER FREQUENCY RESPONSE 
DATA FOR BOTH INCREASING AND DECREASING 
DISTURBANCES IN THE FEED FLOW RATE 
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that determined from the smoothed feedforward operating 
conditions. Reasonable agreement was only obtained in the 
region of the reference steady state. This would be ex- 
pected, since it is generally accepted that most nonlinear 
systems can be approximated in a linear manner for small 
deviations from the reference steady state. It is this 
fact which allows the determination of the two point 


feedforward controller given in Equation (6.5) to be made. 


Figure 6.1 has illustrated a rather unexpected result. 
Despite the nonlinearities exhibited by the open loop 
responses, the feedforward operating conditions are linear 
about the reference steady state, resulting in a constant 
feedforward gain, which is applicable over the range 
of feed flow rate disturbances studied. These results have 
been verified by the tuned feedforward gains measured 
during the evaluation of the on-line closed loop control 
of the column reported in the next chapter. These gains, 
Shown in Figure 6.1, varied only slightly from this con- 
stant value. This result is unexpected, since Luyben (69) 
has pointed out that the feedforward gain would only be 
constant at different feed flow rates, if both overhead and 
bottoms compositions were maintained constant. Clearly, 
additional work is required to explore this type of beha- 
vior that results, no doubt, because of the inherent non- 


linearities in the system. 
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In light of taking into account the nonlinearity 
of the column in calculating the feedforward gain, it 
would appear contradictory to now calculate the feedfor- 
ward dynamics by neglecting these nonlinearities. This 
was assumed possible by taking into account the following 
three considerations. First, as described by Bornard (15), 
that although the process gains exhibit a strong depen- 
dency on the operating conditions, the process dynamics 
would not be expected to possess such a strong dependency. 
Secondly, it is of prime importance to obtain the best 
value possible for the feedforward gain, while larger in- 
accuracies in the dynamic model parameters are tolerable. 
Thirdly, it is hoped that, despite the nonlinear effects, 
the same general form of an approximate model would prove 
adequate. The nonlinearity of the column would only re- 
quire that the parameters evaluated be fine tuned to the 


process. 


Despite the large amount of scatter evident in Figures 
E.2 and E.3, the feedforward controllers, calculated from 
the open loop process models, all tend to exhibit similar 
characteristics, while those calculated using the approxi- 
mate transfer function models exhibit markedly different 
properties, especially at the higher frequencies. This 
indicates that, even though the simple first order plus 


time delay transfer function models adequately predict 
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the open loop transient response of the column, they do 
notiappedr *suttirciently complex .to use In thewca culation 


of the feedforward controller. 


The intersection of the high frequency asymptote of 
the calculated feedforward models with the zero decibel 
low frequency asymptote occurs in the neighbourhood of 0.05 
to 0.15 radians/minute. Since the open loop frequency re- 
sponse values were valid only up to frequencies of 0.2 to 
0.3 radians/minute, the feedforward controller frequency 
response values would be expected to have a similar maxi- 
NM eCUL Otte Trequency.a abnus, che: Majority or Che dyna ¢ 
information concerning the feedforward controller has been 
obscured. In order to recover this information, the open 
loop frequency response of the column must be extended to 


higher frequencies than were obtained during this study. 


ExXaminaciron Of figures Eb. andsk.3 reveals that tive 
main characteristic exhibited by most of the feedforward 
models is the lead function required to fit the resonant 
peaks OCCUTYINg Inthe’ amplsTUderrat1 oVF tin! Tac thowsthin 
the narrow range of valid amplitude ratio values obtained, 
the feedforward controller required to compensate for de- 
creasing feed flow rate disturbances appears to be com- 
pletely described with only a lead function. This ob- 


servation, though, seems to contradict the observed phase 
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angle, which decreases over the frequency range. These 
seemingly contradictory observations could be explained 
if it is assumed that the feedforward controller can be 
approximated with models of the form 


~Tps 


Gee = e (t5s+1) (6.8) 

C a eal T7>T (6.9) 

ff TS tal 

G < 1 c=1..0 (6.10) 
a 


t*s* + 27zcstl 


Since the magnitude of the resonant peaks exhibited 
by the amplitude ratio are relatively small, and the 
phase angle generally decreases over the range of fre- 
quencies studied, it is also possible to use very simple 
teansieyr LUNnETIONS of athe form 


S = D 6.11 
ff ; par 


(6212)) 


Oe: TSt+] 


as approximations to the feedforward controller model. 
The range of the dynamic parameters presented in Tables 
6.1, 6.2, E.4 and E.5 are probably caused by neglecting 
the nonlinearity of the column, and treating it as if it 


were a linear system. 
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Since no specific criteria describing the "goodness 
of fit" are readily available, the choice of which model 
to use is still a subjective decision, based on considera- 
ElOnS  SUCh..as: 

a) is the model realistic? 

b) is the model easily implemented? 

c) is it worthwhile implementing a 

more complex model when a much 
simpler model will perform almost as 


well? 
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CHAPTER 7 
EXPERIMENTAL CONTROL STUDIES 


7.1 4x.Introduction 


In order to make a valid assessment of the effective- 
ness of feedforward control, the performance of a system 
using conventional feedback control was initially deter- 
mined as a base case. A number of previous workers (69, 
91,132) have compared the feedforward response to the 
open loop response of the column. This comparison will 
only indicate if feedforward control is preferable over 
no control at all. No evaluation was made as to whether 
feedforward control can improve the column performance 
beyond that achieved at present using conventional control 


schemes, such as feedback control. 


During this project therefore, the effectiveness of 
the various feedforward controllers will be compared to 
that of a conventional feedback controller, tuned to give 
the best response based on a minimum total deviation from 
the set point, as measured by the integral of the absolute 
error (IAE), a minimum overshoot and a reasonable settling 


time. 


The feedforward models implemented by the computer 
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were purposely kept simple in nature, so that they could 
easily be implemented without a computer, should the re- 
sults justify their usefulness. The three models chosen 
illustrate the simplest possible approximations to the 


frequency response results presented in Chapter 6. 


= 2] 
Gee Kee (7.1) 
G = Je ktaieu 1° fn 
ff ff 
K 
aes oe ee (7.3) 
££ tTStl 


These simple models are easily tuned to the process, 
since at most, only one dynamic parameter must be adjusted, 


in addition to the feedforward gain. 


The combined feedforward-feedback control schemes 
were implemented using the parameters which produced the 
most effective control behavior, as determined for each 
feedforward and feedback control implemented independently. 
The best feedforward-feedback controller was evaluated fur- 
ther to demonstrate its improved performance to a prolonged 


series of disturbances. 


The transient response for each individual test con- 
ducted during this project are included in Appendix F. 
Portions of these figures have been reproduced as composite 


figures comparing the effectiveness of the various control 
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schemes. 
7.2 Procedure 


Proportional plus integral feedback control was im- 
plemented using both a conventional pneumatic analog 
controller and the DDC algorithm available on the IBM- 
1800 digital computer. The controller settings were 
determined using the ultimate sensitivity method of 
Ziegler and Nichols (27). These controller constants 
were also compared with those calculated based on the 
process reaction curve of Cohen and Coon (27). The 
constants determined using the Ziegler-Nichols method 
were initially applied to the feedback controller and 


subsequently tuned to give an acceptable response. 


The gain feedforward controller was initially im- 
plemented using the loop record structure, outlined in 
Figure 3.3, with the gain, defined from Figure 6.2, set 
to 

K = - 0.5 (7.4) 


This value was then subsequently tuned to eliminate as 
much offset as possible, in the response of the overhead 
product composition due to both increasing and decreasing 
disturbances in the feed flow rate. Once the gain was 


determined, the other dynamically compensated feedforward 
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controllers, defined by Equations (7.2) and (7.3), were 
implemented using the loop record structure outlined in 
Figures 3.4 and 3.5. The values of the dynamic para- 


meters, Tp and t, were varied over the range 


Ose (te. .,t) re hE5)) 


D? 
defined in Tables 6.1 to 6.3 in order to determine the 

best value for each term. The feedforward-feedback con- 
trol schemes were implemented using the loop record struc- 
ture defined in Figures 3.6 and 3.7. The feedforward gain 
and feedback control parameters were chosen as those 

'best' values, based on the previous feedforward control 
tests. The value of the time constant was purposely 

Chosen not to be the 'best' determined in the previous 
tests, in order to show the dramatic improvement in the 
control effectiveness by the addition of the feedback trim. 
Once the evaluation of the various control schemes was com- 
pleted, it was decided to compare the effectiveness of the 
best gain plus time lag feedforward-feedback control with 
that of the feedback control alone. This comparison used 

a series of frequently changing step disturbances in the 
feed flow rate in various directions over a period of 
approximately six hours. The feedback controller para- 
meters were set to the tuned values determined previously, 
while the feedforward controller parameters were set to the 


average of the 'best' values that were determined for both 
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the increasing and decreasing feed flow rate disturbances. 


The experimental testing method follows closely that 
described previously during the open loop process model 
determination. When the desired controller had been im- 
plemented, the proper controller constants entered and 
the column was operating at its reference steady state, 
the data acquisition was begun. Data was initially col- 
lected over a period of 10-15 minutes in order to estab- 
lish the initial steady state conditions. These condi- 
tions were also recorded using the material and energy 
balance programs DASS/DATAC/BALNC. After sufficient 
steady state data had been collected, a single step dis- 
turbance in the feed flow rate was introduced. The 
transient response of the system was accumulated until 


the new steady state was attained, (generally within 
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about two hours), after which the new steady state material 


balance was reported. The transient response data were 
transferred from the system disk files and punched onto 
cards and listed by the printer using the program GTDAT. 
When the data from the previous test had been completely 
removed from the system files, a new test could be ini- 
tiated, which would return the feed flow rate to its 
reference steady state value. The documentation of the 
Steady state material balances and the collection of the 
transient response data proceeded as outlined previously. 


The transient response data available on cards, could now 
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be plotted offline using the plotting facilities avail- 
able on the IBM-1800, to examine the effectiveness of the 


control scheme just tested. 


The evaluation of the effectiveness of the various 
control schemes was based on the calculated integral of the 
absolute error (IAE) of the transient response of the con- 
trolled variable (overhead product composition) using the 


following expression: 
TAEs Gee aon Xu Sex ene (74) 


This expression calculates the IAE for the first 93.3 
minutes of the transient response after the initial en- 


trance of the disturbance into the column. 


The magnitude of the step disturbances considered 
during the remainder of this chapter is either +0.45 or 


-0.59 1b/min from the reference value of 2.46 1b/min. 
v3) pip eeaback =Control 


Table 7.1 gives the feedback controller parameters, 
(proportional and integral constants) determined using 
the Ziegler-Nichols method (27), Cohen-Coon method (27) 
and the final tuned values. The Ziegler-Nichols parameters 
were determined based on the gain at which the process 


Started to exhibit an oscillating behavior, while the 
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TAGE 70 FEEDBACK CONTROLLER SETTINGS DETERMINED 
BY VARIOUS METHODS 


Digital 
Analog Digital in same units 
Method as Analog 
PB (%) t, (min) KC KI PB (%) tT, (min) 
Ziegler- -100. 5:0 -1.25 eQer09i - 80. 6.0 
Nichols 
Cohen- -112. SA) =0. 90; 05156 -111. on 
Coon 
Tuned -200. 920 027.5; 92 02.035 aos eo) meee el a ad 6) 


Cohen-Coon parameters were determined from the parameters 

of the approximate open loop model of the column given in 
Figure 5.6. These methods were used to determine. the con- 
troller parameters for both the conventional analog control- 
ler and the digitally implemented controller. Since the 
controller constants used for both controllers are general- 
ly quoted~in different units, Table 7.1 also contains. the 
control lex-constants-for—the+di-g1tal+control-er—converted 


to thessameaunits’'as “the’analog controler. 


The transient response of the column under analog 
feedback control, using the Ziegler-Nichols controller 
constants, 1 1) Mustrdted in Figure 7eane while the transient 
response of the digital feedback controller is illustrated 
in Figure 7.2. The use of the Ziegler-Nichols controller 


settings produced a very good controlled response with 
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INCREASING FEED FLOW RATE 
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FIGURE 7.1 ANALOG FEEDBACK CONTROL: 


VARIATION OF CONTROL EFFECTIVENESS USING 
ZIEGLER-NICHOLS AND TUNED CONTROLLER PARAMETERS 
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INCREASING FEED FLOW RATE 
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FIGURE 7.2 DIGITAL FEEDBACK CONTROL: 


VARIATION OF CONTROL EFFECTIVENESS USING 
ZIEGLER-NICHOLS AND TUNED CONTROLLER PARAMETERS 


141 


[St 


i = ys 
Ae 


as 
7 aie a’ lee - ne 
- i> oe 


TPAD Teena 
Pa mtd ae : : 
2 S3Al 


a 


efodotit-<afents 
-_——_—-- - —— pee ete ent 


I RATRAY_ (LLEAMD 


q HO 1333 SKLZARRO IG 
Sass [3 LEV 108-82 >T2aT 
} on cS. f-"5y4 
- re0.0 = of 


benui 


a 


Sx bay BS 


-JOHTHOD BOARGGSIS JATISIG 
SHL2U 22IMIVITIIAIWT JONTAGD W NOLTATRAY 
‘ASAD SILTJONTAOD OHSABT GHA SIJOHSIN-AIJe3sTh 


142 


reference to only the small values of the IAE obtained. 
However, the response appears to contain very lightly 
damped oscillations. Although the amplitude of these os- 
cillations was very small, ie of the order of 0.1% by 
weight methanol, this type of response would likely be 
unacceptable for industrial control purposes. The con- 
trol parameters were tuned by relaxing the gain and in- 
tegral constants somewhat, in order to produce a response 
in which the oscillations appear to be damped to a 
larger extent. The transient response of the column 
under analog feedback control using these tuned control- 
ler settings are also shown in Figure 7.7. while those 


Under aigital Teedback contrelare: shown in Figure /.2’. 


The response of the controlled column using these 
tuned controller settings generally exhibited larger 
total deviations from the set point, a larger max imum 
deviation and a greater amount of oscillation damping 
than does the response obtained using the Ziegler- 


Nichols settings. 


Due to the nonlinear behavior of the overhead product 
composition to reflux flow rate disturbances, compromise 
controller settings were selected. The controller para- 
meter settings shown inerable 7.1. tend to give tighter 


control for input disturbances, causing increasing ‘rather 
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than decreasing deviations in the overhead product compo- 


Sule Odie 


The effectiveness of the feedback control scheme, as 
measured by the IAE during these tests, is taken as the 
basis of comparison for evaluating the effectiveness of 


the various feedforward control configurations. 


7.4 Gain Feedforward Control 


Initially, the feedforward controller gain was set 
to 0.5, as suggested in Figure 6.2. The gain was then 
tuned in order to eliminate any offset resulting from 
both increasing and decreasing step disturbances in the 
feed flow rate. The tuned gains were determined to be 
0.45 and 0.51 for step disturbances of +0.45 and -0.54 
lb/min in the feed flow rate respectively. It was found 
that slight changes in these values were necessary to 
compensate for small variations in the feed composition 


over long periods of continuous operation. 


Despite the attempts made to maintain as little off- 
set as possible, many runs do contain an appreciable 
offset due to unknown causes. This illustrates the major 
weakness of a feedforward controller; its inability to 


correct for disturbances other than the ones to which it 
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has been related. This contrasts with a feedback con+ 
troller which can eventually correct for any disturbance 


as soon as the disturbance affects the controlled variable. 


Figure 7.3 illustrates the response of the column 
under gain feedforward control for both increasing and 
decreasing disturbances in the feed flow rate away from 
the reference steady state. The overhead product compo- 
sition response can. be seen to deviate from its set 
point in a direction opposite to that expected to be 
caused by the disturbance alone, and returns to the set 
point with no overshoot or oscillations. This indicates 
that the reflux corrective action has been introduced into 
the column too early, demonstrating the need for dynamic 


compensation. 


It should be noted that if the only two alternatives 
available for controlling the column were gain feedforward 
control and conventional feedback control, these tests 
would suggest that feedback control would be the preferred 
configuration. The response of the column under feedback 
control exhibits the smaller IAE value, the smaller maxi- 
mum deviation and the smaller offset in the face of immeas- 
urable disturbances of the two possible alternatives. If 
the gain feedforward controller is compared only to the 


open loop response, then it would be an acceptable control 
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FIGURE 7.3 GAIN FEEDFORWARD CONTROL: 
VARIATION OF THE GAIN WITH THE 
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7.5 Gain plus Time Delay Feedforward Control 


Tables 6.1 and 6.2 present values of the time delays, 
calculated over a range of frequencies, from the feedfor- 
ward controller frequency response data in Tables E.2 and 
E.3. These results suggest that for increasing feed flow 
rate disturbances, a time delay of between three to nine 
minutes is required, while for decreasing disturbances, a 
time delay of the order of zero to six minutes would be 
desirable. The fact that the time delay values calcu- 
lated for increasing disturbances appear reasonably con- 
stant over the frequency interval of interest, suggests 
that a time delay is a fairly good representation of the 
data, whereas the time delays calculated for decreasing 
disturbances exhibit a considerable variation over the 
same frequency interval, suggesting that it is a less 


representative model in this instance. 


Figure 7.4 compares the transient responses of the 
column using a gain plus time delay feedforward controller 
as fa function ofStne time dalay. Examination of this fi- 
GUYS indicates that gihep best value off the time delay) ts: of 
tne orger of five minutes... Examination of. Figures F.3c4d4F.35 


and F.36 reveals that this value appears equally valid for 
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both increasing and decreasing disturbances away from and 


returning to the reference steady state conditions. 


The use of a time delay dynamic compensation cannot 
be expected to totally eliminate the transient response 
in the overhead product composition, unless the feed and 
reflux dynamics differ only by a time delay. The best 
that can be expected is to adjust the time delay such that 
the area under the transient response is approximately 
equally divided on either side of the set point cross 
over, aseshown in Figure 7.4 for a time delay or 4:8 


minutes. 


The addition of time delay dynamic compensation to a 
gain feedforward controller has reduced the maximum devia- 
tion by about 75% and the IAE by approximately 75 to 90% 
from that achieved using a gain feedforward controller 
alone. The time delay feedforward controller has also 
reduced the IAE below that achieved using only a feedback 


controller. 


7.6 Gain plus Time Lag Feedforward Control 


Table 6.3 presents the time constant obtained by 


Titting the feedforward controller trequency) response: va- 
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representation using the method of Levy (63). The re- 
sults suggest a time constant ranging between four and 
ten minutes for increasing disturbances and between two 
and four minutes for decreasing disturbances in the feed 


flow rate away from the reference steady state. 


Figure 7.5 compares the overhead product composi- 
tion response of the column as a function of the time 
constant of the time lag feedforward controller, for an 
increasing step disturbance in the feed flow rate away 
Thomeeene reference steady State. Figures /.6 and’ 7.7 
present similar comparisons for a decreasing step distur- 
bance in the feed flow rate away from and an increasing 
step disturbance in the feed flow rate returning to the 


reference steady state. 


Examination of Figure 7.5 indicates that for an in- 
creasing disturbance away from the reference steady state, 
the best value of the time constant appears to be of the 
order of 6.6 minutes. In fact, the feedforward controller 
using this value of the time constant maintained the over- 
head product composition essentially constant, despite the 
approximately +20% feed flow rate disturbance. Examination 
of Figures 7.6 and /./7 Indicate ‘thatefor a decreasing feed 
flow rate disturbance away from the reference steady state 


and for an increasing feed flow rate disturbance returning 


7 
3 fe 2 


~~. oAT . (£8) yved ° ro oO 5 
pals eae . 
bos wot nsawied patensy Ins: 20 


ows neswied brs esonedyed2tb pntensy: 
ae 
bsat oAt nt 2ee0nedyus2rb valent) 


SD 


.stss2 yosste2 e2na1sts7 | 


» ” 
i 
» 


-tzoqgmoa touborg bsatisvo sat estsqmoa & 


Wi = 
o nottonut 6 26 amyloo eat Fo senoq: 


ss is 7 
- _ 
ng tot .wafforsnos bryswrotbest pst omit edd To 
pa] 


we atsy wolt best ands nat soneadtusretb qete pnt 
t.¥ bas 3.0 2e1weti .stsete ybsste gIns1stay on: 


> ONFZ6ST9ISD § 1o7 onozteques xettmte tn 

(3 

sont as bas mov? yswe sie% wot? Baar eit a 
iJ OF BRITA YNSS ts 1 wolt best saz nt sonsdaudeth 
.siste ybsote onan 
> &® 

7 


bat @.% savet? to nord sntmsx : 
adi movt Yaws sonediusetb enter 


sqq6 tnatenos emit oad To sul ev Be 


isffortnoa biswretbest odd . F567 wl .2etuntm es a to tat 


sa +h noah + 
: 7 ; 


etsi2 yvbsste sansietses sii mov? (ews sonsdtwi2tb dit 


on tnyus97 sdyus2tb ote1 wolt best patzesetont a 


iord anrmaxd sonedaw3 tetb s36% wo i 539 r0S+ sfota 


; - ae 9 ee 
99% pntess1o9b 6 tot ted sdsotbnt 1.4 bos &.¥ sail 


+ ova 
vo ant banned atem InaTteneo sure sn2 to sulsy Rie 


.instenos ylistiasees nots taoqmes tubo 


ia 


7 


e 


XO ( WT PC MECH ) XO ( ¥T PC ELH } 


XD) ( WI PC MECH ) 


TS 


INCREASING FEED FLOW RATE 


oy. TEST tug (FFC=30 CONTROLLED VARTASLE 
Ke p2-0.45 SS =S5-CS 
Ps TAE=11°59 
t =2£040 
qQ7. TES Tc RE PC=36 CONTROLLED VARIABLE 
Ke p=-0.45 SS =Sj5-9Q9 
a TAE= 1°24 
tT = 6.6 
a. TEST: _@FFC-64 CONTROLLED VARTABLE 
Ke p=-0.45 BS =S5*10 
t =10.0 TAE= 6°62 
oO 24 48 ae! 56 oat 
TIME (MIN) 
PaGuUR GE) aD GAIN PLUS TIME LAG FEEDFORWARD CONTROL: 


VARTATION OF THE CONTROL EFFECTIVENESS AS A 
FUNCTION OF THE TIME CONSTANT 


oe — | 


SLBATRAV LOT 
0.825 ce 
GE + bL=SAL 


LSATISAV CSL ae-334. 123) 
‘ 2 a > © = af. pal 
ef =SAl 3.0 « ; _ ae 


“T™ 


— 
Be 
—— ~ —————— — _ SS 
« _ —-~ -+ _ = ~~ ee 


Tyo P£-347  2T2aT acc 
ef. 0-*,44 


9.9 =TA] 0.01 oo 


= ot @p BS 
-JOATHOD GAAWAOIGSSA BAI JMIF 2u49 HIAd 


2A 224N9VIT9IIIa JOATHOD JHT WO WOTTAIAAV 
THAT2HO9 SMIT JHT 30 WOTTINUA 


151 


INVLSNOD 4WIL JHL 40 NOTLONNS 
V SW SSAN3ATLI9449 TONLNOD FHL 10 NOILVIUVA 


*TOMULNOD GYWMYOIGIIS OVI 3WIL SNId NIVO 9° JYNSI4 


(NIN) ATL 


i aR RS °S5 
‘cs 6 
- Bree 5 5 
4 
Eee oc = 
1d: E=3V1 : ey B 
60:S6= SS L$*O-=°"¥ 2 ee 
JBVINVA CSTOMIND 09-944 21831 5 
(NIN) = =3ATL (NIA) 
o2t 95 az 
as) 
ze 
= 
a 
i 5 
ZO*6t=3VI -eT= 
80°S6= SS 1S *0-=t4y Peres 
JEVISVA COTOMIND 29-944 21831 ‘65 = JEVIYVA CSTOYIND 
“oS 
cs 8 
4 
| 
e ne a) 
Ov*S =3VI Cay B 
60:S6= SS L$‘0-=77% 
FEVIVA OSTCYLNGD $-944 *1S31 “5 ~~ JAVINVA COTTON 


J1lvd MO14 0334 ONISV34I9IG 


SAIL 


ve 


99" = 
LS‘0-= 
85-34 


0°0 = 
{S “0-= 
9S-J44 


2 
dty 


*1S31 


1 
$d 


*LS31 


bh # R HW B 


Bb R R R R 


( HOGA Jd 1A) &X 


( HOGAN Od 1K) &X 


a 


i soe = ual ee 
pt > ei — nae eel 
_ a 
a he : > 


4 
¥ ei Sah - 


’ 


56. te SLBATRAY [TEL LOATMO 82-337 +1237 
ee “as de £2.0-my 44 


3.6 * 


$8-939 «T2397 
fe. O-"5 4 


#0 ¢ 


(RL BO OH 
“ 
- 


2 
Ri 
%) 

~ hee 
Mia 


in. oS) a te nr a cc 0 hee grind +e2 
o oS 2 St =o aS a 


(ili) Batt 


SATAY CLOT 68-277 :1237 see 


7 Shck= cz 6 ? 4 — 
PV -EEss Led 7 “e 
+ 6 : Corel 
& : i. , 
“SY a 
- “ 
a —_ a : 4 Te Pe ge - — eee: | 
a , va | _ ; 
: a | ". 
7 — 
s ' ity a 
~ 7 
» 
ae / f 
—— 
Oo - _- oe ele a ee ——_——A QO 
te cm o ‘ c 
sia “27 7 
i} me 18 | 
jORTHO U aAJ : 
— 2A 227H3¥ 4333 JORTHOD rT 40 
rm A } Miy i 4 


162 


INVISNOD AWIL JHL 40 NOTLONAS 
V SV SSINJATL99443 TOULNOD FHL 40 NOTLWIYYA 
>TOULNOD GUVMYOSOIIS IVI JWIL SNId NIWI 


(NIK) FAIL 


ae O2Tt 35 
‘cs 8 
“35 = 
is = = 
| 2 ee 
JEVYISYA COTQRUND 19-944 2 LS39L es oa y 
°S5 
‘cs 6 
“oS = 
O&-TFt=2vI ed i‘ B 
60-S5= 55 1s 0-=7 79 oe 
JBVYIYVA COTOIND 6S-J44 2 1S39L Np 


£°d 3undt4 


(NIN) 
az 


65 
= 
3 
e = 
Ty -OT=3vI © a yy! ci A vO*E2=3VI 
80°S5= SS LS fe =n GT -Ss= 4S 
FSVISVA CITOMIND l$-)44 24831 ‘os ~~ JEVINVA CATON 


(31LVLS AGVILS FON3Y349Y FHL OL ONINUNLIY) 
JLVY MOTS 0334 OINISVIYINI 


Sil 
er 


ve 


ey ar 
1S ‘0-=24y 


Coao 34 


ollisyipt 


B8°p = 1 


33 


{S°O-= 4 


S9-93'3 


t1S4L 


070 a 2 
ts‘o-=? Fy 


SS-944 


sESAL 


( HOSN Jd Ik) OX ( HOSN Jd Lh) &X 


( HN Od 1A) &X 


14-333 77237 2 
FE .G-%4 53 
0.0. r 


YO ( RL OC ELH } 


*JOATHOD GRANROIWISI 2 
A ZA ZZIHIVLTSIAW JOATHOD 
THATCHOD SALT 3H1 


$e 


fh, 


‘ 


a (LOAM 
as B2 
S Al 


Wi 


hj S4TT ZUs9 IAD v.% 390613 
347 


70 HOLITATAAY 
90 KOLTIAUG 


- 
| 
io 


M) 


SAT 


0 ( ab bo eH) 


GO ( 61 GC WET 


to the reference steady state, the best value of the 


time constant appears to be of the order of 4.5 minutes. 


Tables 7.2 and 7.3 compare the controlled responses 
achieved in Figures 7.6 and 7.7 respectively. The IAE 
values have been corrected to remove the bias caused by 
the final offset observed in the response prior to the 
run termination. This would tend to make the comparison 
between the various time constants more dramatic. These 
WeSuUTtS Illustrate that, if 'care 1S not exercised in the 
selection of the time constant, the effectiveness of a 
time lag feedforward controller can be degraded to a level 
below that obtainable using only a gain feedforward con- 
trorier. sUndoupted|y. a isimilar result would be revident 
when using a time delay feedforward controller, had a 
larger range of time delay values been studied. Obvious- — 
ly, too much dynamic compensation can be much more harmful 


than not enough. 


A well tuned gain plus time lag feedforward control- 
ler appears capable of reducing the IAE and the maximum 
deviation of the controlled response by at least 70% below 
those obtainable using a gain feedforward controller, and 
by at least 50% below those obtainable using a conventional 
feedback controller. A gain plus time lag feedforward 


controller and a gain plus time delay feedforward control- 
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ler, when both are properly tuned, appear to be equally 
effective in reducing both the IAE and the maximum devia- 


ETO: 
7.7 Feedforward-Feedback Control 


Figures 7.8 and 7.9 compare the transient response 
of the column obtained by adding feedback trim to the 
feedforward controller with that obtained using the 
feedforward controller alone. The best time constants 
of the gain plus time lag feedforward controller, deter- 
mined in the previous section, were 6.6 minutes for 
increasing and 4.5 minutes for decreasing feed flow 
rate disturbances away from the reference steady state. 
Despite this fact, the values of the time constants 
were chosen as 10.0 and 6.6 respectively, in order to 
illustrate the rather dramatic increase in the effective- 
ness of even a poorly tuned gain plus time lag feedfor- 
ward controller with feedback trim, over that obtained 


using a gain feedforward controller. 


The addition of feedback trim to a gain feedforward 
controller has markedly reduced the effect of the distur- 
bance on the overhead product composition response, as 


measured by the IAE; by at least 60% over that obtainable 
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using a gain feedforward controller; and by at least 20% 
over that obtainable using a conventional feedback con- 
troller. When the feedback trim is added to the gain 

plus time lag feedforward, the effect of the disturbance was 
reduced by at least 80% from that obtained using a poorly 
tuned gain plus time lag feedforward controller alone; by 
at least 20% from that obtained using a well tuned gain 
plus time lag feedforward controller alone; and by at 
least 50% over that obtained using a conventional feedback 
controller alone. The addition of feedback trim has also 
eliminated the majority of the offset due to either model 


errors or immeasurable disturbances. 


This property has been dramatically illustrated in 
Figure 7.10. During the evaluation of the gain feedforward 
controller, the feed temperature controller failed, allow- 
ing the feed temperature to decrease from 163°F to sel sie 
Creating a large offset in the overhead product composition. 
A similar test was also performed when evaluating the gain 
feedforward Contkolien wren feadback trim. A step decrease 
in the feed temperature of -10°F was introduced into the 
column, in addition to the step increase in the feed flow 
rate. The resulting response.contained no offset and a 
marked decrease in the effect of the combined disturbance. 
In fact, the IAE for the gain feedforward controller with 


feedback trim was even reduced well below that obtained 
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with the gain feedforward controller when no immeasurable 


disturbance was present. 


7.8 Comparison of Controller Responses to a Series 


of Step Disturbances 


Since the feedforward-feedback controller produced 
such a significant increase in the effectiveness of the 
control, it was decided to compare the effectiveness of 
this control scheme to that obtained using a conventional 
feedback controller, using a prolonged series of increasing 
and decreasing step disturbances in the feed flow rate 
about the reference steady state. Examination of Figure 
i. Vl indicates a reduction of approximately 702 in both 
the IAE and the maximum deviation about the set point has 
been achieved using a gain plus time lag feedforward 
controller with feedback trim in place of the feedback 
controller. The feedforward controller parameters used 
were the averages of those values reported previously as 
the best values for the increasing and decreasing distur- 


bances in the feed flow rate. 
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reported previously in Table 7.1 were also used. 
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Figures 7.12 and 7.13 compare the effectiveness of 
the various control schemes evaluated during this study. 
Examination of these figures suggests that the order of 
preference of the control schemes discussed throughout 
this project would be 

a) dynamic compensated feedforward 

control with feedback trim 

b) gain feedforward control with 

feedback trim 

c) dynamic compensated feedforward control 

d) feedback control 

e) gain feedforward control 


f) open loop response. 


Feedforward-feedback control schemes, in addition 
to improving the effectiveness of the control, also corrects 
for the effectvof any immeasurable disturbances or of any 
errors in the feedforward control model. These errors may 
include either errors in the feedforward gain or in the 
use of poorly tuned dynamic parameters. Feedforward con- 
trol alone does not exhibit this regulation capability. 
The addition of feedback trim to a feedforward controller 
actually reduces the transient response below that obtain- 


able using feedforward control alone, in addition to the 
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long term elimination of offsets. 


Both time delay and time lag dynamic compensation 
appear of about equal ability in improving the effective- 
ness of the control. However, without access to a digi- 
tal computer, a pure time delay function is quite diffi- 
cult to implement. For this reason, the time lag transfer 
function was used in all the subsequent feedforward-feed- 
back control schemes. ,Lt.isuvinteresting, to note .«that<for 
the present column, the value of the time delay and the 
time constant are both of the order of five minutes. 

These values are contained within the ranges presented 
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CHAPTER 8 
CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 


A series of approximate transfer function models, 
which adequately describe the transient response of 
the column, have been successfully developed using 
a combination of transient response and pulse test- 


ing techniques. 


The process gains should be calculated from the 
smoothed operating conditions, rather than from 
the original experimental data. This suggests 
that any attempt to adaptively tune the gain of a 
controller on-line will require a smoothing or 


filtering algorithm in the procedure. 


The calculated process gains exhibit a strong non- 
linear dependence on the magnitude of the distur- 
bance. This nonlinear behavior is typical of columns 
producing relatively pure products, due to the 
operation of the column within the pinched regions. 
The nonlinear process gains could be adequately 
represented by a least squares straight line over the 


range of interest. 
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The open loop dynamic response of the column was 
successfully represented by a simple transfer 


function model of the form 
SG alpha: side Sale 
| TS+] 
Due to the nonlinear behavior of the column, two sets 
of parameters were determined, one set each for in- 
creasing and decreasing disturbances about the refer- 
ence steady state. No attempt was made to determine 
a functional dependence between the magnitude of the 
time delay or time constant of this model and the 
magnitude of the disturbance. It should be noted 
that a similar model has been successfully used to 


describe the dynamics of large industrial columns 


(iA eile 2))2 


The response of the column was most sensitive to the 
steam flow rate, followed by reflux flow rate, feed 

flow rate and least sensitive to the feed composition. 
This indicates that the column responds most rapidly 
to vapour flow rates, followed by liquid flow rates, 


but responds very slowly for composition changes. 


The feedforward gain could not be calculated directly 
from the open loop process gains due to the nonlinear 


behavior of the column. The feedforward gain 
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was measured from the slope of the feedforward 


Operating conditions, according to 


snp (8.2) 


K poe ute 
Tet: dF Xp 


Despite the nonlinear behavior of the column, the 
feedforward operating conditions turned out to be 
linear, resulting in a constant value of the feed- 
forward gain. The feedforward gain determined 

during the control system evaluation tests varied 


only slightly from this constant value. 


The feedforward controller gain can, however, be 
calculated from the open loop process gains valid 


at the reference steady state. 


The frequency response representation of the feed- 
forward controllers calculated from the experimen- 
tally measured open loop models could not be predicted, 
within satisfactory limits, using the approximate 
open loop transfer function models. This situation 
resulted, despite the fact that the approximate 
models adequately represented the open loop transient 
response of the column. These simple models do not 
appear sufficiently complex to be used to determine 


the feedforward controller. 
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Since the range of the frequency response values 


did not extend to sufficiently high values of fre- 


quency, most of the dynamic information describing 
the feedforward controller has not been obtained. 
However, ranges of dynamic parameters of very simple 


transfer functions of the form 


-THS 
= D 
Gee e Giron) 
| ] 
TSt] 


were determined. 


The implementation of feedforward controllers in- 
corporating these models yielded a significant 
improvement in the effectiveness of the control 

of the column. The experimentally tuned values of 
the parameters of these models were included within 
the range of values determined from the frequency 
response representations. These results agree with 
those of Luyben (69,70) that a first order dynamic 


Compensation should be adequate. 


Despite the nonlinear behavior exhibited by the 
column, simple and effective feedforward and feedfor- 
ward-feedback control schemes were successfully 


implemented. 
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13) Based on the results presented in this project, the 
order of preference for the implementation of the 
various control schemes evaluated during this study 
would include: 


a) dynamic feedforward control with 
feedback trim 


b) gain feedforward control with 
feedback trim 


c) dynamic feedforward control 
d) feedback control 


e) gain feedforward control 


14) Besides eliminating offsets caused by either the 
effect of immeasurable disturbances or errors in the 
feedforward model, due to either the approximate 
nature of the model or the varying column conditions, 
the addition of a feedback trim to a feedforward 
controller actually increases the effectiveness of 


the control. 


15) Although the gain feedforward controller is not as 
effective as a conventional feedback control, the 
combination of the two greatly enhances the effective- 


ness of the overall composition control. 


16) Only a limited knowledge of the column dynamics need 


be known in order to implement a feedforward controller 
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with feedback trim. 
8.2 Recommendations 
8.2.1 Distillation Column 


Some of the flow recorders used during this project, 
including reflux, steam, overhead and bottoms product flows, 
were pperated below approximately 30% of the transmitter 
Span. This reduces the accuracy with which these flows 
Can be measured and reset. Since the column operates very 
Close to the flooding constraint, the span should be re- 
duced to improve the accuracy and repeatability of meter- 


ing these flow rates. 


The cause of the flooding appears to be undersized 
downcomers. If they could be made larger, so that the in- 
ternal reflux flow could be increased, the external flows 
could also be increased, resulting in shorter time con- 
stants for the dynamic response of the column, thus decreas- 


ing.the time required» for the transient response tests. 


During the present study, the column was operated for 
long periods of time, often unattended. Any failure of any 
of the utilities, such as steam, cooling water, instrument 
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to be lost, pressures to rise, pumps to be operated dry, 

or many other potentially hazardous conditions, which 

could result in damage to the equipment. Since the com- 
puter has access to all the process measurements, it would 
be advantageous to have an operations monitor type of 
program , similar to that available on the evaporator (54), 
which could automatically shut down the column when symp- 


toms of any hazardous conditions are detected. 


8.2.2 Distillation Column Model 


Sufficient experimental data has been presented to 
form the basis of a fundamental model which could be 
developed. This model could be used in future studies to 
aid in the study and development of other advanced control 
schemes. A satisfactory model should permit determination 
of the column response about any other column operating 
condition, rather than require the collection of new ex- 
perimental data, since an experimental program is general- 
ly expensive, time consuming and inconvenient when applied 


to an industrial column. 


8.2.3 Distillation Column Approximate Model 


Further pulse tests are required to extend the fre- 


quency response of the open loop responses, so that the 
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feedforward controller frequency response can be calcu- 
lated to frequencies above 0.2 radian/minute. The fre- 
quency response of the feedforward controller was only 
valid up to this frequency, due to the pulse testing 
limitation imposed by the measured open loop frequency 
response. Since the time constant of the feedforward 
controller is of the order of five minutes, the break 
Bcequgencyseqfathermodel¢ersyin «the vicinity of 0.2 radians/ 
minute. Thus most of the high frequency characteristics 


of the model have been lost. 


Further study should be required to determine the 
effect on the pulse test results of the nonlinear behav- 
ior of the column, and the causes of the abnormally high 
values of the normalized frequency content at the limit 
of the valid frequency response data. It is also possible 
that other stochastic testing techniques, such as Kalman 
filtering, pseudo-random binary inputs or an auto and 
cross correlation analysis of process noise, may provide 
better estimates of the dynamic model of the column than 
pulse testing. The existing column, with its data acquisi- 
tion and control capabilities, is an ideal vehicle for the 


evaluation of these alternate methods. 


8.2.4 Transient Response versus Frequency Response 


The determination of the parameters of an approximate 
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model using data in the time or frequency domain requires 
further study. It is necessary to determine criteria 

to decide which of the many methods of fitting the data 
Should be used, and in the case of complex models, to 
determine when a sufficient number of terms have been in- 


cluded. 


8.2.5 Feedforward Control 


The feedforward controllers designed during this 
study are based on conventional continuous methods. The 
Sample time chosen when implementing the controllers on 
the computer was small relative to the time constants 
associated with the process. This, however, does not 
make the best use of the computer. A digital feedforward 
model, based on sampled data techniques, should be designed, 
evaluated experimentally and compared with the control 
obtained using the continuous model. Further study would 
then be required to find the best value of the sampling 


interval for this model. 


Dynamic feedforward control, incorporating only 
simple models such as a time delay or a time lag, has 
reduced the amount of deviation in the overhead product 
composition response in the order of 75% of that obtained 


using conventional feedback control alone. Further study 
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could determine if the use of more complex models, such as 
a first order lead/lag (115), first order lag plus time 
delay (14), a second order lag (66), a second order lag 
plus time delay (66) or a first order lead/second order 
lag (82) would give sufficiently increased effectiveness 


of control to justify their implementation. 


Under certain circumstances, it may be advantageous 
to control both overhead and bottoms product composition. 
Various suggestions have been made (1,9,73,102,103,108) 
to implement non-interacting control, including a few 
based on feedforward concepts (1,73). Further experimen- 
tal study is required to determine the feasibility and 
parameters required to implement the two-point feedforward 


control system. 


8.2.6 Material Balance Control Scheme 


The column material balance is currently controlled 
indirectly, with both overhead and bottoms product flow 
rates on level control. An alternate control scheme, re- 
commended by Shinskey (115) and Nisenfeld (89), controls 
the material balance directly by adjusting either the over- 
head or bottoms flow rate directly, rather than either 
reflux or steam flow rates. Since each control philosophy 


appears to have obtained about equal acceptance, based on 
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literature, references, it would seem worthwhile to experi- 
mentally study the direct method of material balance 
control and compare the results with those presented during 
the current project. This comparison could include the 
column dynamics, the feedforward controllers and the 
effectiveness of the feedforward control schemes. A 

study of this nature may help resolve the controversy 
concerning which material balance control scheme is more 
desirable, or under what circumstances each control scheme 


might be preferred. 
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NOMENCLATURE 


defined by Equation (3.5a) 
analog to digital converter 
amplifier 

amplitude ratio, decibels 
bottoms product flow rate, 1b/min 
defined by Equation (3.5b) 

bias in the DDC loop record 


defined in Equation (3.5c) 


vector of unknown transfer function 
coefficients 


current Output Station 
overhead product flow rate, 1b/min 
defined by Equation (3. 5d) 


m 


(3 ne denominator defined by Equation 
3.9 


degrees of freedom 


error function defined by Equation (3.13) 


error function in frequency domain 
defined by Equation (3.9) 


feed flow rate, Ilb/min 
flow controller 


exponential filter constant in the DDC 
algorithm 


set point of feed flow rate controller 
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steady state value of feed flow rate, 
lb/min 


function defined to fit the experimental 
frequency response data exactly, as 
defined by Equation (3.8) 


flow recorder 
flow recorder/controller 


general transfer function in the 
frequency domain 


general transfer function in the 
Laplace domain 


general feedback controller equation 
in the time domain 


transfer function of the feed flow rate 
measurement transmitter 


transfer function relating the overhead 
product composition to the ‘'jth' input 
variable 


transfer function relating the bottom 
product composition to the ‘jth' input 
variable 
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where the ‘ith' output variable refers to 


B3.a)bottom product flow rate 


D - overhead product flow rate 
QO - any other output variable of 
interest 


and the 'jth' input variable refers to 
F - feed flow rate 


L - any other input variable of 
interest 


R - reflux flow rate 
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S - steam flow rate 
seat feed temperature 
V - boilup flow rate 


Xe - feed composition 


closed loop transfer function of the 
reflux flow rate to a set point change 


transfer function of the reflux flow 
rate controller 


transfer function of the reflux flow 
rate measurement transmitter 


transfer function of the reflux flow 
rate control valve 


transfer function of the overhead 
product composition controller 


transfer function of the overhead 
product composition measurement trans- 
mitter 


transfer function of the feedforward 
controller 


water 


integral of the absolute error 

ues function defined by Equation 
3:8) 

current to, pressure itransducer 

gain of a general transfer function 


feedback controller gain 


feedback controller gain in DDC loop 
record 


gain of the feed flow rate measurement 
transmitter 


integral controller constant for digital 
feedback controller 
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integral controller constant in DDC 
loop record 


gain of the transfer function relating 
the composition of plate N to the ‘jth' 
input variable 


gain of the transfer function relating 


the overhead product composition to the 
'jth' input variable 


gain of the transfer function relating 
the bottoms product composition to the 
'jth' input variable 


gain of the transfer function relating 
the 'ith' output variable to the ‘jth’ 
input variable 

where the 'ith' output variable refers 
to 


B - bottoms flow rate 
D - overhead flow rate 
QO - any other output variable 


and the 'jth' input variable refers to 


F - feed flow rate 


L - any other input variable 
R - reflux flow rate 
S - steam flow rate 


(bes feed temperature 
V - boilup flow rate 


Xe - feed composition 


any other input variable not otherwise 
defined, ie Tp» Ths Tes To 


Laplace transform operator 
level controller 


level indicating controller 
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multiplexer 

measurement in DDC loop record 
general plate number 

number of data points collected 


total number of stages including the 
reboiler 


ea numerator defined by Equation 
Sak 


any other output variable not otherwise 
defined, ie Xw > P 


output in DDC loop record 
pressure, inch H50 


proportional band of feedback con- 
troller, % 


pressure controller 

pressure to current transducer 
pressure indicating controller 
pulse output 

process operators console 

reflux flow rate, Ib/min 

set point of reflux flow controller 


steady state value of the reflux 
flow rate, 1lb/min 


ete function defined by Equation 
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steam flow rate 
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steady state value 
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finite integration limit introduced 
intosequation (3.5) 


atmospheric temperature, deg F 

cooling water inlet temperature, deg F 
feed temperature, deg F 

reflux temperature, deg F 

temperature recorder 

temperature recorder/controller 

vapour boilup rate, 1b/min 

vector defined by Equation (3.15) 
matrix defined by Equation (3.15) 


bottom product composition, weight 
% methanol 


steady state value of the bottom 
product composition, weight % methanol 


overhead product composition, weight 
% methanol 


set point of the overhead product 
Composition controller 


steady state value of the overhead 
product composition, weight % methanol 


feed composition, weight % methanol 


composition on plate N, weight % 
methanol 


coefficients of numerator in polynomial 
i obabllaeh function defined in Equation 
6.8) 


coefficients of denominator in polynomial 
i ee function defined in Equation 
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e - exponential operator, 2.7183... 

i aN ea 

m - degree of denominator in polynomial 
transfer function defined in Equation 
(6.8) : 

n - degree of numerator in polynomial 
transfer function defined in Equation 
(6.8 

r(w) - real part of complex function defined 
by Equation (3. ii) 

S - Laplace transform operator, anh | 

s(w) - imaginary part of complex function 
defined by Equation (3.11) 

t - time, minutes 

ed (br) - transient response of an input 
variable 

at) - transient response of an output 
variable 


iii) Greek 


1 - defined in Equation (3.7a) 

8 - defined by Equation (3.7b) 

Y - defined by Equation (3.7d) 

AB - variation of the bottoms flow rate 
about the reference steady state value, 
lb/min 

AD - variation of the overhead flow rate 
about the reference steady state value, 
1b/min 

AF - variation of the feed flow rate about 


the reference steady state value, Ilb/min 


AR - variation of the reflux flow rate about 
the reference steady state value, I]b/min 
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variation of the steam flow rate about 
the reference steady state value, 1]b/min 


variation of the feed temperature 
about the reference steady state value 


variation of the bottoms product 
composition about the reference steady 
state value, weight % methanol 

variation of the overhead product 
composition about the reference steady 
State value, weight % methanol 

sampling interval, minutes 

error in the process gain determinations 
for feed flow rate disturbances, defined 
by Equation (5.1a) 


Component material balance error of 
closure for water 


component material balance error of 
closure for methanol 


overall material balance error of 
closure for flow rates 


error in process gain determinations 
for reflux flow rate disturbances, 
defined by Equation (5.1b) 

error in process gain determinations 
for steam flow rate disturbances, de- 
fined by Equation (5.1c) 

damping coefficient 

CONS. Canits! 5.01 40... 

density, gm/ml 

defined by Equation (5.1c) 


time constants defined in Equations 
(2h la) at, elic’6p) 


LaMmemconstant OT rirstoonder lag 
transfer function, minutes 


time delay, minutes 
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reset time, minutes 


time constant of a first order lead 
transfer function, minute 


phase angle, degrees 
frequency, radians/minute 


maximum frequency, radians/minute 
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APPENDIX A 
PHYSICAL PROPERTIES, EQUIPMENT DESCRIPTION AND 
CALIBRATION PROCEDURES 


Aas introduction 


This section contains a compilation of data required 
throughout the project, including the pertinent physical 
property data for water, methanol and water methanol solu- 
tions, the methods used to calibrate the column instrumen- 
tation and select the controller constants, a tabulation 
of the results comparing the present column operation 
with that obtained by Svrcek (128), examples of typical 
material balances obtained and listings of the loop re- 


cord formats. 
A.2. Physical Property Data 


The pertinent physical property data for methanol, 
water and methanol-water mixtures were subjected to a 
regression analysis in order to obtain a form which could 
easily be used in the digital computer. 

a) density of liquid water (96) 

DENSITY <ehi00G- 2 dawl0: T= 92.5110) 1. (A.2) 
Range: 40y<_Taesph2ed 
Dnatsi: DENSITY - gm/ml 
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fonsdtom vot Bieb ysreqorq Isoteyda dnsntti9q SAT. 
6 ot bstostdue svew esxutxtm ystew-lonsdgem bas votsw 
bfuoo dofdw mot 6 nkstdo of Y9bi0 nt eteylens nofeesyp9y 
.1atuqmoos Tstiptb odd nt beew od vt bana, 
(a@) yvotew biuptl to yttensb (5 — | 
(S.A) “T° "ofsta.S - T° "OfxAl.S-200.1 = ¥TI2N30 
OSf > T> 0)  sspnait 
ra\me - YTI2WAG = sed tnd 

a 


b) 


fe) 


density of liquid methanol (35) 
DENSITY = 0.808 - Balaar0e'] = ka98POesT* © 0A. 2) 


Range: Oto ELS anc 
Una ts: DENSITY - gm/ml 
1% OG 
density of methanol-water mixture (22,96) 
DBNSUTYCs  ha0c0s-e5ukexu0 «Lh tenleSiehO ° 1° 
(ASS) 
Selene | Oe Xe eo eld Or OX? 
Range: 40.< 1.55 b20 
O~< X°< 400 
Units: DENSITY - gm/ml 
i - OG 
X - wt % methanol 


saturated density of methanol-water solutions (128) 


DENSITY = 0.960 - 1.36x107°X - 7.46x107-°X? (A.4) 


Range: 0 < Ki 100 
Un Ges DENSITY - gm/ml 

X - wt % methanol 
density of steam (96) 
DENSITY. seeQe1]l3,431.862)0° °T - 1.07«x10 °T? 

(A.5) 
SU Nacy Oba MON a 

Range: LO ys be C0 
Units: DENSITY - lb/ft? 

T ~ OF 
heat of vapourization of water (96) 
HEATV = 1080 - 3.99x107!T - 5.62x107°T? (A. 6) 


Range: . 190°0< 7 -< 73500 


Units: HEATV - BTU/1b 
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g) 


h) 


i) 


j) 


k) 


heat of vapourization of methanol (35) 


HEATV = 291.2 - 3.31x10°!T - 1.53x1073T? 


Range: Oy cle cee OU 


Units: HEATV 
ip 


cal/gm 
68 


heat capacity of steam (43) 
HEATC =°7.80 + 3.2x107°T - 4.83x107°T? 
Ranges | 07 <" ) < 35500 


Units: HEATC - cal/gm mole o¢ 
T - C 


heat capacity of liquid methanol (35) 
HEATO =) U200) te luoe Ue Lets O.0Lx 10) o 1 
Range: 0.¢ [.<.120 
Units: HEATC -ecak/gm °C 

T OC 


heat capacity of liquid water (43) 
AeA Gos) 00 Wee leeS G0 we ey, OS le) bala 
Ranges .235 «< [s< 100 
Units: HEATC - cal/gm °C 
iy eit 


heat capacity of methanol vapour (43) 
HEALC (= [Ome clsgSx Leet oaus)..44107 a1" 
143.92%1L00 217 


Range: . 0 yoyT ax 500 


Units: HEATC - cal/gm mole °C 
Tae ac 
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saturated liquid temperature of methanol-water 


liqu 
mixtures (128) 
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TEMP us Sbhe29nen hedS eX gal 985241054%? 


(A.12) 
189;33 2 YOR? 
Range: 0 sXe an 0G 
212 e< 4 EMR <e14s.5 


Units: TEMPie- oF 
X - wt % methanol 


saturated vapour temperature of methanol-water 
mixtures (128) 


TEMP = 211.84 - 2.253100 fume 1279x100" °T? 
(A.13) 
= "92 S300 - > T° 


<ul 0G 
EMP < 148.5 


[AIA 
od >< 


Units: TEMP - °F 
X - wt % methanol 


enthalpy of saturated methanol-water vapour (128) 
ENTHV = 1143 - 5.29 X (A.14) 
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enthalpy of saturated methanol-water liquid (128) 
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X - wt % methanol 


enthalpy of superheated methanol-water vapour (128) 
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q) vapour pressure of liquid methanol (61) 


Loge se VPS 7h 'S7S635 COAG. 10 (Ac 17) 
(T + 230) 
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UA Gees VP - mm Hg 
TT) trajlic 


r) vapour pressure of liquid water (61) 


Log, VP = 7.96680 -_1668.2 (A.18) 
(T + 228) 
Range: 60c cai e00 
Units’: VP - mm Hg 
TonemeGic 


Regression equations could not be obtained which 
adequately predicted the following physical properties 
(94): 

s) equilibrium vapour composition of methanol- 
water system (128) as a function of equili- 
brium liquid composition 

t) equilibrium liquid composition of methanol- 
water system (128) as a function of the 
equilibrium vapour composition 

u) composition of liquid methanol-water mixtures 
(22,96) as a function of both density and 
temperature of mixture 

v) equilibrium liquid composition of methanol- 


water system (25) as a function of temperature 
and pressure of the system. 


A.3 Calibration and Instrument Adjustment Techniques 


This section will outline the procedures that were 


employed to calibrate and/or adjust the control settings 
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on the following instruments: 
i) overhead product composition analyzer 
ii) all flow transducers and their controllers 
iii) overhead pressure transmitter 


iv) both level transducers and their controllers. 
A.3.1 Overhead Product Composition 


The overhead product composition analyzer was cali- 
brated using the pycnometer method. The weight and 
volume of a set of pycnometers were determined from 
replicate measurements. The volume was determined by 
weighing the pycnometers full of water at a known temp- 
erature. Volume corrections were determined for the 
level of the sample in the capillary tube. The composition 
of the methanol-water samples was determined by the den- 
Sity of the solution measured at room temperature. The 
density is calculated from the known volume of the pycno- 
meter, and the weight of the pycnometer filled with sample. 
The sample was initially allowed to come to equilibrium 
with the room temperature in a sealed flask. The filled 
pycnometer was also allowed a few minutes to equilibrate 
with the room temperature before weighing it. The compo- 
sition was then obtained from standard composition-density- 


temperature data available in the literature (22,35,96). 
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The calibration of the overhead product analyzer was 
performed on-line. The column was operated to produce the 
desired range of methanol concentrations. The calibration 
technique followed that outlined in the Foxboro manual 
available for the Dynalog capacitance analyzer (94). The 
analyzer was initially operated with the temperature com- 
pensation disconnected in order to select the desired span. 
After this was accomplished, the temperature compensation 
was reconnected and the column was operated to produce a 
mid-span composition measurement. The temperature of the 
condensate was decreased and the amount of temperature 
compensation was adjusted, returning the composition to 
its previous mid-span measurement. The temperature of 
the condensate was Changed by recycling the cooler reflux 
pump output. back to the reflux accumulaton. NINis» proce- 
dure was repeated until a temperature change in the reflux 
accumulator did not cause any deviations in the measure- 


ment of the composition. 


After the analyzer span and temperature compensation 
have been set up as desired, the column was operated to 
produce compositions covering the complete recorder span. 
Samples of these overhead products were taken at various 
Chart readings and the compositions determined using the 
pycnometers. The results are shown in Figure A.1, complete 
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data. The parameters of the least squares line are given 


in Table A.1. 


The calibration procedure for the bottoms product 
analyzer, an on-line gas chromatograph, has been outlined 


by Berry (9). 
Aeo.e 8 86F lOW-Rates 


The flows are measured using orifice plates with 
flange taps. The pressure drop is transduced to a 3-15 
pSig pneumatic signal using standard Foxboro d/p cells. 

A number of flows, feed, reflux, overhead and bottoms pro- 
duct, use quadrant-edge orifices (130) since the discharge 
coefficients are reported to be constant at the lower 
values of Reynolds number associated with the smaller flow 
rates. Sharp-edge orifices are used for the remaining 
flows, steam and cooling water, since the discharge coeffi- 
cient is constant for the larger values of Reynolds numbers 
associated with these flow rates. The flow calibration 
data are expected to follow the characteristic orifice 
equation (130) 
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where FLOW is the flow rate measured in 1lb/min, AP is the 


K (AP) (A.19) 
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Sg is the specific gravity of the fluid flowing and K 


is the proportionality constant. 


All liquid flows were initially calibrated using 
water. The weight of water at a measured temperature 
passing through the orifice creating a pressure drop in- 
dicated by the chart reading was measured over a known 
time interval. The calculated flow rates, specific gravi- 
ties and chart positions were fit to Equation (A.19) using 
a least squares procedure. For all flows, the intercept 
calculated was not significantly different from zero, so 
the least squares fit was repeated forcing the straight 


line through the origin. 


The flow rates of the methanol-water solutions, namely 
feed, reflux, overhead and bottom product, can now be de- 
termined from Equation (A.19) using the specific gravity of 
these solutions determined from Equation (A.3) and the 


measured composition and temperature of the stream. 


The steam flow orifice was calibrated by introducing 
steam to the reboiler at various chart measurements and 
weighing the condensate collected in a previously weighted 
ice bath in a known time interval. These flow rates were 
fit to Equation (A.19), modified such that the specific 


gravity term was combined in the slope, using the least 
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Squares procedure. The intercept obtained for this flow 
equation was significantly different than zero and was 


thus retained. 


The calibration data and best fit flow equations for 
the different orifices are shown in Figures A.2 to A.7. 
The constants of these equations are displayed in Table 


A.1. 


A.3.3 Column Pressure 


The column pressure transmitter was initially cali- 
brated off-line to determine that the scale readings were 
correct. It was then installed above the vent line so 
that any methanol which condensed in the measurement line 


would return to the process. 


A.3.4 Condenser and Reboiler Levels 


Since a knowledge of the absolute values of the levels 
was not required during this project, the span of these d/p 
cells were adjusted to give a reasonable operating range 
for the levels in the column. Any levels were then mea- 
sured relative to the 3-15 psig output produced by the full 
range d/p cells. The parameters related to these measure- 


ment transducers are summarized in Table A.l. 
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TABLE A. 1 


Process 
Variable 


feed flow 
reflux flow 
steam flow 


overhead 
product 
flow 


bottoms 
product 
flow 


cooling 
water 
flow 


overhead 
product 
composition 


column 
pressure 


condenser 
level 


reboiler 
level 


Equation 
Types 


Orifice 


Si 


(Inch) 


0. 
0. 


SUMMARY OF MEASUREMENT DEVICE PARAMETERS 
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Span 
Orifice Inch 
Type H50 
quadrant 80 
quadrant 254 
sharp 156 
quadrant a 
quadrant 46 
Sharp 250 
= AY X + B 
= A X +B 


Equan 

A B Type 
noe 0.0 1 
42 0.0 ] 
29 On73 1 
naa 0.0 1 
ma 0.0 1 
.90 0.0 1 
O98 Doms Gm ue 
20 -10.0 2 


Y - measurement 


Units 
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lb/min 
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A.3.5 Column Temperatures 


All temperatures on the column were measured using 
iron-constantan thermocouples. Before installation, all 
thermocouples were checked against published thermocouple 
tables at various temperatures between 0 and 100°C. The 
temperatures measured by the thermocouples generally agreed 


to within 1°C of the calibration temperature. 


A.3.6 Analog Controller Constants 


Table A.2 summarizes the controller constants used 
for the various analog controllers. The settings for the 
feed, reflux and steam controllers were chosen to give the 
most rapid response with minimum overshoot for a set point 


change. The level controller settings were chosen to give 


TABLE (A.22 TYPICAL ANALOG CONTROLLER CONSTANTS 


Controller Ke Ty 
feed flow 200 Om 
reflux flow 100 teayans 
steam flow 175 Orlin 
overhead composition -200 ip 
condenser level 50 33 
reboiler level 300 Oe 
column pressure 100 ag 
feed temperature -150 Wes: 
reflux temperature - 20 fag 
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averaging level control, adjusting the product flows 
smoothly while allowing the level to vary between limits. 
The settings of the overhead pressure controller were 
Chosen to maintain as tight a control as possible, without 
Causing oscillations to occur in either the overhead 
pressure or cooling water flow. The reflux and feed pre- 
heat temperature controllers were tuned to give the most 
rapid possible response with minimum overshoot to the 


flow rate disturbances. 


A.4 Comparison of Column Operation with that of 
Svrcek (128) 


Typical results are presented in Table A.3 comparing 
the operation of the column attained during this study with 
that reported by Svrcek (128) before the column was dis- 
assembled and moved. Sufficiently good agreement was ob- 
tained between the runs to conclude that the present 
Operation of the column is consistent. The small differences 
in the temperature and composition profiles may be attri- 
buted to the inaccuracies and insensitivity in manually 
adjusting the pneumatic set points of the pneumatic control- 
lers controlling the input flow rates and the difficulty in 


reproducing the feed composition. 
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TABLE A.3 COMPARISON OF THE COLUMN OPERATING DATA 
WITH THOSE OBTAINED BY SVRCEK (128) 
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A.5 Material Balance Errors of Closure 


Table A.4 gives, as an example, a typical report 
issued by the program BALNC documenting the operation 
of the distillation column. A statistical analysis of 
the errors of closure obtained for each run gives the 


following limits on the average closure obtained: 


Egy (ge: ie al CAL, 
EMeQH SE] F2 + PZ 59 


€ Sr neon eae 
HA0 
This indicates that for 95% of the runs, the over- 


all balance closed within + 3.5% and the component ba- 


lances closed within + 5%. 
A.6 Loop Record Listings 


Tables A.5 to A.9 contain the listings of all DDC 


loop records utilized during the project. 
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TABLE A.4 TYPICAL MATERIAL BALANCE REPORT 
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APPENDIX B 


SCHEMATIC DIAGRAM OF DISTILLATION COLUMN 


This section contains the detailed schematic diagrams 


of the various sections of the distillation column system. 


ie 


suivaae ottemedce ieee RS We 
smet2y2 amulos notaei(idatb aa 


ese 


NOILIAS GVAHYAAD 
~NWN109 NOILWITILSIO 430 WYYSVIG JILVWAHIS L°@. 4aandTg 


| 
Osx Vanes 


6ui6604 


0 
8 
a 
6 
8 
8 
b 


(4eSuapuod yaa 
Quo [ana| aw Oud 
INIA jO13U09 yulod 
s 3eS 10 Hu16H604 


a : (22 
NOY 


: 
: 
aaa 


uleig 


Sik liGcsaadiabiccma | 


I 
CEs ia rath | YaLVM 
46 hi naan te memotiece 
; ? ral) D 6Hul6h 6074 
fie te | 
Qi) 


Cat ' {)-— SHA 
O-o-e 


sMMUJOD MOLTAJJIT2IG 70 MAABAIG DITAMZHI2 
AGITIIZ2 GAINRAVO 


lee 10 gaiggot 
iotines iniey 


$3H4 


See oD oe ey 


8 ED OD 8 ES SF Ob 8 Oe OP OT) 


f.§ 390917 


233 


NOTLOAS GIW 
-NWN109 NOILWITILSICO 40 WYYOVIG JILVWSHOS ¢°@ Junot 


(@u0 jaaa] SHUR} 0}) (2 ulesg 
LONaOUd dOL AN Se aL) Hk 


]04}U09 


iv 
yulod jas 
40 6u166074 ay) 
€9u4 
: j043U09 juUIod @ (me 
(1\ 3820 =@) Uy 
a 


D 
OES) 


tami ine 


2 on er ae er oe ee 


nad 


: “1 
fe2 to pnigg oi < fe 


x 


lovtnes tnies 


¢ © 345 
eon 


19 gaipgel 
inicg ise 
iovines 4 

§ 2A? 


vw 
iW 
ire  *sr 


234 


NOTLIAS WOLLOS 
<NWN109 NOILWITILSIG 40 WVY9VIG IJILVWAHOS €°@ 3ynNgI4 


Ch [ony uoljejnzje9 
i] syur} ee AH Bu660} 
L9NdGOUd 
WOL1LO8 
SS Ses |, See Ae 
oa (1» 
he 
uleig ey TY Bui6bo; 
} ua = 
q SNITOOS a. 
pee ©) 


a}esuapuod oom oom cam e HINA 
: iy = 


\ 

| 
&> q j eee gs eee | LT IP PLP ST SE 

WV3ILS Weay}s c 


— Bisdoz VW 


t ‘, 
e61nd i ‘ BAI], “}24 
Ries 
Daa hitb 
40 BulHbo; Cry, 


~~ 
o 
eB 
@ 
o)) 
be 
Ss 
a 
ag e Le | 
WE 


Ge 
“6 (2 TS = ae eee © gg 

@) aha 
®... 


nolisiusizs 


, MMUIOD MOPTASJITZIG 4O MARDAIG JITAMZHI2 €.9 3RUITI 
S WOITIZ2e MOTTOS 


239 


— 
U 
— 
5 U 
Coles 
a (a) 
> O 
O a ae 
ree ~- Aa z 
& Oo. ua o ws 
us co ta — > 
i VAY V 
s 
Cooling 5 
Water = 
Cooling S 
Water ro) 
HX} |} Hx] Y 
3 4 


FEED 
TANK 


a TANK 4 
BOTTOM TOP eN 
PUMP : 


fo 


FIGURE B.4 SCHEMATIC DIAGRAM OF DISTILLATION COLUMN: 
STORAGE TANKS 


Jat? 


b.8 aAUl 


_ 


nw 
co 
= ©& 
arr mo 
oxz> =< 
-Os <3 
cr + at 
=a 3 
ad < 
o= 
— 
a 
2 
nx 
<= 
= 
= 
Tr 
a 
— 
PF : / a - 
: a." A . ee ws ' ‘ y 2 ™m 
wteW grilood . Pas a ar ‘ ; P Ne 
twav[S>———— - aS \ 2 See ee ee - 
7 \, \ 2 : * ." ‘ace J 4 
, ; : y Vv  " \ erhey = pa 
; ' - } on Cn: a 
1 4 ‘ et 
a e f - . a 
— \ Eo ea \ - “A \ - 
\ 252 3 - ~ 
\ ‘ =} >! . 2 ~~ Co | ‘2 
\ ; . \ = v9 = 
. 4 ot wr sie, —__ 
. ‘ > : 
‘ q ‘ a =—J 4 my iy = 
™ : 
\ f - ) eS cs ar f oO 
\ Yo 4 a9 _ 4 SJ = 
\ 7 = 4 4\ f >) 
‘ i a’ \ j ~ 
. ’ em te ae 
\ \ _ ~ 


ces 


eos 


QN3941 
>NWNT09 NOILVITILSIO 40 WWYOVIG JILVWWIHIS G°@ JYNDTS 
YILNdWOD OL SINITNOISSIWSNWUL VIVQ s= => =m moe 
WVIYLS DISd OL mermmeramemm = SINI1SSID0Ud YILVM ONII00) oem oame 
INA seccceencance SINI7 LNIWNYISNI SS390¥d —————— 


SINI1T SS3I0Ud YILVM/IONVHL3Q ee SINIT LNIWNYISNI OUYWN 


SINI1 $S3908d 


YILNdWOD AG (U-UU JIGVIYVA) TOYINOD ONY 9NID901 dWNd U-d 


YING WOD AG (U-UU JFIGVIYVA) 9NID9OI YIINASNVYL 13A37 


JAIWA QI0N310S YIVIOVLNOI/YOLVIIGN! 13A31 


CD €) 


YIVIOMINOI/YIAYOIIN JUNLVYId WIL JAIWA QNVWH 


> 
Pe 


YICYOIIN JUNIVYIGWIL INIOdILINW (1192 G/P) YJINGSNWYL MOTS 


E) 9 @ & @ 


etn YIINGSNVYL sunssaud(qy (ty) Y3770Y1NOI/YIGYOIIY MOTS 


YIINOSNVYL JYNSSIYd ¥3090939N M014 


2 


HOLIMS JYNSS3IYd = Sd (aqosd arueyiseded) YIINGSNVYL NOILISOdWOD 
JAIWA INILVINDIY IYNSSIYd AYd (@2Uej19EdlD) YIIIOYINOI/YIAYOIIY NOILIS Od WO) 


YIVIOMLNOD/AYOLVIIGN!I 3YNSS3IYd CIP (udesboyewosyd) Y30Y0II9Y NOILISOd WO) 


PDODGH © 


- a = _ = = 


— twas 32022288 pe} Fe 


4 


A soonenta snue 230e('*) be | 


oat 7 = aN 
laa JAUTARIT MIT THIOUTIUM - iie2 Gib) RIDUGZMART WOIT (*) ; 


m 
@ 


a 
| AUI09THODIAIOROIIN 2 }9UTARIOWS ) WVIA¥ GOVAN VHC. : 


3VJAV 01043102 (v2) BIJJOATHOIHOTADIOM! J3V3) Gd) 


Tre =: A3TUIMOD-¥B-ta-nn WAAIRAY) OHIQI00 FS RIQUGZHART 33V9) 


$IWGM0I YR (n-Mn ZIGAIRAV) JOTTMOD GHA DMIAIOI GFy qu 


7 ¢3# 11 2639085 
ei¥#! J 23 ORF AIA AW\ IGHAHTIN CorempyT etter Ee 2713 THIMURTZA JETAMU SUS ee 
- -_ a i 
oF ee : THY eecovsesasaes 2241 J THIMURTEH 2230038 
7 MAZRI2 Dl2S OF sence emcee 23414 22IDORI RITAW QHIICOD wees 
AITUIMOD OT ZIVIS MOILSZIMZHART ATAQ eee we 
sMMUJOD AOITAJJIT2IG 30 MAADAIO JITAMZH32— = ¢ {AUST 
é J gAl ? : Joi 4 


des 


237 


APPENDIX C 
PROCESS GAIN DATA 


This section contains a summary of the experimental 
data used to determine the process gain portion of the 


open loop column model. 
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APPENDIX D 
OPEN LOOP DYNAMIC MODEL DATA 


D.1 Pulse Test Results 


This section contains a summary of the experimental 
data used to obtain the dynamic parameters of a first 
order plus time delay model used to describe the open 
loop response of the distillation column. The frequency 
response diagrams displayed in Figures D.2, D.3, D.5, 
D.6, D.8, D.9, D.11, D.12 and D.13 were prepared from 
the measured transient response of the column to rectan- 
gular pulses in the various input variables using the 
Pulse Testing Analysis Program (PTAP) (135). The dynamic 
parameters of a first order plus time delay model were 
determined using both the ‘Levy' procedure and the proper- 
ties exhibited by first order and time delay transfer 
functions. These parameters were compared to those cal- 
culated from the transient response data using the Rosen- 


brock procedure (110). 


The parameters of higher order polynomial transfer 
function models have been calculated using the 'Levy' 
procedure. The parameters are presented in Tables D.11, 
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FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 
COMPOSITION TO POSITIVE PULSES IN THE 
FEED FLOW RATE 


TABLE D.2 
Feed + 

GAIN 0299 
FREQ MAG PHASE 
RAD/ RATIO ANGLE 
MIN DEGREE 
020020 14000 —-367 
020023 0999 4.3 
020027 04998 561 
00-0031 00997 -5 09 
0°0036 00996 -609 
000042 00994 -8.0 
0e004S 00992 —-9 63 
Oe0057 Ceo9EF —-10e8 
020067 046985 =-12e6 
000078 046979 =14e6 
O0e0091 0606972 -17600 
0520106 02962 —-1906 
000123 00949 =-226¢7 
000144 02933 -26e2 
020167 Oe912 —30el 
00e0195 06885 =-34e5 
Oe0227 02853 —- 3903 
060264 06814 =44.2.6 
020308 04769 =-5003 
000358 Oe719 =-5604 
020417 06664 =-6268 
0°0486 Oe606 ~69e4 
02-0566 04548 —~7602 
020659 00490 -8343 
OoO0767 Co435 —-9006 
020893 06383 -98el 
0el04C 046335 -106el 
0e1210 0e292 -114e6 
0+«1409 02254 -123.7 
021640 04220 -13306 
O0l1910 001990 =-144e6 
002223 00164 -15609 
002588 Del4l -17008 
Ce3012 Oel22 -18625 
003507 042105 -204e5 
0064083 042090 =-2256042 
004753 OeO077 =—248e9 
025534 02066 =-276¢4% 
006442 054057 =-308el 
007499 02049 =344e9 


Feed20 
2093 
MAG PHASE 
RATIO ANGLE 
DEGREE 
1.000 -203 
12000 -2el 
1.987 -36 
0e999 ~366 
00998 —~4e2 
00998 409 
Oe997 528 
00996 ~6e7 
00995 -728 
00993 Gre)! 
00991 ~-1006 
Oe, FBi. =w263 
00984 =—-1404 
0e978 -l60e7 
00970 -19 «4 
00960 =-2265 
00947 =-26020 
00929 -30.1 
Ue906 =34e7 
Oe&77 =3909 
02841 -4526 
0796 =-5128 
00745 —-58e3 
02690 -65el 
06536 -72el 
00585 =-7948 
06540 —“89el 
00484 -101.0 
02«404 -114e4 
006310 -124e4 
00253 -127¢5 
06232 -134e9 
00228 -14246 
0e278 -175-0 
Oel37 =—228e2 
16048 =-215e9 
02089 -157e6 
02086 —-268e2 
02030 -101.0 
00069 -54e9 
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204] 
MAG PHASE 
RATIO ANGLE 
DEGREE 
1-000 -320 
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020999 ~4el 
0«9S9 =Le7 
00998 ~5 05 
00998 =—604 
0.997 ~7e) 
Ge«996 ~8e7 
02995 = 10:02 
02993 -11-8 
O0e991 -1308 
02988 = li6ie,0 
0°984 -1867 
00978 -“2le7 
Oeo970 ~250e2 
06960 -29.03 
02946 -34e1 
00928 =-3906 
006905 ~45e9 
06873 ~—-53el 
06833 ~61¢4 
Oj, 783? «= = 7-Ore, 7 
Oe 215  =Bilye,1 
02648 9265 
00567 —-104e8 
02483 -117.7 
C e403 -131.9 
O0e327 -148e08 
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06109 -—41lel 
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02096 T1409 
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00057 -13300 
OeV075 -166062 
02025 -18104 
02042 -198e3 
02044 =-10109 
02059 88e9 
00490 -23401 
020060 =-94e5 


bSboe4 
$8-é 

32An4 OAM 

BJOWVA OLTAR 
33RO39 
Teh* COOes 
LeE~ 22 CeO 
Let PeL«D 
Seon OY ies 
Gaep~ ee2eeO 
Bee= BPE +O 
,ed~ BRL eO 
Gei= Fee ed 
[a@= | 8CR2d 
deOl[—= AlleO 
SeSie Sead 
Pehl O8C.O 
Bedli-= 682«0 
Cvti~ [8ee0 
,eSS~ evee«0 
eodS~ 8a¥e0 
Te0E&= He oO 

eCE= BEF «0 
Calpe Ti¢sO 
x2 be 86620 
8.ce~ [chs0 
@sad~ E0829 
He hl ~ £7) «0 
Ha Gb= 3020 
SeTe= OF S00 
LeG@Ql=— C1440 
2eOlI=- ESE«O 
Te@Sl- EVSeV 
€eB8SI~ SSS-0 
&«TEL~ 62i+0 
VeGei= Acie 
sav l= alOed 
QsEEL@= TeQs0 
Seddl@= 2@5040 
PelBle 040 
EsHel= SHOe0 
CelOl= a4020 
“?@e86 P2020 
LotES= GPAeD 
2e a= 02000 


esbaat 


taes . 
32AH9 - BAM 
34DKvA OLTAR 

338030 

Oet= 000«! 
@est~ eeCeO 
Lad= CeRaO 
vaa~ Cee O 
eat BeLseG 
feo" BEL2O 
Ge¥= TO. G 
TeS= 62ee9 
SeOl-= e@ead 
BelLi~ ERPR?«D 
Ss€iq~- L@@«0 
Osdl= 88040 
,«Bhi-=~ AB2 sO 
ee Oe BT CeO 
Sets= OT?«Q 
E«CS— Cadel 
[ere aaSeO 
QeQVE~ BFO.0 
CeG@a= O40 
LeE@=- E640 
e#ef[d= EEBs 
YeOV~ Ears 
L.efla- LST e2O 
2.oe~ Sede D 
5BeeOl~ Vdeed 
T.tLli- EGPed 
Ve lElm EGeed 
BeBal= TSE eV 
@eOTi~ EtSeD 
HeAe@l~- BELeO 
SeeT l= €60e0 
E«lTBl~= 38le0 
fela~ 2OLs0 
Cel Teied 
Cent 8@0eS 
@sSli~ LSféel 
E»Sil= 82000 
SeV¥LS= £000 
Eeol EaO«d 
Cele PTO00 


et 


32AHS 
310nA 
339930 


Ee S* 
TeS- 
aeb= 
62t=- 
Sod 
Cra= 
&.é- 
Ted~ 
8. T= 
i.@« 
ae0i- 
GeS{~ 
eaadl~H 
F2dl~ 
oe elm 
@e$S~ 
Oedh~ . 
£«GE= 
Tuabe 
Ga vt= 
deda~ 
66{t-= 
€.ae~ 
I.¢o~ 
[eST= 
BaelT~ 
[+@S~ 
Oe fOLl= 
Pehil~ 
PeeSl= 
éeTSIi- 
Ca HE l= 
oeSei~- 
O.¢Tle 
&-8SS~- 
@.2@fS~= 
aeTel~ 


$e8aS- 


OeidQ{= 
Creat 


Doe % 

Eas 0 i" 

; Lea. oay ates 

OLTAS Bi0MA « SITAR 
33a i 


O0Ose! 
C000! 


~T8Cel 


SVL ed 
BRC «OD 
BSC eG 
Tel«d 
dled 
ee?e0 


ELL 6D 


lee@e0 
BEE«O 
aBleS 


~—6«dBT SD 


OTs20 
0a@eO 
TAeCcO 
eS eed 
aye 
TV Hed 
[#800 
ATO 
feT sd 
Olse0 
GE de0 
&8ee0 
OFeeD 
®6aa0 
POReD 
OLE.0 
EES«D 
SESeO 
8SS20 
&VSe«0 
VELlsG 
GP0el 
28000 


88040 


OEUs0 


2020 


Teeo 
Esa 
Leé- 
VuGe> 
Ream 
0.8~ 
Eet= 
8.0l= 
deSi~= 
deal~ 
Cotl= 
d2el= 
TeSS~ 
52dS~ 
LeQE= 
GepE= 
Esee= 
ose e— 
Ese= 
aeac= 
8aSe~~ 
AaQo~ 
Seat 
EeE8= 
3e92= 
Le8e~ 
{.d0l= 
deal l= 
TeESI- 
@elEL<~ 
ded dfeq 
Gedel[= 
8.0Ti-~ 
éed61i- 
2ee00S= 
$sa@SS= 
Re Besa 
PedTS~ 
Le80E=— 
ReasE~ 


269 


WADLE® D..3 FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 
COMPOSITION TO POSITIVE PULSES IN THE 
REFLUX FLOW RATE 
Reflux + Reflux20 Reflux22 Ref lux23 

GAIN 0099 14280 13029 10°59 
FREQ MAG PHASE MAG PHASE MAG PHASE MAG PHASE 
RAD/ RATIO ANGLE RATIO ANGLE RATIO ANGLE RATIO ANGLE 
MIN DEGREE DEGREE DEGREE DEGREE 
020020 14000 -2e0 14000 -20e3 12000 -2e2 12000 = Velo 
060023 00999 204 00999 “207 00999 “205 0-999 -168 
000027 04999 ~208 06999 -3.1 0.999 —-3e0 0-999 =-Zel 
020031 02999 “302 00999 =357)00999 @—3%5902999 “204 
000036 00998 ~308 00998 “403 06998 “4e0 0-999 -209 
020042 02997 —“4e4 02998 =5e0 02998 —“4e7 00999 363 
020049 040996 = 517009 9H =5%8'00997 505 04998 -309 
0e0057 00995 “600 020996 -6e8 046995 -6e4 02998 -465 
Oe0067 020993 -7e0 00995 ~7e¢9 0995 =705 00997 “503 
060078 04991 -8el 0.993 —-Ge2 0-993 —-8e7 06997 —“6e2 
Oo0091 0-988 “904 006990 -100e7 00991 “1001 02996 all OW 
000106 02983 ~-1lleO 0098G =-1204 OF987 ~“lleB 00994 -8 0% 
000123 De97T7T H120e7 00982 =-14e5 00983 =-13067 00992 -9el 
020144 0970 -14.7 Oe975 -16«8 Oe977 1509 00989 -11e3 
Oe01l67 O0e960 -17e¢0 00966 “1966 04969 —-18e5 00986 -1362 
Ce0195 00946 “1967 04955 ~22e! 02958 -2125 040.981 —-1503 
020227 00929 -2206 00939 =26 04 00943 —-25e0 00974 -1728 
000264 Ce907 =2600 06918 =“3006 00924 —-2900 00965 =—2007 
000308 06880 =2906 0e890 =-3505 06898 =3305 00953 =24el 
C5o0358 CoB47 =3306 00853 -4120 00865 —-38e7 00937 -2709 
Ce0417 Oe807 ~3769 02805 -47.3 Oe821 4405 02916 -32e3 
000486 Oe761 =4204 00743 5402 06766 “5009 06888 =3704% 
000566 Oe710 4701 Oe666 -61l0e6& 00697 “5708 00852 —-4301 
000659 00655 =510e9 00573 =690el 08615 “6407 O0807 4905 
OeO767T 00597 =560e6 00468 =-756¢5 00523 =-7009 04750 ~560e5 
020893 006539 =Sle2 00452 =78e2 00431 =“75e0 00681 —b4el 
001040 00481 6506 00284 =-7306 00357 =-75e¢7 00601 —-7200 
001210 00427 6909 006265 6548 06313 “7409 00513 7907 
021409 06375 =7309 Oe277 “6800 00287 —-770e8 00422 —8607 
001640 006329 =7708 06254 —B8lel 006243 —-84e7 06335 9168 
061910 0+286 —-81.65 00167 —-92e7 0-178 —-84e¢9 02260 —-9306 
002223 06248 9 -SSe) Ool¥S © 65.63 OC L664 “7308 OF2027 ~“9lel 
002588 OsZIS= -8606 O'elh656 27200 Oel71L2 —86585 O619285 —9007 
003013 OF P86" —92e2 OsOPieVIOeF OelOlLS <SLlel OC L609 -97e9 
O0e3507 06160 9508 OolO7? H27e7 Oolll =77e6 Ooll7? 9109 
004083 00138 “9905 02061 —224e5 00064 ~“61e0 00131 —-95e8 
004753 O0llD -1030e5 00462 -17909 02062 “11265 02066 -11301 
005534 050102 “10708 00el69 —1400¢8 02304 -1490e2 046064 -1020e3 
006442 00088 -112+«6 02094 —-9368 0-082 —-1570«2 02-076 -~260e8 
007499 00075 11708 00097 —-8563 0eDG1 “12109 00142 5301 
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TABLE D.4 


FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 


COMPOSITION TO NEGATIVE PULSES IN THE 


STEAM FLOW RATE 


GAIN 


FREQ 
RAD/ 
MIN 


0260020 
020023 
000027 
020031 
020036 
020042 
020049 
020057 
020067 
0.0078 
0.0091 
0.0106 
00-0123 
020144 
020167 
02-0195 
000227 
0e0264 
0290308 
0e0358 
0.0417 
060486 
0«90566 
0.0659 
00767 
0.0893 
0261040 
021210 
021409 
021640 
021910 
002223 
0.2588 
003013 
023507 
024083 
004753 
025534 
006442 
067499 


Steam + 
0099 
MAG PHASE 
RATIO ANGLE 
DEGREE 
1-000 -206 
0.999 = 3 el] 
00999 =-306 
02998 “4.2 
06997 4.9 
02996 = 5el 
00995 -606 
00992 =e 
02990 =i)'9 0 
Ce987 -1064 
006982 -l2el 
00976 -14el 
06968 =-1603 
OeS57 =~-1809 
02943 -21.8 
00925 =25el 
00903 —-28e7 
0e874 -32.8 
06840 -37643 
00799 4201 
02752 “4762 
Oe700 —5265 
00644 -58.20 
OvetoS 65 . 6:3 2G 
O0e527 ~E9el 
0e470 =Thel 
00416 —80602 
00366 —85-8 
00320 =91¢5 
Oe278 ~97¢3 
00242 -10343 
00209 -109e7 
Oe-l18al -11646 
00156 =-124e1 
00134 -122¢5 
00115 -141-8 
0009S -152¢3 
00085 -164¢3 
00073 -178e0 
Ois:016 35 193 on! 


Steam20 
19270 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 -206 
00999 -3.0 
00999 =305 
02.999 -4.1 
00998 -4.8 
02998 -5 26 
00997 -6.5 
02996 -7e5 
00994 ~8.8 
02992 -10.2 
00990 -11.9 
02.986, $13.9 
0-981 -1602 
Oo974 =-1868 
00966 —-2168 
00954 ~2504 
02938 -—29e5 
00916 “3402 
056888 3949 
00852 —-4526 
026805 =-52.5 
006745 60041 
Oe672 6802 
00587 “7645 
0e494 —-84e2 
0 e403 ~89 69 
00330 +930 
O0e282 9565 
00241 -100.7 
O0e189 -105-20 
01:5 95 4100 2 
00165 -10644% 
Oel21 -123.0 
Oe108 =-108.5 
00098 =-14024 
Die OF? 99 1272 onl 
00036 -123-0 
02039 7803 
02096 -173642 
0e051 -188e4% 


Steam22 
28074 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 -1e8 
0999 -2el 
00999 -205 
020999 -209 
0.999 —-304 
00999 “400 
0998 -406 
00998 m5 04 
02997 —-6e3 
02996 -7e.3 
00995 —-8e9 
00993 ALOE) 
Oo991L1 -1106 
O0e987 1365 
Oe 983 ~—15e7 
Oe977 —1Be2 
0«969 ~21le2 
00959 =24e6 
02945 -28e)5 
00926 =3340 
0902 —-38el 
Oe871 -4308 
e881. —S0e«2 
Osa 2 ue) 6:2 
O«/724 ~64e7 
00659 =-7203 
020593 1909 
0% 22 1 “~87e2 
O0e«481 =-95e1 
00437 -104e7 
0e«388 ~116e6 
00332 ~129¢2 
00283 -141.8 
00244 ~15609 
0e215 -~176e7 
Ondo ~—21662 
00074 -192¢8 
02848 408 
02044 —145¢3 
00053 =139e3 
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FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 
COMPOSITION TO NEGATIVE PULSES IN THE 


STEAM FLOW RATE 


GAIN 


FREQ 
RAD/ 
MIN 


020020 
0200023 
020027 
00031 
020036 
020042 
000049 
020057 
00067 
060078 
02-0091 
0e0106 
0e0123 
000144 
000167 
000195 
000227 
020264 
00e0308 
020358 
020417 
02«0486 
020566 
000659 
020767 
020893 
021040 
0e1210 
001409 
021640 
021910 
002223 
02e2588 
003013 
Oe 350-7 
024083 
004753 
025534 
006442 
007499 


Steam24 
21.11 
MAG PHASE 
RATIO ANGLE 
DEGREE 
1-000 -220 
0e999 -203 
02999 ~2e7 
0.6999 =el 
00999 -3e7 
0.998 ~463 
0998 500 
0e«997 528 
00996 —-608 
02995 -7e9 
0993 —-9e2 
00991 -10e7 
06988 “1265 
06984 —-1425 
0.978 -16e8 
O97 le = b0 eid 
02961 -22e7 
02948 —-2603 
00931 =3004 
00«908 35640 
06880 -4063 
006844 4661 
02801 -5265 
Oo751 =5904% 
00695 =66e7 
00637 7405 
0.577 -8209 
0.051139 920d 
00444 =-101le/7 
Oier397 Te i Odie 6 
00339 “11704 
Oe317 $1324 
0249 —-15563 
02149 -16845 
O~d 2 le -1 5509 
00145 -19169 
00074 -1290e4% 
0«l86 -21209 
00377 6920 
002025 -17508 


Steam2/ 
17218 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 “1.9 
00999 -2e3 
00999 —-26 
Oe999 -3el 
00999 -306 
02999 ~4e2 
02998 469 
0.998 “567 
00997 -626 
02996 -T7el 
02995 -9.0 
02993 —-10«5 
O99 3s caleeie2 
0e4¥88 =-—l4e2 
00984 -1626 
0.979 ~19.3 
Oe971 -2204 
0962 —-26el 
00948 =-3003 
0.0931) 352 
00907 =—40¢8 
02876 -47.3 
00836 -54e7 
00784 6361 
Oe718 =7204% 
0«6328 8246 
00544 —-9302 
00441 -10363 
02341 -111.0 
06263 -11405 
06218 1158 
0.191 -121.5 
O00«155 =-130¢9 
O16 79 =13 bicie 
0.00969 -138.3 
02085 -143e1 
000735 -159%2 
00559-15849 
00056 -179e4% 
02064 -191¢5 
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TABLE D.5 


FREQUENCY RESPONSE OF THE BOTTOM PRODUCT 


COMPOSITION TO POSITIVE PULSES IN THE 


FEED FLOW RATE 


GAIN 


FREQ 
RAD/ 
MIN 


02e0020 
020023 
000027 
020031 
020036 
020042 
0e0049 
020057 
020067 
020078 
020091 
000106 
00-0123 
000144 
020167 
000195 
000227 
020264 
02«0308 
020358 
0e0417 
020486 
00566 
020659 
OcO0767 
020893 
0+1040 
021210 
021409 
021640 
021910 
002223 
022588 
023013 
0«3507 
024083 
004753 
025534 
006442 
027499 


Feed +B 
0.99 

MAG PHASE 
RATIO ANGLE 
DEGREE 

1-000 -1.8 
0e999 =Zee 
02999 -2e5 
00999 -209 
00999 ~304 
° 02998 —-4.0 
00998 ~4el 
00997 —5 04 
00996 seal 8 L”, 
02995 —-704 

00993 “Se 

Oo991 -100¢0 
0988 -11.6 
02984 -135 
OCOFE* =1S67 
Oe971 -18e2 
02962 =2i¢O 
02949 -24,3 
0933 -28el 
OTS P2 3205 
0«886 ~-37el 
02854 -42 64 
0-816 “48 e2 
Oe772 5405 
Os72Z1- 6143 
02667 ~6846 
00609 -76e2 
O°e-5/Sid) —-84.2 
0e493 -92e7 
06438 -101¢6 
00386 =-11160 
Cease = 2) Bek 
02e295 = a2 ki 
00256 -144.1 
Or Zi2i2s <= LSet 
Cte IiZe S— ie ok5 
00165 -18943 
00143 -208e6 
00123 -230¢5 
O«l06 =—255e7 


Feed20B 
4063 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12900 eo 
0.999 -2e7 
0.999 -31 
0.999 =3 6 
00998 -4.3 
00997 -5 20 
00996 5 8 
00995 Golo LTA 
0.993 mihi ys) 
0.991 -9el 
00988 =-1046 
06983 1263 
Ceo) S1aes 
06969 =-1646 
02959 ~-19 23 
02945 -2203 
0-926 —~25e7 
02902 =2945 
Oce871 —-3367 
00833 =38e2 
0078S 4249 
Oel737 —~47 64 
02685 oo a) 
Oe641 =-55e2 
00614 —-59-20 
02605 “64.8 
00596 =746¢5 
02552 -88.3 
00455 =-1026¢5 
Caso = 109s! 
00321 —114-8 
00242 -“136 64% 
O-l22 -99.7 
00245 -1U503 
00433 -&7.2 
3-963 -14204 
Oo116 “465 
O«l22 5309 
0-108 14.0 


06311 


-53.3 


Feed24B 
5el2 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 -l1./7 
0e999 —-1e9 
0.999 -2e3 
0999 =2e/7 
00999 —-3el 
00999 —306 
0998 ~4e2 
00998 —-409 
00997 mae 4 I 
0.997 —6e/7 
00995 =-728 
006994 -Jel 
00992 =10e5 
0e989 -12e3 
0986 -14.3 
0e981 —~16e6 
009714 -1903 
0.966 —2204 
02954 —-2620 
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0«919 —-34.8 
Co893 —-4Cel 
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Oehnrt ~59 48 
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00656 7544 
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06538 “92-9 
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Vesee "= hs'reS 
006264 -157e67 
Maat *= POS 6G 
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FREQUENCY RESPONSE OF THE BOTTOM PRODUCT 


COMPOSITION TO POSITIVE PULSES IN THE 


REFLUX FLOW RATE 


TABLE D.6 
Reflux +B 
GAIN 0.99 
FREQ MAG PHASE 
RAD/ RATIO ANGLE 
MIN DEGREE 
020020 1-000 -2el 
020023 0999 -205 
000027 00999 -209 
060031 0.999 =304 
020036 04.999 —-309 
020042 02999" —-406 
000049 0-998 =5 04 
00-0057 0-998 —-603 
020067 00997 -7.3 
0.0078 00996 —-8e5 
02-0091 O<6995 -9e9 
020106 0-993 -11.25 
060123 046991 -1304 
000144 046987 -1506 
00-0167 0-983 -18el 
000195 06977 =-2140 
O0e0227 009790 2404 
000264 04960 =-2863 
020308 06947 =3268 
00e0358 00¢930 =-38e0 
000417 02908 =4348 
020486 0e88l -50 «4 
020566 02848 -5708 
0e0659 Oe809 —66el 
OoO767 Oe764 =—7503 
0.0893 0-713 —-85 04 
0e1C040 02658 —9605 
0e1210 02600 -108e5 
021409 04542 -121le7 
021640 02484 -13602 
Oel1910 00429 =-152.0 
002223 046378 -169e56 
002588 04331 “18902 
003013 GOaZas ~2bb<3 
003507 042250 =-23602 
024083 COe2b “26406 
004753 040187 =2976e0 
025534 000162 -334e3 
006442 02139 =-3770e2 
007499 00120 -426¢8 


Ref lux20B 
5062 

MAG PHASE 

RATIO ANGLE 
DEGREE 
12900 -20e4% 
00999 -208 
00999 -3e2 
00999 -3.8 
0.999 =~4 4 
00998 -5el 
00998 -60 
0e997 -7.0 
00996 -8el 
0«S95 9.5 
00993 -11.0 
00990 -12.8 
0.987 -14.9 
06982 -17.3 
Co976 =2002 
02968 -23 64 
OreS5.7) Caeihee 
02942 = Bho 
00923 —3604 
00898 -42el 
00855 ~=—-45e4% 
Oe824 5563 
Oo775 ~6209 
Oe717 -—700¢7 
00«658 =—78e5 
O60 059:509 
Oe572 =9306 
00565 =-10420 
00567 “1200 
00536 -14360 
00436 -170el 
00293 -188e3 
Oe291 -192e61 
00335 -2356% 
00084 —26307 
Oe B33 2250's 
16435 -176e6 
Ces92.. #2 Py e2 
00124 -219e1 
00255 -119«6 


Reflux22B 
7013 

MAG PHASE 

RATIO ANGLE 
DEGREE 
1-000 -2e2 
00999 =2e6 
0.999 =30 
0.999 ~305 
02999 —4el 
0-998 —-408 
00998 m=5e& 
0-997 —-6e5 
02996 =-7e6 
06995 “809 
0.993 = 1030/3 
0.991 -12-0 
02988 —-140 
0984 ~1603 
00e97T9 =-1Be9 
Oe97L =—22¢0 
06962 —-2545 
02949 -2966 
06932 =34.3 
00909 “3967 
O«8EO —-45e8 
00843 =-5267 
0e796 -SUe3 
Oe741 —-68e5 
Ge678 -77e2 
Ceo61l2 “8601 
06551 -95602 
00499 -105e2 
00«454 —-1l1l7e7 
00402 ~133¢6 
02332 -150e9 
Oe273 -164e7 
00e«260 -182ed 
06229 2 thev 
00205 —226e¢e 
Og262, 263.00 
Oe377 jaa 
00543 -15104 
00195 =-240¢8 
0e118 413 


Ref lux23B 
633 

MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 -le7 
12000 -200 
06999 -203 
00999 Cael at 
0e999 ~Ze2 
0-999 -3e/ 
00998 -4.3 
0e998 =5 20 
0«997 508 
0.997 =6 «6 
00996 -7e9 
00994 “Yee 
046992 ~-10e7 
06990 125 
00«986 -14e5 
0-981 —-1609 
Onee77 5S: 1926 
O06 ~—22e7 
O.8996 =-2663 
06942 -—30¢5 
006924 —-35e01 
Oe9V01 —~40e4 
00875 ~4643 
026844 —5208 
06810 -60el 
Oe776 “S804 
00743 =-78e0 
Oei0e8 =—-9Ce0 
02661 10409 
0.95.36, ~12380 
004380 -14265 
006376 -159e2 
Oe ae& —176e6 
Oagls ~205el 
O0e219 —225e2 
Ore2G ~25863 
00239 4908 
02309 1504 
Oe595 ~-10el 


00235 


—-32el 
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TABLE D.7 


GAIN 


FREQ 
RAD/ 
MIN 


020020 
00023 
020027 
00031 
020036 
020042 
000049 
020057 
020067 
020078 
0°0091 
000106 
000123 
020144 
02-0167 
000195 
000227 
000264 
000308 
000358 
02.0417 
020486 
00566 
000659 
020767 
0.0893 
021040 
On l2l0 
021409 
0261640 
021910 
002223 
022588 
023013 
0.3507 
024083 
0204753 
025534 
006442 
007499 


FREQUENCY RESPONSE OF THE BOTTOM PRODUCT 


COMPOSITION TO NEGATIVE PULSES IN THE 
STEAM FLOW RATE 


Steam +B 
0099 

MAG PHASE 
RATIO ANGLE 

DEGREE 
12000 -1.8 
00999 =2e2 
00999 = 252 
00999 -20e9 
00999 -304 
00998 “420 
0«998 “4e/ 
06997 ~5 04 
00996 -60e3 
0-995 —-7e4 
00993 —8e6 
00991 “100 
026988 =-1lleé 
02984 -135 
0.978 -15e7 
0e971 = 1 8.02 
0962 -21.0 
00949 =-2403 
00933 -28el 
06912 =32e3 
06886 -37el 
02«854 -420e4 
Oe816 =-48e2 
Oe77T2 =5405 
Oc72l..76le3 
O0e667 =-68e6 
0e609 -7662 
0.551 8402 
00493 =-92e7 
00438 =-101¢6 
02386 “111-0 
06338 =-1l2lel 
002554 wl ere) 
00256 -144e0] 
Oe 22 24 e Ice 3 
Oe LO 2, mal 2:05 
Ove 165, 1:8 9:0 3 
0-143 -208e6 
O10 123, 2305 
0«l106 =-25567 


Steam20B 
19-09 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 =-1.9 
00999 “203 
00999 -226 
00999 =-3el 
Oe¥99 -3.6 
00999 “Ge2 
0-998 ~409 
Oe997 -5e/7 
0699 7 -6.6 
00996 -7./7 
00994 -9 0 
056992 =10¢5 
0098 5.5 “2262 
Co986 =-14e2 
De 815 =1 bie 
00974 -19el1 
0965 ~22e2 
020954 —-2508 
00938 ~29e8 
00918 =34e5 
056892 ~-39e7 
O0«860 =4546 
Oe Slo. a =pZ2¢1 
Oeo772 —59el 
02718 —~66e7 
02660 -74.7 
02601 —-8302 
00541 9207 
0047-7 =-10361 
00408 -1136¢5 
06350 -122¢6 
O0e318 -1339 
Oe 2,439 lo le 8 
Co213 -164-5 
Oe 20, 13005 
Ool79 =-20244% 
00194 —-240e60l 
00203 —-26846 
Ore L695, “620 0 
02066 -14766 


Steam22B 
13640 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 -204 
00999 eat RS 
0999 RE) 
00e999 —-328 
ee) “465 
00998 = ee 
00997 =Se 1 
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FREQUENCY RESPONSE OF THE BOTTOM PRODUCT 


COMPOSITION TO NEGATIVE PULSES IN THE 
STEAM FLOW RATE 


Steam24B Steam27B 
GAIN 19.95 14624 


FREQ MAG PHASE MAG PHASE 
RAD/ RATIO ANGLE RATIO ANGLE 


MIN DEGREE DEGREE 
0e0020 1-000 -128 14000 -1.7 
020023 00999 =2el 00«999 “109 
020027 04999 ~205 046999 =203 
020031 0999 226942998 -2e6 
020036 020999 304 00999 2301 
000042 00999 “460 0999 -3 6 
020049 02998 “466 02998 —-4e2 
020057 O0e997 =5e4 02998 409 
020067 0997 “603 06997 =Se7 
0200078 00996 -7.3 00997 -626 
020091 02994 —-8e5 02996 -7e7 
020106 00992 -9 09 00994 =9 0 


020123 04990 —-11e6 00992 ~10.5 
000144 04986 -1365 06990 ~-l2e2 
OceO0l67 00981 =-150e7 00986 =14e0e2 
000195 00975 -182e2' 04981 —=-16e5 
000227 04966 =2lel O8F75 “1902 
000264 00954 =-2405 00967 ~=2203 
000308 0939 —-280e4 02955 ~2508 
000358 046919 =-3209 02940 2909 
000417 02893 —-3709 O0e92l —~3465 
000486 00859 =434e6 00896 =39e8 
020566 050816 “4909 00865 —-45646 
020659 04765 =560e8 00826 =-52el 
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TABLE D.8 


FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 


COMPOSITION TO NEGATIVE PULSES IN THE 
FEED FLOW RATE 


GAIN 


FREQ 
RAD/ 
MIN 


020020 
0.0023 
020027 
0e0031 
0°0036 
02-0042 
020049 
000057 
02-0067 
020078 
00091 
020106 
020123 
00-0144 
0e0167 
00-0195 
000227 
020264 
000308 
000358 
0.0417 
00486 
0290566 
020659 
020767 
0.0893 
021040 
061210 
Ce1409 
021640 
021910 
0e2223 
0.2588 
0+«3013 
023507 
0.4083 
004753 
025534 
026442 
027499 


Feed 
Oe 


MAG 
RATIO 


12000 
06999 
02999 
0-999 
00999 
00999 
00998 
02998 
0.997 
0.997 
0«S995 
02994 
00992 
00989 
00986 
00981 
00975 
026966 
02955 
00941 
00922 
02899 
0«870 
0«834 
0«793 
Oe 745 
00693 
0636 
02578 
02520 
00463 
0410 
02360 
0+«314 
00274 
00237 
00205 
Oel77 
00153 
0e132 


29 


PHASE 
ANGLE 
DEGREE 


-220 
-204 
-2e7 
~3e2 
=3e7 
mgd 
= Sel 
=5.09) 
=609 
-820 
~9 04 
“1009 
“l2e7 
-1408 
~“17e2 
-20020 
—-2302 
-2629 
-31le2 
-36e1 
-41.8 
“46 e2 
5504 
-6306 
-72¢6 
-82.7 
~-93e7 
“10509 
“11962 
-133.8 
—-14909 
~l167e7 
~187e5 
-209e7 
—-23408 
—-2634% 
-29600 
—-33304 
—-37604% 
=“426e1 


Feed21- 
6001 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 =24 
1-000 -2.8 
00999 -323 
00999 -3 8 
00999 ~4e4 
02999 a W4 
02998 —-6 00 
0-998 -7.0 
0.997 —-Be2 
0.995 =-9e5 
00994 -llel 
02992 -129 
00989 -1520 
00985 -17¢5 
00980 =-2023 
Qe Ah2,. s2Zief 
02963 -27¢65 
Oe 95:0 =3200 
02933 -37el 
056910 4361 
020880 “4949 
06842 —57e7 
Oe792 ~=6604% 
00730 ~760e2 
00655 8607 
00569 =-9746 
00479 -10860 
00400 -116e8 
00347 -124e8 
00322 -136e3 
00293 -155640 
00236 -17608 
00186 =-18748 
Oe212 -201¢8 
00182 -224e1 
00441 -166e67 
0.090 9202 
OoN72 -79 09 
Oe072 —-3520 
02060 220 


Feed22- 
72092 
MAG PHASE 
RATIO ANGLE 
DEGREE 
12000 ~2el 
00999 —-204 
02999 =2¢«8 
0999 =3ie3 
00999 -328 
00999 -405 
00998 =5.0,2 
02998 -6el1 
00997 -7el 
00996 -803 
02994 -9 06 
00992 -lle2 
00e«990 =-13640 
0e986 -15062 
00982 -17026 
06975 2005 
00967 =2308 
02956 -2TeTl 
00941 =—32e1 
Oe921 =—37el 
026895 =-42¢8 
0.863 -49e3 
Oo821 —-56e5 
Oe771 =640e4 
Oeo713 =-72¢8 
0+650 —-81e4% 
0-589 “9003 
006535 9948 
0e489 ~llle2 
00437 -125e7 
06360 -141¢8 
00280 -15le7 
Oe282 -157el 
02334 -181¢0 
12049 -179el 
Ool77 =77e2 
00e028 =20606 
02060 —2050e8 
02061 -227¢e3 
Oe061 ©135¢5 
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TABLE DiS FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 


COMPOSITION TO NEGATIVE PULSES IN THE 
REFLUX FLOW RATE 


Reflux - Ref lux21- Ref lux24- 
GAIN 0099 15243 18.92 


FREQ MAG PHASE MAG PHASE MAG PHASE 
RAD/ RATIO ANGLE RATIO ANGLE RATIO ANGLE 


MIN DEGREE DEGREE DEGREE 
0e0020 1000 -1¢8 12000 -le7 12000 -203 
020023 00999 ~2el 0.999 -220 0.999 -2eT7 
Oe0027 00999 “205 06999 —-204 00999 -3e1 
0e0031 0-999 ~209 00999 “208 00999 —-306 
020036 04998 —-304 00999 —~3e2 02998 ~4ed 
020042 02998 -3.9 0.999 -3./7 0.998 —~4e9 
0e0049 02997 ~4e6 02998 “404 06997 -5e/ 
020057 02996 —“503 04997 —~5el 02996 —“6e7 
Oe0067 02995 “602 00997 ~5 09 02994. -7¢8 
020078 02993 ~7e2 040996 609 02992 =-Sel 
02-0091 00991 “B04 02994 “8el 00990 ~-1005 
020106 0-987 #98 026992 “904% 020986 =12e3 


060123 04983 =-1lle4 06989 =-100¢9 06981 -14e3 
000144 0646977 =130e2 00986 -12e7 00975 ~=1606 
020167 0«970 —-15e3 0-981 —-14e8 02966 -1903 
020195 06960 lle? 009TH H17e2 00954 72204 
000227 02946 =2004% 04965 1969 00938 ~26040 
000264 02929 =2365 00954 =—-23461 06917 =-3002 
000308 06908 2740 02938 2608 00889 =—34e9 
000358 04880 —~300e«8 Oe9l7 —-31.0 0+852 -40e2 
000417 00847 =35e0 006890 =3568 06605 ~4602 
000486 02808 =3966 06856 -4lel 06745 —-5207 
020566 O0e762 “4464 04812 -47.0 O«671 ~59e5 
020659 Oe711 4904 06758 =-5303 00584 =66040 
Oe0767 00656 ~5405 00695 ~59 09 00489 —“7le2 
000893 02598 =5906 00625 ~660e5 00401 “7302 
021040 06539 ~-64e7 006554 —72065 00342 -7lel 
001210 00482 —-690e8 00491 =-78e2 00324 ~6904% 
001409 O«427 —-74e7 00439 —-84.7 006319 —~74e5 
021640 04376 “7906 02388 “9345 00284 . ~85e5 
O0l1910 00329 =—-B84e4 00321 -104e3 04213 “9248 
022223 0287 “8902 Oe24l -112.4 OslTZ “8Bed 
002588 00249 —-9440 04189 -114.7 0-158 —9560 
003013 056216 -99el 06155 -124e6 0-103 9205 
003507 040186 =104e3 002104 —-13520 0-119 =-86e7 
024083 Oel6l “109-9 0ell6 -207+43 02-086 —-86e0 
004753 00138 -1160e0 04093 -1026e7 O+101 10243 
005534 O0ell9 -1220e7 0-102 -10005 04-093 9820 
006442 00el02 -1300e2 046093 -11565 02074 =-10263 
007499 02088 -138e¢6 0e097T —-1434¢7 02050 =-12006 
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TABLE D.10 FREQUENCY RESPONSE OF THE OVERHEAD PRODUCT 
COMPOSITION TO POSITIVE PULSES IN THE 
STEAM FLOW RATE 
Steam F Steam2l- Steam23- 
GAIN 0099 32032 25020 


FREQ MAG PHASE MAG PHASE MAG PHASE 
RAD/ RATIO ANGLE RATIO ANGLE RATIO ANGLE 


MIN DEGREE DEGREE DEGREE 
020020 12000 -2e3 12000 “168 1e000 =10e8 
02e0023 02999 -2e/7 00999 -“2062 00999 -2el 
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APPENDIX E 
FEEDFORWARD CONTROLLER MODEL DATA 


This section contains the frequency response repre- 
sentation of the feedforward controllers, in both graphical 
and tabular form, which were calculated from Equation 
(3.23) using the column open loop models determined pre- 


viously. 


Figure E.1] and Table E.1 illustrate the steady state 
open loop feedforward operating conditions exhibited about 
the reference steady state by the column. The feedforward 
Operating Conditions relate the reflux flow rate required 
to maintain the overhead product composition constant at 


various feed flow rates. 


Figures E.2 and E.3 illustrate the frequency response 
of the feedforward controllers calculated according to 
Equation (3.23) using the open loop frequency response mo- 
dels attributed to positive feed flow rate and negative 
reflux flow rate disturbances. Also plotted on these fi- 
gures is the frequency response of the feedforward control- 
ler calculated from the appropriate approximate models 
illustrated in Table 5.6. Figures E.4 and E.5 contain 
similar feedforward controllers calculated from the open 


loop frequency response models attributed to negative feed 
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flow rate and positive reflux flow rate disturbances. 


Similar feedforward controllers were calculated to 
maintain the overhead product composition constant by 
manipulating the steam flow rate correcting for feed 
flow rate disturbances. These controllers are given 


elsewhere (94). 
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CALCULATED FREQUENCY RESPONSE OF THE FEEDFORWARD 


CONTROLLER MANIPULATING THE REFLUX FLOW RATE 


FOR INCREASING FEED FLOW RATE DISTURBANCES 
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TABLE E.3 CALCULATED FREQUENCY RESPONSE OF THE FEEDFORWARD 
CONTROLLER MANIPULATING THE REFLUX FLOW RATE 
FOR DECREASING FEED FLOW RATE DISTURBANCES 
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APPENDIX F 
EXPERIMENTAL CONTROL STUDIES 


F.1 Introduction 


The transient response of the distillation column 
under various control schemes studied during this project 


is represented graphically in this section. 


Since the time period over which the initial steady 
state was evaluated is not constant for each run, the time 
at which the feed disturbance entered the column is noted 
on each figure by the arrow located on the time axis of 
the overhead product composition portion. These times have 


also been summarized in Table F.1. 


The experimental control responses have been presen- 
ted according: to the .types. of «control Ler; studied: 

a) analog feedback control 

b) digital feedback control 

c) gain feedforward control 

d) gain plus time delay feedforward control 

e) gain plus time lag feedforward control 

f) feedforward-feedback control 


g) comparison of digital feedback and feedforward- 
feedback control for a series of disturbances. 
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Each list has also been cross-referenced in Tables 


mem F.3, F.4 and F°5, Gaccordang to directions oi unc 


fed flow rate disturbances, including a listing-of .the 


omtroller parameters used. 


These figures illustrate the response of the major 


vwriables of interest, including: 


I 


i 


7] 


Tie 


et 


a) overhead product composition (controlled variable) 
b) reflux flow rate (manipulated variable) 
c) bottom product composition (output variable). 


atdition, certaifa fitiures Ure 10m F-clweteao dna! 49) 


sO present the response of other variables of interest, 


cluding: 


a) overhead product flow rate (output variable) 

b) bottom product flow rate (output variable) 

c) feed flow rate (disturbance variable) 
response of other process variables is afso avaziable 


ewhere (94). 
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F.2 Analog Feedback Control 


The transient response of the distillation column 
was determined under a conventional feedback control scheme 
using a Foxboro pneumatic Dynalog controller. The control- 
ler attempted to maintain the overhead product composition 
constant by manipulation of the set point of the reflux 
flow controller. The controller constants of the propor- 
tional plus integral controller were initially evaluated, 
using the Ziegler-Nichols criteria. The measured tran- 
sient responses of the controlled column using these con- 
troller constants were observed to exhibit only slightly 
damped oscillations, which were judged to be undesirable. 
The controller constants were tuned such that the dampening 
of the oscillations was increased. The controller con- 


stants are compared in Table F.6. 


TABEE F..6 ANALOG FEEDBACK CONTROLLER CONSTANTS 


Ziegler- 
Nichols Tuned 
Method Method 
Analog Representation 
PB -100% -200% 
Ty 5.0 9.0 
Digital Representation 
Ke oF 120 bt 0.5 
K O.. 1.07 0.059 
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The transient response of the column was measured 
for both increasing and decreasing step disturbances in 
the feed flow rate both away from and returning to the 
reference steady state conditions. The initial and 
final steady state operating conditions of the column 


are presented in Table F.7. 
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n OUTPUT _VARTASLE 
“| = 0-58 

a 

1«O06- 

O» 


Pele MANIPULATED VARIABLE 
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FIGURE F.1 ANALOG FEEDBACK CONTROL: (Z-N) 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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2-50 OUTPUT to ae: 
= 0-44 
407 
a 
Oe 
=(¢ 
Pod MANIPULATED VARIABLE 
SS = 2-04 
oy. TEST: FB-90A CONTROLLED VARIABLE 
Kxc= -0-8 SS =95-44 
tT) = 9.0 TAE= G74 
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FIGURE F.2 ANALOG FEEDBACK CONTROL: (TUNED) 


STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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OUTPUT ee 


MANIPULATED VARTABLE 
SS = 1-93 


——EE—— 
o7.oor TEST: FB-94A CONTROLLED VAR 
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FIGURE F.3 ANALOG FEEDBACK CONTROL (Z-N) 
STEP DECREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.4 ANALOG FEEDBACK CONTROL: (TUNED) 
STEP DECREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.6 ANALOG FEEDBACK CONTROL: (TUNED) 
STEP DECREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 
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FIGURE F.8 ANALOG FEEDBACK CONTROL: (TUNED) 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 
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F.3 Digital Feedback Control 


The transient response of the distillation column 
under a conventional feedback control scheme was imple- 
mented using the DDC program available on the IBM-1800 
digital computer. A standard proportional plus integral 
controller was used to maintain the overhead product com- 
position constant by manipulating the reflux flow rate. 
The controller constants were initially determined using 
the Ziegler-Nichols procedure. Again, the use of these 
controller constants resulted in slightly damped oscil- 
lations in the transient response. These constants were 
also tuned to increase the dampening of the oscillations. 


Table F.8 compares the values obtained. 
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The response of the column has been measured for 
both increasing and decreasing steps in the feed flow 
rate both away from and returning to the reference steady 
state. The initial and final steady state operating con- 


ditions are outlined in Table F.9. 


The effectiveness of the various feedforward control 
schemes, as measured by the integral of the absolute er- 
ror (IAE), has been compared with that obtained by feed- 
back control in order to establish the most effective 


control scheme of those studied. 
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FIGURE F.9 DIGITAL FEEDBACK CONTROL: (Z-N) 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.10 DIGITAL FEEDBACK CONTROL: (TUNED) 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.18 DIGITAL FEEDBACK CONTROL: (TUNED) 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 
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F.4 Gain Feedforward Control 


The transient response of the distillation column 
under gain feedforward control is presented for both 
increasing and decreasing step disturbances in the feed 
flow rate away from and returning to the reference steady 
state. The object of the control system was to maintain 
the overhead product composition constant by manipulating 
the reflux flow rate. Due to the nonlinear behavior ex- 
hibited by the column, it was necessary to determine the 
feedforward gain separately for the positive and negative 
feed flow rate disturbances. A feedforward gain of 0.5, 
as determined from Figure 6.2, was initially employed, but 
this value was subsequently tuned to eliminate the offsets 
which resulted. The tuned values were found to be -0.45 
for positive steps and -0.51 for negative steps in the feed 


flow rate. 


1 a 


amufo> mere en ts v0 8 es a 
dtod vot betnszeq et Tor 109 Br re 


\ 


best ond nt 2oanadwwd2tb ¢ ‘ 
ybssd2 sonsiste1 odd of Sansa brs 1 
nfsintem of 2 5W metey2. Torvtnoa ais ase 
patisluqtnsm yd tastenos nottteoqmoo ae 
-~x8 yworvedsd issatinon odd oF oul .5t84 ued ” 
ads antmretoab of yi1seeooen esw sf a ont | bette 
evitsesn bas svitteog sAt +o0t lh ah $e 5p b> Sea : 
«2.0 Yo nteo byswrotbest A 2sonsdyusetb 9367 Lose 


dns 
a 


oP a 


tud ,bsyolfqms yllststat esw ,.8$.d oyuptd mort bentmrsdab in | 
ets2tto odd otantmtis ot benus ylinsupsedue esw 9 Tsv if. 

é&.0- sd ot biuot s1ow esulsy bonus onT ‘pss tue04 Aotite 
best eds mt eqsat2 avtispen az0ot f2.0- bas eqos2. avisteag, 


332 


Pee 9] ay. 
rl Roc 92 
Oez7ulco- = 0 | 
Qaeees:2- Ll 
pif 20°2- ol 
g°€ FS°2- 3°0 
tL MOC’ y'2 
Onc EDC g°2 
Gays es eT 
9:5 B°0-- 1°2 
Saye eee l= SL 
ve ae) PZ 
1) a), (er 
Gcomosw-— 2° 1 
Pee 1 6 | 
Cece 0 le 6 Ss". 
face ote 80 
ice Bo 30- ec | 


0°H, HOPW, AD; 


aoueteg [eLuareW 


aul. 0050 > 71.296 
UL 0050" 71-296 


[ee OSes ae 0 189.6 
Lo. 0 ert 9.0296 
elee-0S5":02 6.0296 
be sO0C “£0 296 


GUS 0020-80596 
L OVO #tL‘96 


Olga 0 Vet ata 6 
CC eet Vet O0en 5 1629.6 
Col 005 60.596 
Vc 10S 0" 00496 
Oeecles Ol l= =S0296 
Oil ee-0 505, 907296 
Verte 10:9) *1..£30°496 
Cel 20S 0 20) 206 
Vi lemed ls 80 290 
6.02 10550 00 496 


sand3zno 


SUOLZLPUOD [eULY 
SUOLZLPUOD LPLZLUL - 


SiS? 
Sac? 


0°9P 
0°9P 
039P 
0797 


0°8t 
0°8P 


0°8P 
0°8t 
0°8P 
0°8P 
0°9t 
0°9 
0°8P 
0°8P 
0°8t 
0°8P 


X 


VlVd INILVYdSd0 NWN109 
>TOYLNOD GYVMYOIdSd4 NIV 


t 
ey LL 


sOsee V6 12-8)" 
Sums. cc cs 76 1 
S022 - £0. 2s Siac 
30° ¢ . VS th 262 
SOc 00° ch 8Y Zz 
a0 ¢c— Le e ¢c6 2 
SOS Ce 0 Co Ce a VO-| 
S02) VOWS ec 
00° 2: Ul) ae -16 2 
O032 06. ea Sic 
OOS o> 01 lea Gc 
LOStCa 0 Ont amc Ya 
SQ-¢ €8 1f-tc6 < 
G02 C220 2a eoy<c 
LOG, Oe See 
(O02 2e 76 Sale ace ec 
iWec— oL97t a Ga8c 
toe ¢ 6SC- 4 Lape 
S | = 
Sindu] 
OW 4S agi 


m— LL 


ALL eae LL 


me Li 


wt ee et LY Li 


Sa; J.44 


bemoda 
A IE 


9S~Jad 


él-Jd4 
Oia 
OCao 4a 
COS 43 
LO=J44 


7S9L 


| 


. 


1°30 1° 


a Si 


1°52 


“3°91 
13S 


0°20 
1*\e2 


¢8° 
6} #28°0 ae10 
#8°0 82°08 
t2°0. 38°00 0°20 


i7eh §701 
g2 §°0! 


H . 


- 


XB ( WT PC MEOH ) 


R ( LBAMIN ) 


XO ( WI PC MeQH ) 


333 


DUTPUT ae 
= O51 


TEST: FFC-01 CONTROLLED oS ee 
Kee=-0. 43 SS =25-03 
TAE=10°OF 


@) e4 43 72 = 120 


Tike (MIN) 


FIGURE F.19 GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.22 GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.23 GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
EFFECT OF SIMULTANEOUS DISTURBANCE IN 
THE FEED TEMPERATURE (AT. =-10) 
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FIGURE F.24 GAIN FEEDFORWARD CONTROL: 


STEP DECREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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GAIN FEEDFORWARD CONTROL: 
STEP DECREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 
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FIGURE F.26 GAIN FEEDFORWARD CONTROL: 


STEP DECREASE IN THE FEED FLOW RATE 
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GAIN FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 
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F.5 Gain plus Time Delay Feedforward Control 


The transient response of the column has been deter- 
mined under a gain plus time delay feedforward control for 
both increasing and decreasing step disturbances in the 
feed flow rate away from and returning to the reference 
steady state. The control objective was to maintain the 
overhead product composition constant by manipulating the 
reflux flow rate. The feedforward gains, determined pre- 
viously, were used in conjunction with various values of 
the time delay chosen from the range suggested in Tables 
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FIGURE F.28 GAIN PLUS TIME DELAY FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
AWAY FROM THE REFERENCE STEADY STATE 
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FIGURE F.36 GAIN PLUS TIME DELAY FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 
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F.6 Gain plus Time Lag Feedforward Control 


The transient response data of the column under 
gain plus time lag feedforward control was measured for 
both increasing and decreasing step disturbances in feed 
flow rate away from and returning to the reference steady 
state. The feedforward gain used was that determined 
previously during the gain feedforward control tests. 
The values of time constant of the time lag function 


were chosen from the range illustrated in Table 7.4. 
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FIGURE F.45 GAIN PLUS TIME LAG FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 


SLBATAAV CSTAUATIMAM 
IS-5 = <2 


Ne Warf vanh yo 


¥ 
i 


plane 


AJZATRAVY OUGLUDATVOO 3-974 :Teat 


Of-ece= 22 fe O-*544 
ON+ ESAT i — ay 


| 
ost a 


BS oO 


St Se 


Wh) Sar 


:JONTMOD GRAWROIGIZA OAL AMIT 2019 WIAD 2.7 3nUaTA 
ITAA WO C332 BHT WI FASAML 9aT2 
3TAT2 YGAIT2 3QW3AZISA ANT OT BWIMAUTAA 


7 


XB ( WT PC MeQH ) 


R ( LEAMIN ) 


) 


XO ( WT PC Me 


366 


ee , QUTPUT VARIAELE 
| SS = 0-00 
1°75 
1°00-- 
O- 
-—<Os 
Le MANIPULATED Vi RIAELE 
S595 = eel 


ce 
Mn AU nl Any Vd NA, 
4. 
ee TEST: FFC-59 CONTRE_LLED VARTASLE 
Kep=-0.51 SS =55-03 
tT = 6.6 ribsahiigm 
> 
O5- 
S5-00- t ’ 
O c4 495 Jeo Ss 420 


Tite (MIN) 


FIGURE F.46 GAIN PLUS TIME LAG FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 


SSSATRAY CSTAJUAIMAM 
iS-S = & 


hia ania J aa a] 
ae ee 
BRATBAY OUSELLORTVOO @@-2472 :T23aT 
COsacs ad fa 10-8594 
Of+£f=3AT 


ic 
Wi) = SIT 


b 


ITAA WOST OFF IHT AL A2ARAOWT GATS 


| | 2 
:JONTHOD ORAWROWIIIA DAY IMIT ZWI9 ae 30.4 38U9T4 
ATAT2 YOASTE JOMIRIIIA ANT OT BHIMAUTIW 


XB ( WI PC MECH ) 


R ( LEMIN ) 


XD ( WT PC MeQH ) 


367 


OUTPUT VA ea 
SS = 


MANIPULATED VARIASLE 
SS = e-c0 


TESES 


Kee=-O. 
tT = 4, 


FFC-65 eae VARTAGLE F 
5] =o °4{3 
5 TE = 233 


FIGURE F.47 


24 43 72 SS 420 
TIKe (MIN) 


GAIN PLUS TIME LAG FEEDFORWARD CONTROL: 
STEP INCREASE IN THE FEED FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 


LEATAAY lias 
05-5 = 


A aa a 
JY AW YAK yoni Mahan | 
————-_——} —_ —}— 
@RAIRAY CLOTTED 20-298 <feet aan 
Efsce= a2 id Ontagd . 


EES =SA 
“> AT aa — 


EN 
Wi) QMIT 


> JOATHOD QRAWAOVGIII BAS BMIT 2019 WIAD \O.9 anuatd 
ITAA WOT G3F9 JHT WI J2ARARIMI 93T2 . 
JTAT2 YOAIT? JIMIAIIIA FHT OT ONIMAUTIA 


XB ( WT PC MECH ) 


R ( LBAMIN ) 


) 


7 


f 
eas 


XO ( WT PC! 


368 


2.5 CUTPUT = VARTAELE 
SS = O-CO 
41°75 
1°00 
0-2 
-O- th tt 
>. MANIPULATED D_VARIAELE 
2 
Cet 
ee 
_ iy van tian VMAS AAR othe tla Mal 
4. 
o.oo. TEST: FEC- 57 CONTROLLED __ VARIABLE 
haces TAEeiOsdt 
15 = e 


O 24 


FIGURE F.48 


48 72 oS 420 
Tike (MIN) 


GAIN PLUS TIME LAG FEEDFORWARD CONTROL: 
SvGPeINCREASE. IN THE? FEED: FLOW RATE 
RETURNING TO THE REFERENCE STEADY STATE 


SUSATRAY OST AUALMAM 
ISS = 22 


WV deans Nay giiiat Ap AAW venti AM | 


$—— 


SLBALSAV _CELLOATMED {2-399 <T23T sog.ep oo 
= ea fa 0-8, 54 _ 
i. OESaAT ret. a 
4 a a eee 
a ve siti _ L002 
OS 2s sf Qh AS . @- . 
(MER = BAIT 


eS tn 
:LONTHOS GRAWHOIGISI DAL AMIT 2449 “ies 86.7 3AValT 
3TAR WO G339 FHT WI F2AFROMI GAT 
ITAT2 VOAST? JOMANTISA GUT OF SHLRAUT aA 


F.7 Feedforward-Feedback Control 


The transient response of the column operating 
under combined feedforward-feedback control was deter- 
mined for both increasing and decreasing step distur- 
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utilized included both gain and gain plus time lag models. 
The feedforward controllers used were those determined 
from previous tests which produced the most effective 
control for the individual control schemes. The time 
delay feedforward controller was not implemented since a 
pure time delay cannot be implemented easily without the 
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F.8 Series of Disturbances 


The effectiveness of a conventional feedback con- 
troller and a gain plus time lag feedforward-feedback 
controller were compared for a series of increasing and 
decreasing step disturbances in the feed flow rate about 
the reference steady state. The feedforward parameters 
were chosen as the average of the best values determined 
previously for each disturbance direction. The tuned 
feedback controller parameters were also used in each 


controller. 
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APPENDIX G 
PROGRAM DESCRIPTIONS 


G.1 Introduction 


Several programs were written during the course of 
Cols, project. Inese programs fall into’ the following ‘two 
categories: 
i) user interfaces to DDC 
ii) data reduction. 


Listings of these programs are given elsewhere (94). 
G.2 User Interfaces to DDC 
a) ONOFF 


The execution of this program, queued from the key- 
board, allows a number of options to be performed. These 
options include: 

- turn specified variable loop records on 

- turn specified variable loop records off 

- turn specified data accumulation ring buffers on 

- turn specified data accumulation ring buffers off 

- initiate buffering to disk 

- start gas chromatograph cycle 


- stop gas chromatograph. 
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This program was most useful since it allowed the 
bulk turning on and off of loop records, rather than the 
tedious procedure of initiating them individually from the 


process operator's consol, as is normal procedure. 
b) DATAC/DASS/BALNC 


This set of programs, initiated from the keyboard, per- 
forms the steady state material balance calculations. These 
programs were executed before and after each dynamic or 
control test. DASS initializes parameters required by DATAC 
and BALNC. DASS also issues the command to queue program 
DATAC at a regular time interval, such as every 30 seconds, 


Tor a specified length of time," te 15 minutes. 


DATAC is queued at a regular interval, as defined to 
the system, to read the analog input points of all the column 
variables, converts the measured values to engineering units 
and calculates a running sum and sum of squares. These va- 
lues for each column variable are stored on file in order to 
be saved for the next time interval. When sufficient data 
has been accumulated, DATAC queues program BALNC and then 


cancels the repetitive queuing of itself. 


BALNC initially calculates the average and standard 


deviation of all the measured column variables from the 
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running sums and sums of squares which were stored on the 
file. The program then performs the standard material ba- 
lance calculations about the column and issues a report 


shown in Table A.4. 
G)enGTDAL 


This program provides the interface to the system loop 
record buffers on the disk. The program will sort the data 
contained in the system file based on the specified loop 
record number. This data is then stored based on this 
loop record number in a user file. When all the data for 
up to 10 specified loop records has been found, the data is 


output .on both .the printer and cards for future reference. 
G.3 Data Reduction Programs 
a) PTAP 


This program, as described by Wildman (135), reduces 
the transient response produced from a pulse test on the 
column into a frequency response diagram. The transient 
Seance of the desired output and disturbance variables is 
read into the program PTAP from cards produced by GTDAT at 
the end of a pulse test. The program execution and subse- 


guent graphical output is carried out according’ to the pro- 
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cedures outlined by Wildman (135). 
b) GEN 


This program produces a table of frequency response 


values for a given complex transfer function. 
c) LEVY 


This program fits transfer functions to frequency 
response data using the method of Levy (63). The frequency 
response data are initially read in from cards produced 
during the execution of PTAP or FFCl. The program will 
now fit the frequency response data to a specified trans- 
fer function model. A series of models may be fit in 
sequence to the same frequency response data, in order to 


compare their degree of fit and appropriateness. 
dy), FFG 


This program calculates the frequency response of a 
feedforward controller based on the open loop frequency 
response values of the manipulated and disturbed values. 
The open loop frequency response values used are those 
punched on cards during the pulse testing analysis using 


PTAP. 
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e) MLR 


MLR is a standard multiple linear regression program 
which was used to correlate the various properties of the 
methanol and water system, as outlined in Appendix A. 
Various on-line programs, such as BALNC and GTDAT, used 
these correlations to generate the properties required 


in the calculations. 
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